Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,066,671 C→A 6.7% G175V (GGA→GTA)  yihP ← putative 2,3‑dihydroxypropane‑1‑sulphonate exporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,066,6710CA 6.7% 293.1 / 5.5 86G175V (GGA→GTA) yihPputative 2,3‑dihydroxypropane‑1‑sulphonate exporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base C (41/39);  new base A (3/3);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGTTTCGACATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTAAAG  >  NC_000913/4066533‑4066813
                                                                                                                                          |                                                                                                                                              
tGGGTTTCGACATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGca                                                                                                                                      >  4:354711/1‑149 (MQ=255)
tGGGTTTCGACATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGca                                                                                                                                      >  6:370382/1‑149 (MQ=255)
tGGGTTTCGACATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGca                                                                                                                                      >  6:353948/1‑149 (MQ=255)
  ggTTTCGACATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGca                                                                                                                                      <  3:354711/147‑1 (MQ=255)
        gACATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCAcagcag                                                                                                                                >  7:250641/1‑147 (MQ=255)
         aCATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCAcagcagcag                                                                                                                             <  5:44508/149‑1 (MQ=255)
              aCGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGccc                                                                                                                         <  5:353948/148‑1 (MQ=255)
              aCGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCa                                                                                                                        <  7:119501/149‑1 (MQ=255)
               cGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCTAAGATATAGCCAAGTTGCTGATTACCTTCGATATGATTCATAACTGCACTAATACCACCGTGCACAGCAGCAGGCCCAGc                                                                                                                      <  5:164694/149‑1 (MQ=255)
                gCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGc                                                                                                                      <  2:306570/149‑1 (MQ=255)
                  tctTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCg                                                                                                                     >  7:151479/1‑148 (MQ=255)
                                   ccacctCCACATAAACAGCATTACAAACAGTGCAAACAGAGCGTCTGCTACGATATCGCCAACTTGCTTATTACCCTCGATTCGATTCCTAACTCGCTCTAATCCCAGGTTTCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGcc                                                                                                   <  4:252292/144‑1 (MQ=255)
                                    cAGATCCACATAAACAGCAGGCCAAACAGGGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCa                                                                                                  <  6:124815/149‑1 (MQ=255)
                                          cacaTAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATgccgc                                                                                            >  5:190049/1‑149 (MQ=255)
                                           acaTAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCAAGCGTAGCGCCGCCCTGACGCCATgccgc                                                                                            >  8:40627/1‑148 (MQ=255)
                                             aTAAACAGTAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATACCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATgccgc                                                                                            >  4:302741/1‑146 (MQ=255)
                                             aTAAACAGTAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATACCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATgccgc                                                                                            <  3:302741/146‑1 (MQ=255)
                                              tAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGt                                                                                        >  3:244867/1‑149 (MQ=255)
                                                      agGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGcgcg                                                                                >  6:216139/1‑149 (MQ=255)
                                                       ggCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGcgt                                                                               >  2:210548/1‑149 (MQ=255)
                                                               aGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTggg                                                                       <  4:244867/149‑1 (MQ=255)
                                                               aGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTggg                                                                       >  7:62562/1‑149 (MQ=255)
                                                                 tGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCAcagcagca                                                                                                                              <  2:262256/92‑1 (MQ=255)
                                                                 tGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCAcagcagca                                                                                                                              >  1:262256/1‑92 (MQ=255)
                                                                 tGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATgc                                                                                               >  4:33421/1‑123 (MQ=255)
                                                                 tGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATgc                                                                                               <  3:33421/123‑1 (MQ=255)
                                                                    aaaCAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGtttt                                                                  >  3:338230/1‑149 (MQ=255)
                                                                      aCAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTg                                                                <  3:79895/149‑1 (MQ=255)
                                                                      aCAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTg                                                                <  7:226483/149‑1 (MQ=255)
                                                                        aGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGt                                                              <  3:180549/149‑1 (MQ=255)
                                                                        aGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGt                                                              <  5:97889/149‑1 (MQ=255)
                                                                        aGCGTGGCGGCGAAGATATAGCCAAGTGGCTGATTACCTTCGATAAGATTCATCACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGt                                                              <  3:185461/149‑1 (MQ=255)
                                                                         gCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGcgttcgt                                                                           <  7:172078/135‑1 (MQ=255)
                                                                         gCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGcgttcgt                                                                           >  8:172078/1‑135 (MQ=255)
                                                                         gCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTgg                                                               <  8:254178/147‑1 (MQ=255)
                                                                         gCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTgg                                                               >  7:254178/1‑147 (MQ=255)
                                                                           gTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGAt                                                           <  1:210548/149‑1 (MQ=255)
                                                                               aggcgAAGATATAGCCAAGTTGCTGAGTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGAAGCGCCGCCATGACGCCATGACTCCAGTGAGGCGCGTTCTTGGCGGTTTTTGGTGATATCg                                                       >  6:295710/2‑149 (MQ=255)
                                                                                   gAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGaa                                                   <  6:190049/149‑1 (MQ=255)
                                                                                        tataGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCAtt                                              <  8:151479/149‑1 (MQ=255)
                                                                                         ataGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGTAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGACCCATTg                                             >  3:324448/1‑149 (MQ=255)
                                                                                            gCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTgcgc                                          >  2:367554/1‑149 (MQ=255)
                                                                                             ccAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTg                                                                         >  2:374527/1‑117 (MQ=255)
                                                                                             ccAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTg                                                                         <  1:374527/117‑1 (MQ=255)
                                                                                             ccAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGa                                                            <  3:156621/130‑1 (MQ=255)
                                                                                             ccAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGa                                                            >  4:156621/1‑130 (MQ=255)
                                                                                                aGTTTCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCTCCCTGACGCCATGCCGCCAGTGATCCGCGTTCGTTGGGGTTTTTGGTGATAGCGTGAACCATTGCGCCGTa                                      <  5:295710/149‑1 (MQ=255)
                                                                                                aGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTa                                      <  5:216139/149‑1 (MQ=255)
                                                                                                   tGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAg                                     >  8:283193/1‑147 (MQ=255)
                                                                                                       gATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAg                               >  1:185002/1‑149 (MQ=255)
                                                                                                         ttACCTTCGATAAGATTCATAACTGGCACGAATACCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGa                                   <  3:106305/143‑1 (MQ=255)
                                                                                                         ttACCTTCGATAAGATTCATAACTGGCACGAATACCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGa                                   >  4:106305/1‑143 (MQ=255)
                                                                                                             cTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCgg                                                      <  8:175811/120‑1 (MQ=255)
                                                                                                             cTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCgg                                                      >  7:175811/1‑120 (MQ=255)
                                                                                                             cTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTtcat                          >  4:317001/1‑148 (MQ=255)
                                                                                                              ttCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGCGCAGTtcatca                        >  6:252378/1‑149 (MQ=255)
                                                                                                              ttCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTtcatca                        >  1:292026/1‑149 (MQ=255)
                                                                                                                cGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGCGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTcatcat                       <  5:370382/148‑1 (MQ=255)
                                                                                                                 gATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCAt                                                                                                 <  5:193025/73‑1 (MQ=255)
                                                                                                                 gATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCAt                                                                                                 >  6:193025/1‑73 (MQ=255)
                                                                                                                  aTAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCAtt                                              <  2:226005/123‑1 (MQ=255)
                                                                                                                  aTAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCAtt                                              >  1:226005/1‑123 (MQ=255)
                                                                                                                  aTAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCt                    >  4:181858/1‑149 (MQ=255)
                                                                                                                   tAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCt                    >  1:253539/1‑148 (MQ=255)
                                                                                                                      gATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGgc                                                                                   <  4:130427/82‑1 (MQ=255)
                                                                                                                      gATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGgc                                                                                   >  3:130427/1‑82 (MQ=255)
                                                                                                                      gATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGc                                                        <  8:25181/109‑1 (MQ=255)
                                                                                                                      gATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGaa                 >  7:260455/1‑148 (MQ=255)
                                                                                                                        ttCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGaaaaa              <  3:181858/149‑1 (MQ=255)
                                                                                                                           aTAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGtttt                                                                  <  8:300033/94‑1 (MQ=255)
                                                                                                                           aTAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGtttt                                                                  >  7:300033/1‑94 (MQ=255)
                                                                                                                           aTAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAg           <  4:338230/149‑1 (MQ=255)
                                                                                                                            tAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGt          >  3:86514/1‑149 (MQ=255)
                                                                                                                              aCTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGACCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTcc        <  4:86514/149‑1 (MQ=255)
                                                                                                                               cTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGcgttcgtt                                                                          >  3:164008/1‑82 (MQ=255)
                                                                                                                               cTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGcgttcgtt                                                                          <  4:164008/82‑1 (MQ=255)
                                                                                                                                tGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCg       >  4:179535/1‑148 (MQ=255)
                                                                                                                                   cACGAATCGCTCCGTGCACTGCAGCAGGCCCAGCGTCGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATCCTGAAAAACAGTCCGTaaa   <  6:86513/149‑1 (MQ=255)
                                                                                                                                   cACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTaaa   <  2:253539/149‑1 (MQ=255)
                                                                                                                                   cACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTaaa   <  8:260455/149‑1 (MQ=255)
                                                                                                                                   cACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTa     >  4:196115/1‑147 (MQ=255)
                                                                                                                                    aCGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGt                              <  3:585/121‑1 (MQ=255)
                                                                                                                                    aCGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGt                              >  4:585/1‑121 (MQ=255)
                                                                                                                                    aCGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTAAAg  >  7:194086/1‑149 (MQ=255)
                                                                                                                                    aCGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTAAAg  >  1:183246/1‑149 (MQ=255)
                                                                                                                                    aCGAATACCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTAAAg  >  5:43483/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                              
TGGGTTTCGACATAACGCTCTTTCACGCCCGAGTAGCAGATCCACATAAACAGCAGGCCAAACAGTGAAAACAGCGTGGCGGCGAAGATATAGCCAAGTTGCTGATTACCTTCGATAAGATTCATAACTGGCACGAATCCCACCGTGCACAGCAGCAGGCCCAGCGTAGCGCCGCCCTGACGCCATGCCGCCAGTGAGGCGCGTTCGTTGGGGTTTTTGGTGATAGCGGGAACCATTGCGCCGTAGGAGCAGTTCATCATGCTGAAAAACAGTCCGTAAAG  >  NC_000913/4066533‑4066813

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: