Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,620 C→T 5.7% intergenic (+86/‑110) 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,6200CT 5.7% 407.1 / 10.2 126intergenic (+86/‑110)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (47/69);  new base T (3/4);  total (53/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 129 total aligned reads displayed.

ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164473‑164769
                                                                                                                                                     |                                                                                                                                                     
aTGTCTGGCCGGACGCCCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTATTCCAACGGTTGGTTTACACCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                     <  2:424138/150‑1 (MQ=38)
 tGTCTGGCCGGACCACCCGCAAAATACGCACCCGACGCGACGTACAAAAAATTTTTCTTATGTTTTTTTTTCCCCCTCCCCGGGGGGGAGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                    <  1:1287468/150‑1 (MQ=255)
  gtgtGGCCGGCCGACCCGCCAAATACTGCACCGACGCGCGGGACTAAAATTTAGTCGTTTGGTTGTTTTTCTCCCTCTCCCTGGGTGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                   <  1:863919/147‑1 (MQ=255)
  gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCCCGCCCCAAACCACACCCACaacg                                                                                                                                                   >  1:3299418/1‑148 (MQ=37)
  gTCGGGCCGGACGACCCGGCAAATACGGCACCGACGCGAGGGACTAAAATTTTTTGGTATGTTTGTTTTTCCCCCTCGCCGTGGGGGTTTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                   <  1:864020/150‑1 (MQ=255)
   tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACggg                                                                                                                                                  >  2:2346451/1‑150 (MQ=255)
   tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCAATCCCACAAggg                                                                                                                                                  >  1:1672150/1‑150 (MQ=37)
       gCCGTACTACCCTGAAAATATTTCACCGAGGCTACGGACTAAAAATTTTTCTTAAGTTTGTTTTTATCCCTCTCCCTTTGGTATTGGCCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                              <  2:2992687/150‑1 (MQ=255)
       gCCGCACTACCCGGAAAATAGTGCACCTACGCGACTGACTAAAAATTTTTCTTACGGTTGTTTTTCTCCCTCTCCATGTGGGTGAGGGCCGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCATCTCCCTTTAGGGagag                                                                                                                                              <  2:1924299/150‑1 (MQ=255)
        ccGGCCGCCCCGGCAAATCCGGCCCCGCCGCGCCGGCCGAACACTTATTCGTAGGGTGGGTTGGCGCCCTCGCCGGGGGGGTGTGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                             <  2:2490384/150‑1 (MQ=255)
          ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCCCGTGCCCCCCCACCCAAACCCCCGCCCCCCAGGGGGTgggg                                                                                                                                           >  2:1312561/1‑150 (MQ=37)
              gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCCCCCCCCCCCCAACCCCCAACCCCAAGGGGGAGGGGGCCg                                                                                                                                       >  2:2460728/1‑150 (MQ=25)
                  cGGCAAATACTGCACCGACGCGACGGACGAAAATTTATTGGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTTAGGGGGGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGTGAGGGGACCGATCg                                                                                                                                   <  1:2555745/150‑1 (MQ=37)
                   ggaaaaTCCTGCACCCCCGCGCGGGACCAAAAAGTTTTCCTATGGTTTGTTTTCGCCCCCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                                                                  <  1:2568859/147‑1 (MQ=255)
                   ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCCAACCCACGCCCACCAGGGCGGGGGGGCCGAtagg                                                                                                                                  >  1:235409/1‑147 (MQ=25)
                       aaTACTGCACCGACGCGAGGGACTATAAAGTTTTGGTTCGGTGGGTTTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                                                              <  1:3979752/150‑1 (MQ=255)
                          aCTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                                                                           <  1:3059213/150‑1 (MQ=255)
                          aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACCCCCCAAACCACACCCACAAAGGGGAGGGGGCCCCTCGCGCACcaa                                                                                                                           >  2:3544571/1‑150 (MQ=21)
                             gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAACAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                                                        <  2:1008722/150‑1 (MQ=255)
                               aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCGAAACCACACCCGCCCGGGCGGGGGGGCCGCTCGCGCCCaaaagttt                                                                                                                      >  1:2738035/1‑145 (MQ=21)
                                 ccgccgCGACGGACTAAAATTTTTGCTTATGGTTGTTTTCCTCCCTCTCCGTGGGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  1:149421/150‑1 (MQ=255)
                                 ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCAACCGAAACCACCACCAACAGGGGGTGGGGGCCGATCCAGCCCCAAtatttt                                                                                                                     >  2:1299635/1‑148 (MQ=18)
                                         ccggacgAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                                                             <  1:1342319/149‑1 (MQ=255)
                                         acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCAAACCCGCAAGGGAGAGGGGGCCGATCCAGCCCAAATTTTGGGAGatata                                                                                                             >  2:607407/1‑149 (MQ=37)
                                                     cttttCTTACGGTGGTTTTCCCCCCTCGCCGGGTGGGTGTGGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                 <  2:3125265/147‑1 (MQ=255)
                                                       aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCCACACCCGCACCGGCGAGGGGACCCATCGAGCACAAAAATTGCGAGATATATTCTTCATTCTTGt                                                                                               >  2:2077266/1‑150 (MQ=21)
                                                         tcttacggttggttTTCTCCCCCGCCCGGGGGGTGAGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGTGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                             <  1:2313049/144‑1 (MQ=255)
                                                             aaGGTTGGTTTTCTCCCTCACCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCATGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                         <  2:2983639/150‑1 (MQ=255)
                                                                 ttgtttTTCGCCCGCGCCGTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                                     <  1:1964254/150‑1 (MQ=255)
                                                                 ttggttTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                                     <  1:1002772/150‑1 (MQ=255)
                                                                   gtttttCTCCCTCTCCGGGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                   <  1:215174/148‑1 (MQ=255)
                                                                   ggttTTCTCCCTCTCCGTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                   <  2:1430688/150‑1 (MQ=255)
                                                                   ggttTTCTCCCACGCCCTCTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                   <  1:4112936/150‑1 (MQ=255)
                                                                          gccctctccgtGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                                            <  1:3886965/149‑1 (MQ=255)
                                                                            cctctccctGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                          <  1:2856379/150‑1 (MQ=255)
                                                                                 ccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                                     <  1:3769649/150‑1 (MQ=255)
                                                                                  ccggtgGGAGGGGGCCGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACCGAAGAACAGAAAATCGGGCtttt                                                                    <  1:1575843/147‑1 (MQ=255)
                                                                                   gtgtgggtgtgggCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                                   <  1:2846230/149‑1 (MQ=255)
                                                                                   ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCGCACCCGAACCCACACCCGCAAAGGAGAGGGGACCGCTCGAGCACAAATTTTGTGAGATATAATCTTATTTCTTTTTACAGAAGAACAGACAAAAGGGCgtttt                                                                   >  2:3647755/1‑149 (MQ=255)
                                                                                   ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCAAACCCACAACCGCAAGGGGGAGGGGGCCGATCGGGCCCAAATTTTGAGAGATATATTATTCTGTCGTTTTACAGAAGAACAGGAAATGGGGCgtttt                                                                   >  2:4051447/1‑149 (MQ=19)
                                                                                   ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTCCCACAAGGGAGTGGGGACCCATAGAGCACAAAAGTTGAGTGATATAATCTTAAGTATTGTAACACAAAAACAGAAAAAAGGGGttttt                                                                   >  2:4051601/1‑149 (MQ=255)
                                                                                    tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCGCCAGGGCGTGGGGGCCGCACGTGCCCAAAAGTTTAGTGATATATTATTATGTCGTGTAAAAGAAAAAAAAAAAAAAGGGCGTTTgc                                                                  >  2:2044135/1‑150 (MQ=255)
                                                                                     gtgGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                                                                 <  2:682178/150‑1 (MQ=255)
                                                                                     gtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                                                                 <  1:1726931/150‑1 (MQ=255)
                                                                                       gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAAACACACCCAAAAAGGAGGGGGGCCCGATGGAGCACAAAATTTGAGAGAAATATTAATTTGTGTTTGAAAAGTAGAACAGGAAATAGGGCtttttgggg                                                               >  2:3726507/1‑145 (MQ=255)
                                                                                        ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                              <  1:3108117/150‑1 (MQ=255)
                                                                                         gagcGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                                                             <  1:2067073/150‑1 (MQ=255)
                                                                                         gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCAAAACCACTACCGCAACGGGGGGGGGACCCATCGAGCACAAAATTTGAGGGATATATTATTAAGTCTTTTAAAAGAAGAACAGAAAAAAGGGCtttttagccag                                                             >  2:1019019/1‑143 (MQ=255)
                                                                                            agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCCAACCCACACCCAAAAGGGGGAGGGGACCGATCGCGCACAAAATTTTTGAGATAAATTAATAATTGTTGTTACAGAAGAACACAAAATACGGCTTTTTGGCCGGAat                                                          >  2:1697249/1‑150 (MQ=255)
                                                                                              ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTCCCACAAGGGAGAGGGGGCCGATCGAGCACAAATTTTGAGAGATATCATCATATGTCGTGTAACAGTAGAACACACAATCGGGCGTTTGCGCCGGAaaat                                                        >  2:1391039/1‑150 (MQ=255)
                                                                                              ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCAAAACCAAACCCACACCCGCAAGGGCGAGGGGGCCGATCGAGCACAAATTTTTAGAGATATAATCTTCAGTCTTGTAAAATAAGAACAGAAAATACGGCTTTTTCGCCGGAatat                                                        >  2:3445941/1‑150 (MQ=255)
                                                                                               gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCCCTCCCACGAGGGGGAGGGGCCCCAAAGCGCCCCAATTGTGTGATAAAAAATCTTCTGGGTTGTAAGACAACAACCGACAAACGGGCGTTTTGGGCGGAATAtt                                                       >  2:3026520/1‑150 (MQ=255)
                                                                                                 cGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTCGCGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGACAATCGGGCTTTTGCGCCGGAATATTGc                                                     <  2:3240225/150‑1 (MQ=255)
                                                                                                  ggggggCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:42854/150‑1 (MQ=255)
                                                                                                  ggggTGGAGGCCCTCAGCGAGCACGTTCACCCTCAACCTAACCCTCTCCCTCAAGGGAGAGGGGACCGACCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTATCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:3796051/150‑1 (MQ=255)
                                                                                                  ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:2897497/150‑1 (MQ=255)
                                                                                                  ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:2346451/150‑1 (MQ=255)
                                                                                                   gggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg                                                   <  2:1357092/150‑1 (MQ=255)
                                                                                                     ggtttgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                  <  2:789661/145‑1 (MQ=255)
                                                                                                     gTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCTCACCCACAAGGGAGAGGGGGCCCATCCCGCCCAAAATTTTTGGGTTATATTATTAAGTGGTTGAACAGAAGAACAGCAAAAAGGGCGTTTTCGCCCGGAAATTTCggcg                                                 >  2:3325314/1‑148 (MQ=255)
                                                                                                           ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                            <  2:2494200/150‑1 (MQ=255)
                                                                                                           gccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                            <  1:3201739/148‑1 (MQ=255)
                                                                                                                           ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  2:302558/150‑1 (MQ=255)
                                                                                                                            tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAc                           <  1:1173682/150‑1 (MQ=255)
                                                                                                                              accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                         <  1:1302820/150‑1 (MQ=255)
                                                                                                                               ccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGTCCGATCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGTAAAAGCATGGCCGGGAATGACcgc                        <  2:2727752/150‑1 (MQ=255)
                                                                                                                                cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgcg                       >  1:2828505/1‑150 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:1641908/150‑1 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:1753129/150‑1 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCGGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  2:2708203/1‑150 (MQ=255)
                                                                                                                                  tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAg                     >  1:2669597/1‑150 (MQ=255)
                                                                                                                                  tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     <  2:1780613/150‑1 (MQ=255)
                                                                                                                                   caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  2:2188273/150‑1 (MQ=255)
                                                                                                                                     ccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCa                  <  2:1403191/150‑1 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCac                 >  2:2685535/1‑149 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:1206220/1‑150 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAAATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:3283628/1‑150 (MQ=255)
                                                                                                                                       ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCact                >  2:548951/1‑148 (MQ=255)
                                                                                                                                       ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt                >  2:548796/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               >  1:2512305/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGAATGGCCGGGAATGACCGCGAGCCAAtt               >  1:2512360/1‑150 (MQ=255)
                                                                                                                                         tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              >  1:2387500/3‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  1:3473611/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  1:663092/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTCCCGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCACTTggcc           >  2:831648/1‑148 (MQ=255)
                                                                                                                                             ctctCCCTCAAGGGAGGGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  2:3540799/1‑150 (MQ=255)
                                                                                                                                             ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  2:2341100/1‑150 (MQ=255)
                                                                                                                                               ctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGACAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCa        >  1:1681592/1‑150 (MQ=255)
                                                                                                                                               ctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCa        >  2:109427/1‑150 (MQ=255)
                                                                                                                                                 cccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaaa      <  1:1019019/150‑1 (MQ=255)
                                                                                                                                                  ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  2:1760347/150‑1 (MQ=255)
                                                                                                                                                  ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  1:1222080/150‑1 (MQ=255)
                                                                                                                                                  ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAACAACACAAAATCGGGCTTTTGCGCCAGAATATTGCGGACAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  1:2394456/150‑1 (MQ=255)
                                                                                                                                                   cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    <  2:2412512/150‑1 (MQ=255)
                                                                                                                                                   cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  1:762517/1‑150 (MQ=255)
                                                                                                                                                   cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    <  1:3503293/150‑1 (MQ=255)
                                                                                                                                                   cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    <  1:2125531/150‑1 (MQ=255)
                                                                                                                                                   cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATAGCGGAGAAAAATCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  2:1011954/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:3045671/1‑150 (MQ=255)
                                                                                                                                                     cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa  >  1:71213/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164473‑164769

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: