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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC15.3% 195.3 / 10.8 116intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (27/69);  new base C (2/18);  total (29/87)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                      |                                                                                                                                                     
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                       <  2:2719185/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:4142420/150‑1 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  2:3782485/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  2:1984728/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCcgcgcg                                                                                                                                               >  1:2490508/1‑150 (MQ=255)
              atgtgaTGAAAGGGCGTTATCCCAAACTTGTCTGGCGGGACGACCCGGAAAATACTGCACCGACGAGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                                         <  2:2700444/146‑1 (MQ=38)
              aaGAGATGATTGGGCTTTATCCCAAACTTGTGGGGCCGTACGACCGGGAAAATACTGCACCGACGCGACGGACTAAAATGTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGTGGGCATCAGCGCGCACGtt                                                                                                                                         <  2:3866992/150‑1 (MQ=21)
              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCggg                                                                                                                                         >  1:9558/1‑148 (MQ=255)
              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                         >  2:3912156/1‑149 (MQ=255)
                 agaTGACCGGGCGTTATCCCCAACATGTCTGGCGGGACGCCCCGGCACATACTGCACCGAGGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                      <  2:3055815/150‑1 (MQ=37)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggcc                                                                                                                                      >  1:2003930/1‑146 (MQ=255)
                  tattaatgggCGTTATCACAAACTTTTGTGGCCGGACGACCCGGCAAAAACTGCACCGAGGCGCCGTACTAAATTTTTTTCTTTTGGTTGTTTTTCCCCCTCCCCGTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:3785614/143‑1 (MQ=255)
                  gatgatTGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGAGGGACGAAAAAGTTTTCGTACGGTTGGTTTTCGCCCTGTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:893818/150‑1 (MQ=255)
                  gaTGCACGGGCTTTATCCCAAACTGGTCTGGCCGGACGACCCGGAAAATCCTGCACCGCCGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:5310542/150‑1 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccg                                                                                                                                   >  1:1911418/1‑143 (MQ=255)
                      aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                 >  2:3968180/1‑150 (MQ=255)
                       acgggCGTTATCCCAATCTTGTCGGGCCGGACGCCCCGGCAAATACTGCCCCGACGCGAGGGACTAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                <  2:4180123/148‑1 (MQ=255)
                        tgggCTTTAACCAAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTATTGTTAGGGTTGGTTTACCCCCACGCCCTGTGGGTTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                               <  2:2102292/149‑1 (MQ=255)
                        cgggCTTTACCCCAAACAGGTCTGGCCGGACGACCCGCCAAATACTGCCCCGACGCGCCGGACGAAAATTTTTTGTTTTGGTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                               <  2:2109126/149‑1 (MQ=25)
                             gTTACCCCAATCTTGTCGGGCGGGCCGACCGGTAAATGCCGGCCCCGACGGGACGTACTAAAATGTATTGGTAGGGTTGGTTTTCTCCCCCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTatc                                                                                                                          <  2:4210093/150‑3 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGTCCCACACCCaaaac                                                                                                                         >  1:965308/1‑148 (MQ=38)
                                aacccAAACATCACTGCCCGCACGACCCCCAAAAAACTCCACCCACCCCACGCACGAAAAACTATTCGTAAGGTTGCTTTACACCCACACCCTCTGGGACAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                       <  2:3207817/148‑1 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACGAAAAAGTATACTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:2774557/150‑1 (MQ=38)
                                         ttttCTGGCCGGCCCCCCCGCAAAATAGTGCACCGCCGCGAGGTACCAAAATTTATTCTTACTTTTGTTTTCCCCCCCCCCCGTGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  1:1720012/147‑1 (MQ=255)
                                         ttGTCGGGCCGGACGACCCGGAAATTACTGCACCGAGGCGACGTACTAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCCCCCTGGGGTTGAGGGCCGGGTTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:1773130/149‑1 (MQ=255)
                                         atgtggggCCGGACGACCCGGCAAATACGGCACGGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:4376797/144‑1 (MQ=38)
                                         aTGTCGGGCCGGACGACCCGCCAAATACTGCCCCGACGCGCCGTACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGCTCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                              <  2:4376748/150‑1 (MQ=255)
                                          ttgggggCCGGACGCCCCGGAAAATCCTGCACCGCGGCGAGGGACGAAAATGTTTTCGTTGGGTTGGTTTTCGCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  1:2648209/145‑1 (MQ=25)
                                          tGTCGGGCCGGACGACCCGGAAAATACTCCACCTACGCGACGGACTAAAATGTATTCGTATGGTTGGTTTTCCCCCTCTCCCTTTGGGATTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:3893837/150‑1 (MQ=37)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAACCCACAccccccacgg                                                                                                           >  2:641800/1‑143 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCGCACCCGAAACCACGCCcacaccggcg                                                                                                         >  1:3001512/1‑144 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCTCTCCCGCCAGGgag                                                                                                         >  2:792226/1‑150 (MQ=255)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAACCCGCTCCCCCAAgggtg                                                                                                         >  1:4397130/1‑148 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGAGCCCACACCCAAACCCACTCACACAAgggtg                                                                                                         >  1:3001511/1‑148 (MQ=37)
                                                gCCGGACGACCCGGAAAATACTGCACCGCCGCGACGGACGAAAATTTTTTCGTATGGTTGGTTTCCCCCCCCTCCCTGGGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  1:4899699/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                      <  2:5427290/150‑1 (MQ=255)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                      <  2:897620/150‑1 (MQ=255)
                                                 ccGGACGACCCGCAAAATACGCCACGGACGCGCGGGACTAAAATTTTTTCGTCTGGTTGTTTTTACCCCGCGCCCTGTGGGAGCGGGGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                      <  1:5226318/150‑1 (MQ=255)
                                                  cGGCCGCCCCGGCAAATCCGCCCCCGCCCCGCCGGAGGAAAACGTCTTCGTACGGTGGGTTTCCTCCCTCGCCCGGGGGGCGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  2:2922165/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCTCTCCCGCAAGGGAGgggg                                                                                                     >  2:4135861/1‑150 (MQ=38)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGACGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGGGCGCACGGTCACCCTCACCCGCACCCTCTCCCGCCCGGGAGAggg                                                                                                     >  2:4135392/1‑150 (MQ=37)
                                                    gccgccccGGAAAAAAGTGCACCTAGGCGCCGGCCTAAAATTTTTTCTTATGGTTGTTTTTCTCCCTCTCCGGGGGGGTGTGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                   <  1:1139393/145‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCACCCAAACCCCCCCCCCAAAGGGGGTGGGGGcccaac                                                                                             >  2:1710160/1‑146 (MQ=25)
                                                             gCAAAGACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  2:5042921/150‑1 (MQ=255)
                                                               aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACCCCCAAACCCCCGCCCACACGGGGGGGGGGGCCCCTCGAGcc                                                                                        >  2:2603040/1‑149 (MQ=21)
                                                                aaTACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGGTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGTGAGGGGACCGATCGAGcac                                                                                       <  2:3530845/150‑1 (MQ=255)
                                                                     cTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                   <  2:823151/150‑1 (MQ=255)
                                                                      tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAGCCACACCCCAAACCACACCCACCACGGGGAGGGGCCCGCTCGGGCCCaaaat                                                                                  >  1:809734/1‑148 (MQ=21)
                                                                      tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACAACCAAACCCACTACCACCAGGGAGGGGGGCCCGGGAGAGCACAAAta                                                                                  >  2:266111/1‑149 (MQ=21)
                                                                        cACCGACGCGACGGACTAAAATGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  2:543414/150‑1 (MQ=255)
                                                                         aCCGACGCGACGTACCAAAATGTTTTCGTATGTTTGTTTTCCCCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                               <  2:5149604/150‑1 (MQ=38)
                                                                          ccGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:3968700/150‑1 (MQ=255)
                                                                          ccGACGCGACGGACGAAAAATTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:4555868/150‑1 (MQ=38)
                                                                          ccGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:4455943/150‑1 (MQ=255)
                                                                           cgacgcgaggGACGAAAATGTTTTCGTACGGTTGGTGTTCTCCCTCCCCCGGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                             <  1:1238378/150‑1 (MQ=38)
                                                                            gaggcgacgGACGAAAATGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:2662313/147‑1 (MQ=255)
                                                                             acgcgccgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:3929264/150‑1 (MQ=255)
                                                                                            aattATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                            <  1:5466108/147‑1 (MQ=255)
                                                                                              gTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:2209148/150‑1 (MQ=255)
                                                                                                  tCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  2:3026337/150‑1 (MQ=255)
                                                                                                      atgtttggttTCCTCCCTCCCCGTGGGGTGGTGGGCCGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  1:2165767/146‑1 (MQ=255)
                                                                                                      aatgttggttTTCCCCCTCTCCCTGTGGGAGAGGGCCGTCGTGAGGGCATCACCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGCTGACCGATCGAGCACAAATTTAGAGAGATATATTCTTCTGTCTTGTAACaca                                                  <  1:1830041/147‑3 (MQ=38)
                                                                                                       tggttggttTTCCCCCCCTCCGGGGGGGTGCGGGCGGGGTGGTGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  2:9558/149‑1 (MQ=255)
                                                                                                       aGGTTGGTTTTCACCCTCACCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:4788362/150‑1 (MQ=255)
                                                                                                           tggttTTCCCCCTCGCCGTGTGGGTGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  1:304733/150‑1 (MQ=255)
                                                                                                                  ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGAGCCCAAACCCACACCCAACAAGGAGAGGGGCCCGCTCGAGCACAAATGTTGAGAGAAATCATATTCAGTGTTGTAACAGAACAACAGaaaaa                                      >  2:2285069/1‑149 (MQ=16)
                                                                                                                   tccctctccgtGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:1760805/150‑1 (MQ=255)
                                                                                                                   tcccccccccgGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:2030071/150‑1 (MQ=255)
                                                                                                                   cccccccccctGCGGCAGAGGCCCGGGGAGCGCCCCTCAGCGCCCACGTTCACCCACACCCTAACCCTCACCCACAACGGAGAGGGGACCGATCGAGCACAAAATTTGAGAGATATCTTCTTCTGACTTGTAACAGAAGAACAGAAAATc                                     <  2:1208798/143‑1 (MQ=255)
                                                                                                                    ccctctccgtGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:2252226/150‑1 (MQ=255)
                                                                                                                    ccctctccctGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:2252437/150‑1 (MQ=255)
                                                                                                                      ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                  <  2:3479033/150‑1 (MQ=255)
                                                                                                                       tcgcccgGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                                 <  1:3949615/147‑1 (MQ=255)
                                                                                                                           cctGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  1:1216789/150‑1 (MQ=255)
                                                                                                                           ccacagggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAACCCACACCCCCAAGGGGGCGGGGCCCCAAAGAGCCCAAAAGTTTAGAGATAAAATAGTCTGTCTTGTAAAACAATAACAGCAAAAAGGGCtttt                             >  2:1800096/6‑150 (MQ=255)
                                                                                                                             tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                           <  2:3501762/150‑1 (MQ=255)
                                                                                                                             tgggggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                           <  2:5175163/147‑1 (MQ=255)
                                                                                                                               tgGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  1:3444659/150‑1 (MQ=255)
                                                                                                                               ggggAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  2:3931816/149‑1 (MQ=255)
                                                                                                                                gggggggggCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  2:924817/144‑1 (MQ=255)
                                                                                                                                 ggtgtgggCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:1029852/145‑1 (MQ=255)
                                                                                                                                 ggcgagGGCGGGGGTGTGGGCCTCATCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACAGATCGAGCACAAATTTTGAGAGATATCTTCGTCTGTCTTGTAACAGCAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:496026/147‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:129068/150‑1 (MQ=255)
                                                                                                                                 ggAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:496067/150‑1 (MQ=255)
                                                                                                                                  gagcGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                      <  2:4274172/150‑1 (MQ=255)
                                                                                                                                   agagGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                     <  2:5325084/150‑1 (MQ=255)
                                                                                                                                    gtgggCCGGGTGGCGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                    <  2:4831338/148‑1 (MQ=255)
                                                                                                                                     cgggCCGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                   <  1:1324637/149‑1 (MQ=255)
                                                                                                                                       ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAGGGAGTGGGGGCCGCTAGTGCCCAAAAATTTAGTGATATAATATTAAGTGTTGTTAAAGAAGTACACAAAATCGGGCGTTTTCGCCGGAatat                 >  2:2746333/1‑150 (MQ=255)
                                                                                                                                          gggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:4808749/149‑1 (MQ=255)
                                                                                                                                          cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACCGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  2:3770657/150‑1 (MQ=255)
                                                                                                                                           ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:4376681/150‑1 (MQ=255)
                                                                                                                                           ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:2603040/150‑1 (MQ=255)
                                                                                                                                             ggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  1:196616/150‑1 (MQ=255)
                                                                                                                                             ggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCAATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:1523559/150‑1 (MQ=255)
                                                                                                                                                   ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:3491620/150‑1 (MQ=255)
                                                                                                                                                    gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  1:437586/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:3324774/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:446817/150‑1 (MQ=255)
                                                                                                                                                      aTCAGCGCGCACGTTCACCCACAACCAAAACCAAAACCACAAGGGAGAGGGGACCCCTAGAGCCCAAAATTTTAGAGAAATAATATTATGTCTTGTAACAGAAGAACACAAAAACGGGCATTTGCGCCAGAAAATTTCTGATAAAAAGCa  >  2:2284516/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730

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: