Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,155,255 C→T 5.2% R414C (CGC→TGC)  mdtA → multidrug efflux system, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,155,2550CT 5.2% 532.5 / 13.2 155R414C (CGC→TGC) mdtAmultidrug efflux system, subunit A
Reads supporting (aligned to +/- strand):  ref base C (73/74);  new base T (4/4);  total (77/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 155 total aligned reads displayed.

AGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGC  >  NC_000913/2155113‑2155401
                                                                                                                                              |                                                                                                                                                  
  gTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGa                                                                                                                                            >  2:45797/1‑149 (MQ=255)
  gTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGa                                                                                                                                            <  2:658943/149‑1 (MQ=255)
   tATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCGTGCTCCtttt                                                                                                                                           <  1:229643/149‑4 (MQ=255)
     tttCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGc                                                                                                                                         >  4:699063/1‑149 (MQ=255)
       tCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAg                                                                                                                                       <  3:627506/149‑1 (MQ=255)
        cTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAgg                                                                                                                                      <  2:329507/149‑1 (MQ=255)
            gggCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGtt                                                                                                                                  >  4:1598163/1‑149 (MQ=255)
            gggCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGtt                                                                                                                                  >  4:1182782/1‑149 (MQ=255)
            gggCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGCAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGtt                                                                                                                                  <  3:1100668/149‑1 (MQ=255)
                gATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTAccc                                                                                                                              >  3:1250409/1‑149 (MQ=255)
                   cgcgTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGAAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCg                                                                                                                           <  2:222406/149‑1 (MQ=255)
                     cgTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGag                                                                                                                         <  3:84397/149‑1 (MQ=255)
                      gtggtgACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGagc                                                                                                                        <  3:1266075/149‑1 (MQ=255)
                         gtgACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGCAGCCCAAAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGtgcagc                                                                                                                     >  1:1861068/1‑149 (MQ=255)
                          tgACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAgcagca                                                                                                                    >  1:810742/1‑149 (MQ=255)
                             cAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAg                                                                                                                 <  3:699063/149‑1 (MQ=255)
                              aGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAgg                                                                                                                >  1:671370/1‑149 (MQ=255)
                               gACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCGTGCTCCttt                                                                                                                                            >  3:1509321/1‑118 (MQ=255)
                                aCGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggcg                                                                                                              >  1:1359677/1‑149 (MQ=255)
                                 cGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAggcgg                                                                                                             <  1:49557/149‑1 (MQ=255)
                                     aTTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCg                                                                                                         >  2:1291794/1‑149 (MQ=255)
                                     aTTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCg                                                                                                         >  2:422835/1‑149 (MQ=255)
                                        gATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCGTGCTCCttt                                                                                                                                            >  4:92234/1‑109 (MQ=255)
                                        gATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCGTGCTCCttt                                                                                                                                            <  3:92234/111‑3 (MQ=255)
                                        gATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTcg                                                                                                      >  2:1705133/1‑149 (MQ=255)
                                         aTCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCTGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTcgc                                                                                                     >  1:123155/1‑149 (MQ=255)
                                         aTCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGgt                                                                                                                                     <  1:295171/117‑1 (MQ=255)
                                         aTCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTcgc                                                                                                     >  1:1245647/1‑149 (MQ=255)
                                           cGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCgcgc                                                                                                   <  1:422835/149‑1 (MQ=255)
                                           cGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCgcgc                                                                                                   <  3:1684949/149‑1 (MQ=255)
                                            gCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGcc                                                                                                  >  4:1264482/1‑149 (MQ=255)
                                             ccTGACCGAAGGGGCTAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCt                                                                                                 >  4:1402013/1‑149 (MQ=255)
                                             ccTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCt                                                                                                 >  3:723928/1‑149 (MQ=255)
                                             ccTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCt                                                                                                 >  2:251539/1‑149 (MQ=255)
                                                gACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGtt                                                                                              >  2:1240367/1‑149 (MQ=255)
                                                   cGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttat                                                                                           <  1:1396762/149‑1 (MQ=255)
                                                    gAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatt                                                                                          >  2:1254703/1‑149 (MQ=255)
                                                    gAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatt                                                                                          <  4:723928/149‑1 (MQ=255)
                                                        gggCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGc                                                                                      <  3:1569778/149‑1 (MQ=255)
                                                        gggCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGc                                                                                      <  1:491296/149‑1 (MQ=255)
                                                          gCGACAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGt                                                                                    >  4:1773825/1‑149 (MQ=255)
                                                           cGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTc                                                                                   <  3:1598163/149‑1 (MQ=255)
                                                           cGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTc                                                                                   <  3:1182782/149‑1 (MQ=255)
                                                            gAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGagca                                                                                                                       >  2:1552218/1‑112 (MQ=255)
                                                            gAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGagca                                                                                                                       <  1:1552218/112‑1 (MQ=255)
                                                                gtggaagtggTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCtgtg                                                                              <  4:1131165/149‑1 (MQ=255)
                                                                gtggaagtggTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCtgtg                                                                              >  3:1344715/1‑149 (MQ=255)
                                                                    aagtggTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTCTATTATGCGTCCTGTGGcca                                                                          >  3:187789/1‑149 (MQ=255)
                                                                     agtggTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGccac                                                                         <  1:1705133/149‑1 (MQ=255)
                                                                           ggAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAg                                                                                                                 <  1:1050087/103‑1 (MQ=255)
                                                                           ggAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAg                                                                                                                 >  2:1050087/1‑103 (MQ=255)
                                                                           ggAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACgct                                                                   <  1:748904/149‑1 (MQ=255)
                                                                            gAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatta                                                                                         <  4:1675043/126‑1 (MQ=255)
                                                                            gAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatta                                                                                         <  3:376229/126‑1 (MQ=255)
                                                                            gAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatta                                                                                         >  4:376229/1‑126 (MQ=255)
                                                                                 ccAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACgc                                                                    >  2:139539/1‑142 (MQ=255)
                                                                                  cAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGAtg                                                            <  1:1254703/149‑1 (MQ=255)
                                                                                    gagCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGAtggt                                                          <  3:1402013/149‑1 (MQ=255)
                                                                                     agCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGAtggtg                                                         <  4:1344715/149‑1 (MQ=255)
                                                                                         ccACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGa                                                     <  1:811782/149‑1 (MQ=255)
                                                                                          cACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCtgt                                                                               >  2:336581/1‑122 (MQ=255)
                                                                                          cACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACgctg                                                                  >  3:23849/1‑135 (MQ=255)
                                                                                          cACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACgctg                                                                  <  4:23849/135‑1 (MQ=255)
                                                                                          cACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGAt                                                    <  2:627593/149‑1 (MQ=255)
                                                                                          cACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGAt                                                    >  3:1547292/1‑149 (MQ=255)
                                                                                           actactCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATc                                                   >  4:1119822/1‑149 (MQ=255)
                                                                                           actactCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATc                                                   <  3:1063480/149‑1 (MQ=255)
                                                                                           actactCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGGTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATc                                                   <  2:810742/149‑1 (MQ=255)
                                                                                               ctCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCt                                                  <  1:625218/146‑1 (MQ=255)
                                                                                               ctCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTAc                                               <  4:187789/149‑1 (MQ=255)
                                                                                                tCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTtattat                                                                                        <  4:504275/107‑1 (MQ=255)
                                                                                                tCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGAtg                                                            <  2:1389191/135‑1 (MQ=255)
                                                                                                tCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGAtg                                                            >  1:1389191/1‑135 (MQ=255)
                                                                                                tCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTa                                                >  2:504369/1‑147 (MQ=255)
                                                                                                  cGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCg                                            >  2:1602544/1‑149 (MQ=255)
                                                                                                        agaAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGa                                      >  3:664744/1‑149 (MQ=255)
                                                                                                        agaAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGa                                      >  2:78680/1‑149 (MQ=255)
                                                                                                        agaAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGa                                      <  1:251539/149‑1 (MQ=255)
                                                                                                         gaAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCtgtg                                                                              <  2:1164184/108‑1 (MQ=255)
                                                                                                          aaaGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCt                                                                                                 >  4:748338/1‑88 (MQ=255)
                                                                                                          aaaGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCt                                                                                                 <  3:748338/88‑1 (MQ=255)
                                                                                                           gtaacaccaGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCg                                         <  1:1083965/139‑1 (MQ=255)
                                                                                                           gtaacaccaGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCg                                         >  2:1083965/5‑143 (MQ=255)
                                                                                                             gCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTa                                                >  3:968601/1‑134 (MQ=255)
                                                                                                             gCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATc                                 <  1:1047795/149‑1 (MQ=255)
                                                                                                             gCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATc                                 <  3:1773825/149‑1 (MQ=255)
                                                                                                                accaGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGt                              <  4:1547292/149‑1 (MQ=255)
                                                                                                                accaGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGt                              >  2:1438677/1‑149 (MQ=255)
                                                                                                                  caGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTa                            >  1:1881295/1‑149 (MQ=255)
                                                                                                                     ccGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGa                                      <  4:753936/136‑1 (MQ=255)
                                                                                                                     ccGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGa                                      >  3:753936/1‑136 (MQ=255)
                                                                                                                       gcgAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTa                            >  4:1559707/1‑144 (MQ=255)
                                                                                                                       gcgAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATcgcg                       >  1:113514/1‑149 (MQ=255)
                                                                                                                             aCGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGc                 <  1:1602544/149‑1 (MQ=255)
                                                                                                                              cgcgAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTttatta                                                                                         >  2:1243320/1‑76 (MQ=255)
                                                                                                                                   aaaaaGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGttt           <  1:1240367/149‑1 (MQ=255)
                                                                                                                                    atcaGGAGCGTGCTCCTGATGAAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGc                                                       <  1:347996/101‑1 (MQ=255)
                                                                                                                                    aaaaGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTc          >  1:508122/1‑149 (MQ=255)
                                                                                                                                      aaGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGc                                                                            >  2:1547105/1‑81 (MQ=255)
                                                                                                                                      aaGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGc                                                                            <  1:1547105/81‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
AGGTATTTCTGCGGGCGATCGCGTGGTGACAGACGGCATTGATCGCCTGACCGAAGGGGCGAAAGTGGAAGTGGTGGAAGCCCAGAGCGCCACTACTCCGGAAGAGAAAGCCACCAGCCGCGAATACGCGAAAAAAGGAGCACGCTCCTGATGCAGGTGTTACCCCCGAGCAGCACAGGCGGCCCGTCGCGCCTGTTTATTATGCGTCCTGTGGCCACCACGCTGCTGATGGTGGCGATCTTACTCGCCGGGATTATCGGTTATCGCGCCCTGCCCGTTTCGGCGCTGC  >  NC_000913/2155113‑2155401

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