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.

ATGTCTGGCCGGACGACCCGGCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164473‑164769
                                                                                                                                                     |                                                                                                                                                     
aTGTCTGGCCGGACGCCCCGGCAAATCCTGCA‑CCGACGCGACGGACGAAAAAGTATTCCAACGGTTGGTTTACACCCTCTCCCTGTGGGAGTGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                     <  2:424138/150‑1 (MQ=38)
 tGTCTGGCCGGACCACCCGCAAAATAC‑GCACCCGACGCGACGTACAAAAAATTTTTCTTATGTTTTTTTTTCCCCCTCCCCGGGGGGGAGGGGGCGGGGGT‑GCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                    <  1:1287468/150‑1 (MQ=255)
  gtgtGGCCGGCCGACCCGCCAAATACTGCA‑CCGACGCGCGGGACTAAAATTTAGTCGTTTGGTTGTTTTTCTCCCTCTCCCTGGGTGTGTGGGCGGGGGT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                   <  1:863919/147‑1 (MQ=255)
  gTCTGGCCGGACGACCCGGCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGCGCCCCCGCCCCAAACCACACCCACaacg                                                                                                                                                   >  1:3299418/1‑148 (MQ=37)
  gTCGGGCCGGACGACCCGGCAAATACGGCA‑CCGACGCGAGGGACTAAAATTTTTTGGTATGTTTGTTTTTCCCCCTCGCCGTGGGGGTTTGGGCGGGGGT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                   <  1:864020/150‑1 (MQ=255)
   tCTGGCCGGACGACCCGGCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACggg                                                                                                                                                  >  2:2346451/1‑150 (MQ=255)
   tCTGGCCGGACGACCCGGCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCAATCCCACAAggg                                                                                                                                                  >  1:1672150/1‑150 (MQ=37)
       gCCGTACTACCCTGAAAATATTTCA‑CCGAGGCTACGGACTAAAAATTTTTCTTAAGTTTGTTTTTATCCCTCTCCCTTTGGTATTGGCCCGGGGT‑GTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                              <  2:2992687/150‑1 (MQ=255)
       gCCGCACTACCCGGAAAATAGTGCA‑CCTACGCGACTGACTAAAAATTTTTCTTACGGTTGTTTTTCTCCCTCTCCATGTGGGTGAGGGCCGGGTT‑GTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCATCTCCCTTTAGGGagag                                                                                                                                              <  2:1924299/150‑1 (MQ=255)
        ccGGCCGCCCCGGCAAATCCGGCC‑CCGCCGCGCCGGCCGAACACTTATTCGTAGGGTGGGTTGGCGCCCTCGCCGGGGGGGTGTGGGCCGGGGG‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                             <  2:2490384/150‑1 (MQ=255)
          ggACGACCCGGCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCCGCGCGCCCGTGCCCCCCCACCCAAACCCCCGCCCCCCAGGGGGTgggg                                                                                                                                           >  2:1312561/1‑150 (MQ=37)
              gacCCGGCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCCTCAGCGCGCACGGGCCCCCCCCCCCCAACCCCCAACCCCAAGGGGGAGGGGGCCg                                                                                                                                       >  2:2460728/1‑150 (MQ=25)
                  cGGCAAATACTGCA‑CCGACGCGACGGACGAAAATTTATTGGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTTAGGGGGGGGGT‑TAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGTGAGGGGACCGATCg                                                                                                                                   <  1:2555745/150‑1 (MQ=37)
                   ggaaaaTCCTGCA‑CCCCCGCGCGGGACCAAAAAGTTTTCCTATGGTTTGTTTTCGCCCCCTCCCTGTGGGTGTGGGCGGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                                                                  <  1:2568859/147‑1 (MQ=255)
                   ggCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGGACCCGCACCCCAACCCACGCCCACCAGGGCGGGGGGGCCGAtagg                                                                                                                                  >  1:235409/1‑147 (MQ=25)
                       aaTACTGCA‑CCGACGCGAGGGACTATAAAGTTTTGGTTCGGTGGGTTTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                                                              <  1:3979752/150‑1 (MQ=255)
                          aCTGCA‑CCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                                                                           <  1:3059213/150‑1 (MQ=255)
                          aCTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGGGCCCACCCCCCAAACCACACCCACAAAGGGGAGGGGGCCCCTCGCGCACcaa                                                                                                                           >  2:3544571/1‑150 (MQ=21)
                             gCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCAACAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                                                        <  2:1008722/150‑1 (MQ=255)
                               a‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGCCCCCACACCCGAAACCACACCCGCCCGGGCGGGGGGGCCGCTCGCGCCCaaaagttt                                                                                                                      >  1:2738035/1‑145 (MQ=21)
                                 ccgccgCGACGGACTAAAATTTTTGCTTATGGTTGTTTTCCTCCCTCTCCGTGGGGTAGTGGGCGGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  1:149421/150‑1 (MQ=255)
                                 ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTGACCCGCAACCGAAACCACCACCAACAGGGGGTGGGGGCCGATCCAGCCCCAAtatttt                                                                                                                     >  2:1299635/1‑148 (MQ=18)
                                         ccggacgAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGAGCGGGCGGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                                                             <  1:1342319/149‑1 (MQ=255)
                                         acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCAAACCCGCAAGGGAGAGGGGGCCGATCCAGCCCAAATTTTGGGAGatata                                                                                                             >  2:607407/1‑149 (MQ=37)
                                                     cttttCTTACGGTGGTTTTCCCCCCTCGCCGGGTGGGTGTGGGCGGGGTT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                 <  2:3125265/147‑1 (MQ=255)
                                                       aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCCACACCCGCACCGGCGAGGGGACCCATCGAGCACAAAAATTGCGAGATATATTCTTCATTCTTGt                                                                                               >  2:2077266/1‑150 (MQ=21)
                                                         tcttacggttggttTTCTCCCCCGCCCGGGGGGTGAGGGCGGGGGG‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGTGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                                             <  1:2313049/144‑1 (MQ=255)
                                                             aaGGTTGGTTTTCTCCCTCACCCTGTGGGAGTGGGCCGGGGT‑GTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCATGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                         <  2:2983639/150‑1 (MQ=255)
                                                                 ttgtttTTCGCCCGCGCCGTGGGGGTGTGGGCGGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                                     <  1:1964254/150‑1 (MQ=255)
                                                                 ttggttTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                                                     <  1:1002772/150‑1 (MQ=255)
                                                                   gtttttCTCCCTCTCCGGGTGGGTGTGGGCGGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                   <  1:215174/148‑1 (MQ=255)
                                                                   ggttTTCTCCCTCTCCGTGTGGGAGAGGGCGGGGGT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                   <  2:1430688/150‑1 (MQ=255)
                                                                   ggttTTCTCCCACGCCCTCTGGGAGAGGGCCGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                   <  1:4112936/150‑1 (MQ=255)
                                                                          gccctctccgtGGGGTTGTGGGCGGGGGT‑GTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                                            <  1:3886965/149‑1 (MQ=255)
                                                                            cctctccctGTGGGTGTGGGCGGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                          <  1:2856379/150‑1 (MQ=255)
                                                                                 ccctGTGGGAGAGGGCCGGGGT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                                     <  1:3769649/150‑1 (MQ=255)
                                                                                  ccggtgGGAGGGGGCCGGGGT‑TCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACCGAAGAACAGAAAATCGGGCtttt                                                                    <  1:1575843/147‑1 (MQ=255)
                                                                                   gtgtgggtgtgggCGGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                                   <  1:2846230/149‑1 (MQ=255)
                                                                                   ctGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCGCCCGCACCCGAACCCACACCCGCAAAGGAGAGGGGACCGCTCGAGCACAAATTTTGTGAGATATAATCTTATTTCTTTTTACAGAAGAACAGACAAAAGGGCgtttt                                                                   >  2:3647755/1‑149 (MQ=255)
                                                                                   ctGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCGCCCCCAAACCCACAACCGCAAGGGGGAGGGGGCCGATCGGGCCCAAATTTTGAGAGATATATTATTCTGTCGTTTTACAGAAGAACAGGAAATGGGGCgtttt                                                                   >  2:4051447/1‑149 (MQ=19)
                                                                                   ctGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTCCCACAAGGGAGTGGGGACCCATAGAGCACAAAAGTTGAGTGATATAATCTTAAGTATTGTAACACAAAAACAGAAAAAAGGGGttttt                                                                   >  2:4051601/1‑149 (MQ=255)
                                                                                    tgtgGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCGCCAGGGCGTGGGGGCCGCACGTGCCCAAAAGTTTAGTGATATATTATTATGTCGTGTAAAAGAAAAAAAAAAAAAAGGGCGTTTgc                                                                  >  2:2044135/1‑150 (MQ=255)
                                                                                     gtgGGAGAGGGCGGGGGT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                                                                 <  2:682178/150‑1 (MQ=255)
                                                                                     gtgGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                                                                 <  1:1726931/150‑1 (MQ=255)
                                                                                       gggAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAAACACACCCAAAAAGGAGGGGGGCCCGATGGAGCACAAAATTTGAGAGAAATATTAATTTGTGTTTGAAAAGTAGAACAGGAAATAGGGCtttttgggg                                                               >  2:3726507/1‑145 (MQ=255)
                                                                                        ggAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                              <  1:3108117/150‑1 (MQ=255)
                                                                                         gagcGGGCCGGGGT‑GTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                                                             <  1:2067073/150‑1 (MQ=255)
                                                                                         gagaGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCACAACCAAAACCACTACCGCAACGGGGGGGGGACCCATCGAGCACAAAATTTGAGGGATATATTATTAAGTCTTTTAAAAGAAGAACAGAAAAAAGGGCtttttagccag                                                             >  2:1019019/1‑143 (MQ=255)
                                                                                            agGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCACAACCCAACCCACACCCAAAAGGGGGAGGGGACCGATCGCGCACAAAATTTTTGAGATAAATTAATAATTGTTGTTACAGAAGAACACAAAATACGGCTTTTTGGCCGGAat                                                          >  2:1697249/1‑150 (MQ=255)
                                                                                              ggCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTCCCACAAGGGAGAGGGGGCCGATCGAGCACAAATTTTGAGAGATATCATCATATGTCGTGTAACAGTAGAACACACAATCGGGCGTTTGCGCCGGAaaat                                                        >  2:1391039/1‑150 (MQ=255)
                                                                                              ggCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCAAAACCAAACCCACACCCGCAAGGGCGAGGGGGCCGATCGAGCACAAATTTTTAGAGATATAATCTTCAGTCTTGTAAAATAAGAACAGAAAATACGGCTTTTTCGCCGGAatat                                                        >  2:3445941/1‑150 (MQ=255)
                                                                                               gCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCCCTCCCACGAGGGGGAGGGGCCCCAAAGCGCCCCAATTGTGTGATAAAAAATCTTCTGGGTTGTAAGACAACAACCGACAAACGGGCGTTTTGGGCGGAATAtt                                                       >  2:3026520/1‑150 (MQ=255)
                                                                                                 cGGGGG‑GAGGGCCTCAGCGCGCACGGTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTCGCGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGACAATCGGGCTTTTGCGCCGGAATATTGc                                                     <  2:3240225/150‑1 (MQ=255)
                                                                                                  ggggg‑gCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:42854/150‑1 (MQ=255)
                                                                                                  ggggTGGAGGCCCTCAGCGAGCACGTTCACCCTCAACCTAACCCTCTCCCTCAAGGGAGAGGGGACCGACCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTATCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:3796051/150‑1 (MQ=255)
                                                                                                  ggggT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:2897497/150‑1 (MQ=255)
                                                                                                  ggggT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:2346451/150‑1 (MQ=255)
                                                                                                   gggT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg                                                   <  2:1357092/150‑1 (MQ=255)
                                                                                                     ggtttgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                                  <  2:789661/145‑1 (MQ=255)
                                                                                                     gT‑GAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCTCACCCACAAGGGAGAGGGGGCCCATCCCGCCCAAAATTTTTGGGTTATATTATTAAGTGGTTGAACAGAAGAACAGCAAAAAGGGCGTTTTCGCCCGGAAATTTCggcg                                                 >  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)
                                                                                                                                                     |                                                                                                                                                     
ATGTCTGGCCGGACGACCCGGCAAATACTGCA‑CCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  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: ‑