Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,620 C→T 6.6% 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.6% 832.3 / 32.5 259intergenic (+86/‑110)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (78/164);  new base T (9/8);  total (87/172)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 265 total aligned reads displayed.

ACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGA  >  NC_000913/164471‑164769
                                                                                                                                                     |                                                                                                                                                     
aCATGTCTGGCCGGACGACCCGGAAAATCCTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCCCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                                                                       <  1:8351824/150‑1 (MQ=255)
  ttttCTGGCCGCACGCCCGGCCAAATACTGCCCCGACGCGCGGGACTAAAATTTTTTCTTATGTTTGGTTTTCTCCCTCCCCGTGTGGTCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                     <  1:8320468/147‑1 (MQ=255)
  ttgtgtgGCCGGCCGACCCGCCAAATCCTGCACCGCCGCGACGGACGAAAATTTTTTCGTTTGGTTGGTTTACTCCCTCGCCCTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                     <  1:2978453/149‑1 (MQ=255)
  ttgtggggCCGCACGACCCGCAAAATAGTGCCCCGCGGGGCGGGACTAAAATTTTTTCTTTTGTTTGTTTTTCCCCCTCTCCGTGTGGGTGTGGGGGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                     <  2:2529311/144‑1 (MQ=255)
  aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAACTGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGGCCCCCACCCGGTCCCTCACCCTccc                                                                                                                                                     >  2:6249335/1‑148 (MQ=25)
    gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                                   <  2:5873264/150‑1 (MQ=255)
        ggccggccGCCCCGGCAAATACGGCACCGACGCGCCGGACGAAAAAGTATTCGTACGTTTGGTTTTCCCCCACTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                                                               <  2:5745192/150‑1 (MQ=37)
        ggCCGCACTACCGGCAAAATACGGCACCGACGCGAGGTACTAAAATTTTTTCTTACGGTTGTTTTACTCCCTCCCCCTTTGGGTTTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                                               <  1:3122354/150‑1 (MQ=255)
         gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCGAAACCACTCCCTCAAggggggg                                                                                                                                              >  1:3856213/1‑146 (MQ=37)
         gCCGGACGACCCGCCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGTTTTTCTCCCCCTCCCTGTGGTAGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                              <  2:7265758/150‑1 (MQ=255)
          ccGGACGACCGGTAAAATCCTGCACCGACGCGCCGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                             <  1:7704719/150‑1 (MQ=37)
           cgtacgacCCGGCAAATACTGCACCGAGGCGACGTACTAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                                            <  1:8475060/147‑1 (MQ=255)
            ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTGGCCCCCCACCCCAACCCACACCCACAACGGCgggggg                                                                                                                                           >  1:8130761/1‑150 (MQ=25)
            ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCACCCAAAACCAAACCCACAAGGGGGCgggg                                                                                                                                           >  1:583176/1‑150 (MQ=25)
              ccgacCCGGCAACAACGGCCCCGACGCGCCGGCCGAACAAGTTTTCGTACGGTTGGTTTTCCCCCTCTCCCTGGGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggtta                                                                                                                                         <  2:8325345/149‑4 (MQ=255)
              ccgacCCGGCAAACACTGCACCGACGCGACGGACGAAAATGTATTCGTACGGTTGTTTTTCTCCCTCCCCGTGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggtta                                                                                                                                         <  2:4468551/149‑4 (MQ=255)
              acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGAAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggtta                                                                                                                                         <  2:8545718/150‑4 (MQ=255)
                cacCCGCAAAATACCCCCCCGCCGCGCGGGACTAAAATGTATTCGTTCGGTGGGTTTTCCCCGTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                                                                       <  1:5214362/149‑6 (MQ=255)
                   ccgcaaaaTACTGCCCCGACGCGAGTGACTAAAATTTTGTGTTTTGTTTGGTTTTCCCCCTCGCCCGGTGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                                                    <  1:4896911/145‑1 (MQ=255)
                             cTGCCCCGACGCGCCGGACGAAAAAGTATTCGTACGGTTGGTTGTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                                                          <  2:7718288/150‑1 (MQ=39)
                             cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGCCCCGCACCCGAACCCACACCCACCAGGGCGGGGGGGCCCAGCGAGcccccaa                                                                                                                          >  1:4533980/1‑144 (MQ=21)
                             cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGTCACCCTCACCCGAAACCACTCCCACAAGGGGGAGGGGACCGATCGAGCACACAt                                                                                                                          >  2:6896208/1‑150 (MQ=37)
                                cACCGACGCGAGGGACGAAAATTTTTGCTTACGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                       <  1:3938032/150‑1 (MQ=37)
                                  ccGACGCGCCGGACGACAACTTTTTGGTATGTTTGTTTTTCTCCCTCTCCGGGTGGTCGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                     <  2:1864697/150‑1 (MQ=255)
                                   cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                                                    <  2:5999492/150‑1 (MQ=255)
                                       cctacgcaccAAAAATTATTCTAAATTTTGTTTTACTCCCTCTCCCTGTGGGAGAGGCCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                                                <  2:1309144/147‑1 (MQ=37)
                                             gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCAACCCAAACCCAAAACCACAAGGGAGAGGGGACCCATACCGCACAAAATTTTAGATATAtatta                                                                                                          >  1:6929490/1‑149 (MQ=18)
                                               cgAAAAAGTTTGCGTATGGTTGGTTTTCTCCCTCTCCGGGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                                                        <  1:5983761/150‑1 (MQ=255)
                                               cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                                                        <  2:9187416/150‑1 (MQ=255)
                                                 aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACACCCACAACGGGGAGGGGACCCATCCAGCACAAAATTTTATAGATAtcttaataa                                                                                                      >  2:4679162/1‑145 (MQ=21)
                                                   aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                                    <  1:9183101/150‑1 (MQ=255)
                                                    acgTATGCGTAGGGTGGGTTTTCGCCCGCCCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                                                   <  2:7507718/148‑1 (MQ=255)
                                                      gTATTCGTATGGTTGGTTTTCTCCCTCGCCGTGGGGGTGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                 <  1:533941/150‑1 (MQ=255)
                                                      gTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                 <  1:6661782/150‑1 (MQ=255)
                                                               tgtttgtttTACCCCCTGGCCGTGGGGGAGGGGGGGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                        <  1:2567434/147‑1 (MQ=255)
                                                                gtttggttTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGTTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                                                       <  1:4803143/148‑1 (MQ=255)
                                                                ggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                                                       <  1:6552254/150‑1 (MQ=255)
                                                                 gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAAACCACACCCACAAGGGCGTGGGGCCCCATCGAGCACCAAATTTGTGAGATATCTTCTTATTTCGTGTAACagaaga                                                                                      >  2:3735443/1‑150 (MQ=21)
                                                                         tctccctgtccgTGTGGGATAGGGCGGGGTTGGGGGCCTCAGCGCGCACGTTCACCCTCTCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                                              <  2:3582767/150‑1 (MQ=255)
                                                                             cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                          <  2:1456709/150‑1 (MQ=255)
                                                                                ctccgtGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                                                       <  1:6178150/150‑1 (MQ=255)
                                                                                   cctGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                                    <  2:4533980/150‑1 (MQ=255)
                                                                                         ggAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                              <  1:1896138/150‑1 (MQ=255)
                                                                                         ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                              <  1:1069103/150‑1 (MQ=255)
                                                                                         ggAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                              <  2:901485/150‑1 (MQ=255)
                                                                                                  cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                                     <  1:950443/150‑1 (MQ=255)
                                                                                                   ggggggggggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:3292240/150‑1 (MQ=255)
                                                                                                   ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:8707281/150‑1 (MQ=255)
                                                                                                   ggggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:3338625/150‑1 (MQ=255)
                                                                                                   ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:694673/150‑1 (MQ=255)
                                                                                                   ggggTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  2:1041365/150‑1 (MQ=255)
                                                                                                   ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                                    <  1:8363227/150‑1 (MQ=255)
                                                                                                        gAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                                               <  2:4241837/150‑1 (MQ=255)
                                                                                                         aGGGCATCAGCGCGCACGTTGACCCACAACCGAAACCACACCCACAAAGGGGAGGGGACCCAACGAGCACAAAATTTTAGAGATAAAATAATAAGTCTTTTTACAAAAAAAAAGAAAAAACGGCGTTTTCGCCGGAATAATTCGGGGaaa                                              >  2:3802493/1‑150 (MQ=255)
                                                                                                           ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                            <  1:3009625/150‑1 (MQ=255)
                                                                                                           ggCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCACAAGGGAGAGGGGACCCAACCAGCACAAAAAATTAGATAAAAAATAATAAGTCATTTAAAAAAAAAACACAAAAAAGGGCTTTTTCGCCAGAAAAATTCGGAGaaaaa                                            >  2:9345931/1‑150 (MQ=255)
                                                                                                            gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                           <  1:229463/150‑1 (MQ=255)
                                                                                                            gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                           <  2:806771/150‑1 (MQ=255)
                                                                                                             cctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                                          <  1:3289124/148‑1 (MQ=255)
                                                                                                              ctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                         <  2:300753/149‑1 (MQ=255)
                                                                                                              aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                         <  2:3037034/150‑1 (MQ=255)
                                                                                                                  gctatggcGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGc                                     <  1:5449830/143‑1 (MQ=255)
                                                                                                                     cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCNTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                                  <  1:3275969/150‑1 (MQ=255)
                                                                                                                      gcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCggg                                 <  2:541593/150‑1 (MQ=255)
                                                                                                                           gtcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  2:525791/149‑1 (MQ=255)
                                                                                                                           ctcaccctcacccGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                            <  1:1555745/149‑1 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGa                      >  2:5075352/1‑150 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:3192781/1‑150 (MQ=255)
                                                                                                                                 ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                      >  1:7512357/1‑150 (MQ=255)
                                                                                                                                  tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAg                     >  2:2236743/1‑150 (MQ=255)
                                                                                                                                  tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg                     >  2:616667/1‑150 (MQ=255)
                                                                                                                                   caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc                    <  2:7702615/150‑1 (MQ=255)
                                                                                                                                     ccatAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCa                  >  2:5149452/4‑150 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  1:6745143/1‑150 (MQ=255)
                                                                                                                                      cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                 >  1:3245331/1‑150 (MQ=255)
                                                                                                                                       ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt                <  2:4968993/150‑1 (MQ=255)
                                                                                                                                       ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt                >  2:401599/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               >  2:5423718/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               >  1:9020294/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt               <  1:4358864/150‑1 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGAGGAGAAAAGGCATGGCCGGGAATGACCGCGAGCCCAtt               >  1:7235578/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATGGCGGAGAACCAGCATGGCCGGGACTGACCGCGCGCCAAtt               >  1:5267560/1‑150 (MQ=255)
                                                                                                                                         aaCCCTCTCCCTCCAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              <  1:1823781/150‑1 (MQ=255)
                                                                                                                                         aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg              >  2:2940834/1‑150 (MQ=255)
                                                                                                                                           cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgca            >  1:4879972/1‑148 (MQ=255)
                                                                                                                                           cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  2:3465323/150‑1 (MQ=255)
                                                                                                                                           cccTCTACATCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGACTATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa            <  2:8388119/150‑1 (MQ=255)
                                                                                                                                            ccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           >  2:8715333/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc           <  2:8130761/150‑1 (MQ=255)
                                                                                                                                             ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  2:6596451/1‑150 (MQ=255)
                                                                                                                                             ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg          >  1:5668620/1‑150 (MQ=255)
                                                                                                                                                tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaa       >  1:2834470/1‑150 (MQ=255)
                                                                                                                                                 cccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaaa      >  1:6433092/1‑150 (MQ=255)
                                                                                                                                                  ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTCCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAg     <  1:2466138/150‑1 (MQ=255)
                                                                                                                                                   cTCACGGGAGAGGGGACCGCTCGAGCACAAATCTTGAGAGATATCTTCTTCTGTCTTGTAACAGCAGAACAGAAAATCGGGCGTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    <  1:7915820/150‑1 (MQ=255)
                                                                                                                                                   cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  1:9194796/1‑150 (MQ=255)
                                                                                                                                                   cTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg    >  2:5313414/1‑150 (MQ=255)
                                                                                                                                                    tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAggg   >  2:222010/1‑150 (MQ=255)
                                                                                                                                                     cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa  <  1:8183221/150‑1 (MQ=255)
                                                                                                                                                     cAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGa  <  2:1706411/150‑1 (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: