Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 19.5% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC19.5% 401.9 / 55.3 245intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (70/123);  new base C (3/49);  total (73/172)
Fisher's exact test for biased strand distribution p-value = 5.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164432‑164727
                                                                                                                                                     |                                                                                                                                                  
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCag                                                                                                                                                   >  1:6249569/1‑149 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                   >  1:6651900/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                   >  1:2884728/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                 <  2:7067781/150‑1 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACCTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAgcgcg                                                                                                                                            <  2:6828957/150‑1 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                         >  2:3799305/1‑150 (MQ=255)
                 agaTGAAAGGGCGATATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTATTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGtgaac                                                                                                                                   >  1:6447298/1‑146 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccgc                                                                                                                               >  1:4372536/1‑148 (MQ=255)
                      aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccaca                                                                                                                              >  2:5392020/1‑150 (MQ=255)
                      aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgca                                                                                                                              >  2:8180735/1‑150 (MQ=255)
                        ggggCGTTACCCACAACCGGTCGGCCCGGACGACCCGCCAAATCCTCCACCGACGCTAGGGACCAAAAATTATTCGTACGGTTGTTTTTCCCCATCGCCGTGGGGGTGCGGGCGGGGTTGAGGGCCTCAGCGCGCACGTtcaccctcacc                                                                                                                            <  2:4451285/149‑1 (MQ=255)
                        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcccccgcacc                                                                                                                            >  1:10184034/1‑150 (MQ=255)
                            cGTTAGCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                        <  2:1483505/150‑1 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGGTTTCTCCCTCTCCCAGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCAAAcccaca                                                                                                                  >  2:8550684/1‑147 (MQ=37)
                                   ccAAACATGTCTGGTCGGACGACCCGGAAAATACTGCACCGACGGGACGTACGAAAATGTTTTGGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                 <  2:5720308/150‑1 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                 >  2:6244560/1‑150 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAcccacac                                                                                                                 >  2:6351426/1‑146 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCCCCCCAAcccacac                                                                                                                 >  1:2188196/1‑146 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCCAAcacacac                                                                                                                 >  1:4724308/1‑144 (MQ=37)
                                    cAACCATGTCGGGCCGGACGACCCGGCAAATAGTGCACCGCCGCGAGGGACGAAAATGTTTTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGTAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                <  2:9151552/150‑1 (MQ=255)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCACACCCAAAACCGCAccc                                                                                                               >  1:2459162/1‑150 (MQ=37)
                                       acctGTCTGGCCGGCCGCCCCGGCAAATCCTGCCCCGACGCGACGGACGAAAAAGTATTCTAATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                             <  2:3212728/147‑1 (MQ=37)
                                          tGTCTGGCCGTACGACCCGGAAAATAGTGCACGGCCGCTAGGTACTAAAATTTTTTCGTATGGTTGTTTTTCTCCCTATCCCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                          <  2:4651175/150‑1 (MQ=255)
                                          tGTCTGGCCGGACGACCCGTTAAATACTGCACCGACGCGACGTACTAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCATCACCCTAACCCTCTCCCTCAAg                                                                                                          <  2:10255683/150‑1 (MQ=37)
                                           gcggggCCGGACGACCCGGAAAATACTGCACCGCCGCGAGGGACGAAAATGTTTTCGTATGGTTGTTTTTCGCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                         <  2:5207047/146‑1 (MQ=255)
                                           gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCAAACCCGCTGCCGCCAgg                                                                                                         >  2:9889288/1‑150 (MQ=38)
                                               ggCCGGACGCCCCGGCAAATCCTGCACCGACGCGCGGGACGACACAGTATTCGTACGGTTGGTTTTCTCCCTCCCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  2:601650/150‑1 (MQ=255)
                                               ggCCGGACGACCCGGCAAATCCTGCACCGCCGCGCCGGACGAAAAAGTATTCGTAAGGTGGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  2:9309788/150‑1 (MQ=39)
                                               ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCAAAACCACACCCACCAGGgaga                                                                                                     >  1:9104064/1‑150 (MQ=37)
                                                gccggccgCCCCGGCAAATACTCCACCGACGCGACGGACTAAAATTTTTTCTTACGTTGGGTTTACTCCATCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:8908211/150‑1 (MQ=255)
                                                gCCGGACGCCCGGGCAATTACTCCACCGCCGGGCCGGACTAAAAATTTTTCGTACGGTTGTTTTTCCCCCCCCCCCTGGGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:2851973/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAATTAGTGCCCCGACGCTCCGGACGAAAAATTATTCGTCCGGTGGGTTTTCTCCCTCGCCCTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:7897109/150‑1 (MQ=255)
                                                gCCGGACGACCCGCAAAATAGTCCACCGACGCGTCGCAATAAAAAGTTTACTTATGGTTGTTTTTCCCCCACTCCCTGTGGTTTTGGGGGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:1874161/150‑1 (MQ=255)
                                                  cGGCCGCCCCGCCAAATACTGCACCGACGCGACGGACTAAAAATTATTCTTTAGGTTGTTTTTCTCCCTCCCCGTGTGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                  <  1:10006227/150‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCGAACCCACTCCCACAAGGGAgggggg                                                                                                 >  1:9414026/1‑150 (MQ=38)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCCACACCCGCCCGGGAGAgggg                                                                                                 >  1:9440309/1‑150 (MQ=37)
                                                    gccgccccGCCAAATACTGCACCGCCGCGCCGGACTAAAACTTTTTCTTACGGTTGGTTTCCTCCCTCTCCGTGTGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                <  1:4943092/145‑1 (MQ=255)
                                                     acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCCCCCGAACCCCCACCCGCAAggggggggggcg                                                                                               >  1:1002375/1‑148 (MQ=25)
                                                      ccccccGTCATATCCTCCCCCGCCGCGACGTAACAAAATGTTTCGGTCGGGTTGGTGCACCCCCTCGCCGTGGGGGCGGGGGGGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                                              <  1:3517386/147‑5 (MQ=255)
                                                          ccGGCAACTACTGCCCCGCCGCGACGGACGAATAAGTTTTCGTAAGGTTGGTTGTCTCCCTCGCCCTGGGGGTGAGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGagc                                                                                          <  2:7157760/150‑3 (MQ=255)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCGAACCCACTCCCACAAGGGAGAGGGGACCgcgcg                                                                                         >  1:7667892/1‑146 (MQ=38)
                                                                      gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                             <  2:7574559/119‑1 (MQ=255)
                                                                      gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                             >  1:7574559/1‑119 (MQ=255)
                                                                       cACCGCCGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  1:2263069/150‑1 (MQ=255)
                                                                       cACCGACGCGACGGACTAAAAATTATTCGTATGGTTGGTTTACTCCCTCTCCGTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  1:5592513/150‑1 (MQ=255)
                                                                         cggccgcgactGAATAAAAATTTTTCGTATGGTTTTTTTCCTCCCACTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:5501515/146‑1 (MQ=255)
                                                                         ccGACGCGACGGACTAAAATGTATTCGTTTTTTTTTTTTTCTCCCTCTCCCTGTGGTTTTGGGCGGGGGTGTGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:133771/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:6614011/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:4276064/150‑1 (MQ=39)
                                                                          cgccgccacgTACTAAAATTTTTTCGTATGTTTGGTTTTCCCCCCCTCCGGGTGGGCGCGGGCGGGGGTGAGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                          <  1:8175320/147‑1 (MQ=255)
                                                                           gacgggacggCCGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                         <  2:3298810/150‑1 (MQ=38)
                                                                            acgcgacgGACCAAAAAGTATTCTTATGTTGGGTTTTCTCCCTCGCCCTGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                        <  1:3959626/150‑1 (MQ=255)
                                                                                    gacgAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                <  2:5768985/150‑1 (MQ=255)
                                                                                    gacgAAAAAGTATTCGTACGGTGGGTTTCCTCCCTCGCCCGGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                <  2:3338753/150‑1 (MQ=255)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCACACCCGTAACCACTCCCGCCAGGGTGTGGGGCCCCAAAGAGCACAAACTTTGGGGGGTAtctacttc                                                             >  2:8434627/1‑150 (MQ=18)
                                                                                           atgTATTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  1:9806866/148‑1 (MQ=255)
                                                                                             ttagtCGTATGGTTGTTTTTCCCCCCCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  1:8000603/146‑1 (MQ=255)
                                                                                             gttttCGTAAGGTTTTTTTACTCCCTCTCCCTGTGGGAGTGGGCCTGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  2:4938815/147‑1 (MQ=39)
                                                                                             gTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  1:6204665/150‑1 (MQ=255)
                                                                                             gTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGGGGGTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  1:1391866/150‑1 (MQ=39)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                       <  2:6651900/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCACAACCAAAACCACTACCACCAGGGGGAGGGGGCCCAAAGAGCCCAAAAATTGTGATATAtattattatttagt                                                       >  2:4074820/1‑143 (MQ=12)
                                                                                             ctcttCGTAAGGTGGGCCTCCGCCCGCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCtt                                                       <  2:6631246/147‑1 (MQ=255)
                                                                                              tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                      <  2:5018211/150‑1 (MQ=255)
                                                                                               aTTCGTATGGTTGTTTTTCTCCCACTCCCGGTGGGTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:5048770/150‑1 (MQ=255)
                                                                                               aTTCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGTGTGGGCCGGGGTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  2:5570088/150‑1 (MQ=255)
                                                                                                 tggTACGGTTGGTTTTCCCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  1:8550684/148‑1 (MQ=255)
                                                                                                   gTACGGTTGGTTTTCTCCCTCGCCGTGTGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:2237114/150‑1 (MQ=255)
                                                                                                   gTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGACTTGTAACa                                                 <  2:6802128/150‑1 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTGACCCGCACCCGAACCCGCTACCGCAACGGGGGGGGGACCGCGACAGCACCAAATTTTAGTGGTATATTATTGTGtgttgttaaag                                                >  1:9967196/1‑145 (MQ=12)
                                                                                                     atggttgtttTTGGCCCTCGCCTGGTGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                               <  2:4492279/148‑1 (MQ=255)
                                                                                                     acggtgggttTTCTCCCTCTCCCTGTGGGTGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                               <  1:2832611/148‑1 (MQ=39)
                                                                                                         ttgtttTTCACCCCCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgcagaa                                           <  1:1059934/150‑1 (MQ=39)
                                                                                                         tgggtgTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTGGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  2:9718199/148‑1 (MQ=38)
                                                                                                          tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGGTCACCCACACCCTAACCCCCGCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  2:5101291/150‑1 (MQ=38)
                                                                                                          tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:8972779/150‑1 (MQ=255)
                                                                                                          tggttTTCTCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:7270362/150‑1 (MQ=255)
                                                                                                                tctccctcccccTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                    <  2:4000339/150‑1 (MQ=255)
                                                                                                                    cctctcccaCCGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:1024440/150‑1 (MQ=255)
                                                                                                                     ctctccctGTGGTCGGGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  1:5256324/150‑1 (MQ=255)
                                                                                                                     ctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  1:6791549/150‑1 (MQ=255)
                                                                                                                          cctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCCCACCCCAAACCACCCCCACCAGGGGGTGGGGGCCCATCGAGCCCAAAATTTTTTAGATATAATATTATGTGGTTTTACATAATAACACGAAATAGGGCgttt                          >  2:799170/1‑150 (MQ=255)
                                                                                                                          cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCTATCGAGCACAAATTTTGAGATATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                          <  1:4057155/150‑1 (MQ=255)
                                                                                                                            ggtgGGGGGGGGCGGGGGGGCGGGCCTCCGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCGTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                        <  1:7902527/149‑1 (MQ=255)
                                                                                                                               gggcgcgggCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                     <  1:1885270/144‑1 (MQ=255)
                                                                                                                                ggtgtgggCCGGGGGGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:3874615/145‑1 (MQ=255)
                                                                                                                                ggAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:8366254/150‑1 (MQ=255)
                                                                                                                                   gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCAAACCCTCAACGGGGAGGGGACCGCTAGCGCCCAAAATTTTAGAGATATCTTATTATGTCTTTTAACAGAATAAAAGAAAAAAGGGCTTTTTAGCCAGaa                 >  2:5953532/1‑150 (MQ=255)
                                                                                                                                     gggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata               <  2:2545036/150‑1 (MQ=255)
                                                                                                                                         cGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTGTGCGCCTGAATATTGc           <  2:1961451/150‑1 (MQ=255)
                                                                                                                                          gggtTGTGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCATGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACTGAAGAACAGAATATCGGGCTTTTGCGCCTGAATATTGCg          <  2:4780756/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:2663426/150‑1 (MQ=255)
                                                                                                                                            ggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:3115667/150‑1 (MQ=255)
                                                                                                                                            ggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:1019464/150‑1 (MQ=255)
                                                                                                                                               gcgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa     <  1:7424044/148‑1 (MQ=255)
                                                                                                                                               gAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa     <  1:6914021/150‑1 (MQ=255)
                                                                                                                                                aGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAATTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa    <  1:2285833/150‑1 (MQ=255)
                                                                                                                                                 gggCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCGTGTAACAGAAGAACAGACCATCGGGCATTTGCGCCTGAATATTGCGGAGaaaa   <  2:1142598/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCACACCCTAACCCACTCCCGCAAGGGCGAGGGGCCCGAACGCGCCCAACATTGTAGAGAAATCTTCATATGTCTAGTCACACGCGAACAGCAAAACGGGCTGTTTCGCCTGCAGAATTAGGAGaaaag  >  2:7487401/1‑149 (MQ=255)
                                                                                                                                                  ggCAGCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:5837780/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164432‑164727

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