Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,620 C→T 6.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 6.7% 571.5 / 20.1 182intergenic (+86/‑110)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (48/118);  new base T (9/3);  total (60/122)
Fisher's exact test for biased strand distribution p-value = 2.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

ACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164471‑164769
                                                                                                                                                     |                                                                                                                                                     
  aTGTCTGGCCGGACGCCCCGGCAAATACTGCCCCGACGCGACGGACGACAAAGTAGTCGTCCGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                     <  2:4887186/150‑1 (MQ=37)
   tGTCGGGCCGGCCGACCCGGAAAATCCTGCCCCGACGCGACGGACTAAAAAGTATTCTTACGGTTGGTTTTCACCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                    <  2:332692/150‑1 (MQ=255)
   ggtgtGGCCGGACGACCGGGAAAATACTGCACCGCGGCGACGGACTAAAATGTATTCTTATGTTTGGTTTTCTCCCTCTCCGTGTGGTAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                                    <  2:1711019/146‑1 (MQ=255)
       tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACAACCACAAGGgag                                                                                                                                                >  2:4504978/1‑150 (MQ=38)
       tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACACCCACAAGGgag                                                                                                                                                >  2:4505439/1‑150 (MQ=38)
        ggCCGGACGACCCGGCAATTACTGCACGGACGCGACGGACGAAAAAGTTTTGGTAAGGTTGGTTTCCTCCCTCGCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                                               <  2:4814256/150‑1 (MQ=255)
        ggCCGGACGACCCGGCAAACACTGCCCCGACGCGACGGACGAAAAATTATTCGTAAGGTGGGTTTTCTCCCTCCCCGTGTGGGTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                                               <  2:4156640/150‑1 (MQ=255)
         gcggtacgccCCTGAAAATACGGCACGGCCGCGCCGGACGAAAATTTTTTCTTTTGTTTGTTTTTCACCCTCTCCGTGTGGTTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                              <  1:3684220/145‑1 (MQ=255)
         gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCGCACCCCAAACCACGCCCGCAAggggggg                                                                                                                                              >  1:1478369/1‑146 (MQ=37)
         gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                              <  2:4864286/150‑1 (MQ=255)
         gCCGGACCCCCCGCAAAATCCTGCACCCACGCGACGGACTAAAAATTTTTCGTATGTTTGGTTTTCCCCCTCTCCCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                              <  2:5647137/150‑1 (MQ=255)
           cgtacgacCCCTAAAATATTCCACCTACGCGCCTCACTAAAATTTTTTCTTATGTTTGTTTTTGCCCCTCTCCTTGTGGTTGTGGGCGGGGGTTTGGGCATCAGCGCGCACGTTCACACTCACCCTAACACTCTCACTCAAGGGAGAggg                                                                                                                                            <  2:5054183/147‑1 (MQ=255)
            ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACACCCAAAACCACTCCCACAAGGGAGAgggg                                                                                                                                           >  2:1498538/1‑150 (MQ=38)
             gacgccCCGGAAAATACGGCCCCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCGCGCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                          <  2:729266/150‑1 (MQ=38)
             gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCCCCCCAACCCGCGCCCGCCACGGGGAggggg                                                                                                                                          >  2:5470533/1‑149 (MQ=37)
                   ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCTAACCCTATCCCTAAATGGTGGGGGGGGCGAtg                                                                                                                                    >  2:5479791/1‑149 (MQ=37)
                   ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAACCAAAACCACAAGGGGGAGGGGccccaac                                                                                                                                    >  2:3833282/1‑146 (MQ=25)
                    cGGCAAATACTGCACCGACGCGACGGACGACAAAGTATTCGTCAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCGCACCCCAACCCGCGCCCCCCCGGGCGCGGGGAccccgcg                                                                                                                                   >  2:2491752/1‑145 (MQ=21)
                    cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCAAACCCACACCCACAACGGAGCGGGGGCCGAACg                                                                                                                                   >  1:5740558/1‑150 (MQ=25)
                      gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                                                 <  2:4435827/150‑1 (MQ=255)
                       aaaaTCCTCCACCGACGCGACGGACTAAAATTTATTCTTCTGTTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                                                <  1:2872994/149‑1 (MQ=255)
                        aaaTAGTGCACCGCCGCTACGCACGAAAATTTTTGCGTCTGTTTGTTTTTCGCCCTCGCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                                               <  2:200292/150‑1 (MQ=255)
                        aaaTACTTCACCGAGGCGCCTGACTAAAATTTTTTCTTATGTTTGTTTTTCTCCATACCCCTGTGGGAGTGGGCCGGGTTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                                               <  2:200355/150‑1 (MQ=255)
                             cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAAACCACACCCTCCCGGGAGCGGGGGCCGGTCGCGCCCCAAt                                                                                                                          >  1:5880434/1‑150 (MQ=25)
                                cACCGCCGCGACGGACGAAAATGTTTTCTTATGGTTGGTTTTCGCCCTCTCCGGGTGGGAGTGGGCGGGGGGGTGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                       <  1:1362674/150‑1 (MQ=255)
                                 acggccgcgacgGACGAAAATTTATAGGTATGGTGGGTTTCCCCCCGCTCCGGGGGGGATCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                                      <  1:3263290/145‑1 (MQ=255)
                                  ccgccgCGACGCACTAAAAAGTTTTCGTACGGTGGGTTTACGCCCTCCCCCTGGGGGTGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  2:2246253/150‑1 (MQ=255)
                                  ccgccgCCACGGACCAAAAATTTTTCGTCCGGTGGTTTTTCCCCCCCTCCCGGGGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  1:4646319/150‑1 (MQ=255)
                                  ccGACGCGACGGACTAAAATGTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  2:3654856/150‑1 (MQ=38)
                                  ccGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  2:4451034/150‑1 (MQ=255)
                                  ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCGGGGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  2:3654599/150‑1 (MQ=255)
                                   ggccgcgacgGACCAAATTGTAGTCGTACGGTTGGTTTTCCCCCTCTCCGTGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                                                    <  1:4757748/147‑1 (MQ=255)
                                    gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCCCACCCGAAACCACACACACAAGGGAGGGGGGGCCCAGAGAGCACCAAAATTgaga                                                                                                                   >  1:5817471/1‑150 (MQ=21)
                                            ggacgAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                                                           <  2:3395320/150‑1 (MQ=255)
                                             gacgAACCAGTATGCGTACGGTGGGTTTTCTCCCTGGCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                                                          <  2:3525724/150‑1 (MQ=255)
                                                    aattATTCGGTAGGTTGGTTTGCCCCCTCACCCTGGGGTTGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                                                   <  2:4111279/147‑1 (MQ=255)
                                                      gTATTGGTTGGGTTGGTGTTCTCCCTCGCCGGGGGGGTGTGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                 <  2:609880/150‑1 (MQ=255)
                                                      gTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                 <  2:3576614/150‑1 (MQ=255)
                                                            gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCTGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCCAAACCGCGACCCCAAGGGAGAGGGGGCCGCGCGCGCCCAAAATTTGTGAGCTATCGTCTTCTGTCGTGTGACa                                                                                           >  2:3495373/1‑150 (MQ=16)
                                                               cggtgggttTTCGCCCGCGCCGTGTGGGCGCGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                        <  2:1704955/149‑1 (MQ=255)
                                                                 gttgtttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                                      <  1:4505283/150‑1 (MQ=255)
                                                                 gttggttTTCTCCCTCTCCCTGTGGCTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                                      <  1:4504978/150‑1 (MQ=255)
                                                                    ggttTTCGCCATCTCCGGGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                   <  1:2760826/150‑1 (MQ=255)
                                                                         tctccctctcccGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                                              <  2:1446265/150‑1 (MQ=255)
                                                                           cccctcaccccGGGGGAGCGGCCCGGGGCGCGGCCATCAGCGCGCACCTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTGTGAGAGAGATCTTCCTCTGTCTTGTAACAGAAGAACAGAAAATc                                                                            <  2:243296/149‑1 (MQ=255)
                                                                             cctctccgtGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                          <  1:3277824/150‑1 (MQ=255)
                                                                                  ccctGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                                     <  1:2847954/150‑1 (MQ=255)
                                                                                   cctGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                                    <  2:5817471/150‑1 (MQ=255)
                                                                                      gggggCGAGGGCCGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                                                                 <  2:2384660/148‑1 (MQ=255)
                                                                                       ggggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                                <  1:1386496/149‑1 (MQ=255)
                                                                                              gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                                                         <  2:4376385/150‑1 (MQ=255)
                                                                                                 cgggggggTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACATAAAATCGGGCTTTTGCGCCTGAATATTg                                                      <  1:734269/148‑1 (MQ=255)
                                                                                                  cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  2:1478369/150‑1 (MQ=255)
                                                                                                   ggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:5550300/150‑1 (MQ=255)
                                                                                                   ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:1950152/150‑1 (MQ=255)
                                                                                                   ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:5995372/150‑1 (MQ=255)
                                                                                                          gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                                             <  2:4623265/150‑1 (MQ=255)
                                                                                                            gcctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                           <  2:2565258/147‑1 (MQ=255)
                                                                                                                cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg                                       <  2:4779581/150‑1 (MQ=255)
                                                                                                                    gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                                   <  1:3694863/150‑1 (MQ=255)
                                                                                                                     cgcACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                  <  1:4132333/150‑1 (MQ=255)
                                                                                                                          gTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTATTGTAACACAAGAACAGAAAATCGGGCTTTTGCGCCAGAATATTGCGGAGAAAAAGCATGGCCGGGAATg                             >  2:4903144/1‑150 (MQ=255)
                                                                                                                              accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                         <  1:3502273/150‑1 (MQ=255)
                                                                                                                               ccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgc                        <  2:5140642/150‑1 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  2:318018/1‑150 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  2:3057702/1‑150 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:1811995/1‑150 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:884742/1‑150 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGACGGGACCGCTCGAGCACAAATTTTGAGAGCTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGTGAAAAAGCATGGCCGGGAATGACCGCGa                      <  1:5345142/150‑1 (MQ=255)
                                                                                                                                   caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  2:5702485/150‑1 (MQ=255)
                                                                                                                                   caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  2:5702380/150‑1 (MQ=255)
                                                                                                                                   caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  1:3363177/150‑1 (MQ=255)
                                                                                                                                   caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  1:3544596/150‑1 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:3562917/1‑150 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  1:1778293/1‑150 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:261209/1‑150 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  1:2200576/1‑150 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  2:4957678/1‑150 (MQ=255)
                                                                                                                                       ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt                >  2:3739635/1‑150 (MQ=255)
                                                                                                                                       ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt                <  2:3448643/150‑1 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  1:4090540/150‑1 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  1:25518/150‑1 (MQ=255)
                                                                                                                                         tccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              >  1:180958/3‑150 (MQ=255)
                                                                                                                                         tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTg              >  2:854584/3‑150 (MQ=255)
                                                                                                                                         tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              >  2:2923503/3‑150 (MQ=255)
                                                                                                                                           cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            >  2:1733547/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAATGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  2:5568424/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTGGAc           >  2:4062265/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  1:3902196/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  2:2313152/1‑150 (MQ=255)
                                                                                                                                             ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          <  1:1337007/150‑1 (MQ=255)
                                                                                                                                             ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          <  1:1337003/150‑1 (MQ=255)
                                                                                                                                              tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc         >  1:421002/1‑150 (MQ=255)
                                                                                                                                                acccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTAGAGAGATATCTTCTTCTGTCTTGAAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa       <  2:710523/149‑1 (MQ=255)
                                                                                                                                                 cccTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGGAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAAGTGGACGCaaa      >  2:5362818/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATTGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:2360722/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:5017281/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:1798870/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   <  1:491047/150‑1 (MQ=255)
                                                                                                                                                    tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACACAAAggg   >  1:189001/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164471‑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: