Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,518,501 C→T intergenic (‑15/+42) aroB ← / ← aroK 3‑dehydroquinate synthase/shikimate kinase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,518,5010CT95.4% 209.4 / NA 67intergenic (‑15/+42)aroB/aroK3‑dehydroquinate synthase/shikimate kinase 1
Reads supporting (aligned to +/- strand):  ref base C (0/2);  major base T (31/31);  minor base G (0/3);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCTTTTATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGA  >  NC_000913/3518255‑3518745
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ggCTTTTATACTGCCCGCCGTCAGGGAGGATAACGCTATCGAAGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGAGCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTAAGGGAAAAGAGAGGAACGACAATCCTCTCCATAACGCGACATCCACTTTaa                                                                                                                                                                                                                                                  <  1:334291/251‑1 (MQ=255)
 gCTTTTTTCCTGCCCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAAt                                                                                                                                                                                                                                                 <  1:464943/251‑1 (MQ=255)
 gCTTTTATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAAt                                                                                                                                                                                                                                                 <  1:167478/251‑1 (MQ=255)
  cTTTTATACTGCTCGCCGTAAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACACACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCCCTTTAAtt                                                                                                                                                                                                                                                >  1:429286/1‑251 (MQ=255)
        taCTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCACGTACGCCGCGGACCTTATCGAGATACAGAGGCGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGCTGCGTTGGTAATTGGGTAACTACGTTCTCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGt                                                                                                                                                                                                                                           >  2:512084/1‑250 (MQ=255)
            gCTCGCCTTCAGGGGGGATAACGCTATCGCCTTTAACCGCCGCCTGTTCAAGTACGCCGCGGACCTTATCGATATACAGAGGAGCCAGGGTTTCGTTGGTTACCAACATGACCGGCTCGCCCGAGTTCAGCGGTAAGAATGAAGCTGTTTCATTAAACAAACCAGATGCGACTGTTATTGGGTAACTACTTTCCCCGATAGTAACGACAATCCTCTCAATAGCGCGACATCCACCTTAATTACTGTTccc                                                                                                                                                                                                                                       <  2:483768/250‑1 (MQ=255)
                gCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGGAATTGGGTAACTACGTTCCCCGAGAGTACCGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGc                                                                                                                                                                                                                                     >  1:324053/1‑248 (MQ=255)
                 ccGTCAGGGGGGAGAACGGTATCGAGGTGACCACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGGGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACGTTAATTACTGTACCCGCAGAc                                                                                                                                                                                                                                 <  2:494532/251‑1 (MQ=255)
                 ccGTCAGGGAGGATAACGCTATCGACGTTAACACACGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGGGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAAACCTCTCCATAACGCGACATCCACTTTAATTAGTGTACCCGCAGAc                                                                                                                                                                                                                                 >  2:378509/1‑251 (MQ=255)
                   gtcagaacggctaaagCTATCGTCGTTAAAATCCGCCTATTCACGCACGCCGCGGACATTATCGAGACACAAAAGCGCCCGGCATTCGTAGCTGACCAACAACAACTGCTCGCCCGATTTCAGCGGTGAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTACCCGAGAGTAACGACAATCCGCTCCATAACGCGACATCCCCTTTACTTACTGTCCCCACAGACg                                                                                                                                                                                                                                >  2:536016/16‑250 (MQ=255)
                     cAGGGAGGATAACGCTATCGACGTTAACCCCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTCCCCGCAGACGAg                                                                                                                                                                                                                              <  1:461492/250‑1 (MQ=255)
                       gggAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAgtgt                                                                                                                                                                                                                           <  1:146440/251‑1 (MQ=255)
                           ggaTAACGATATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATAGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGGGAGtatat                                                                                                                                                                                                                       >  2:404968/1‑251 (MQ=255)
                                 cGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGcc                                                                                                                                                                                                                 <  1:7578/251‑1 (MQ=255)
                                  gCTATCGACTTAAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCa                                                                                                                                                                                                                <  2:79735/251‑1 (MQ=255)
                                  gCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCACAATGACCTGCTCGCCCGATTTCAGCGGTAAGCATGCAGCTGGTTCATTACACAAACCAGATGCGATGGTAATTGGGCAACTACGTTCCCCGAGAGTAACGTCAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGcc                                                                                                                                                                                                                 >  1:574312/1‑250 (MQ=255)
                                               acacCCGCCTGTTCAAGTACGCCGCGTACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGCCCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGATAGTAACGACAATCCTCTCCATAACTCGACATCCACGTTAAGTACTGTTCCCGCAGACGAGTGTATATAAAGCCAGAGTTAGTTGCtg                                                                                                                                                                                                   <  2:40308/251‑2 (MQ=255)
                                               acacCCGCCTGTTCAAGTACGCCGCGGACCTGATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAACCCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCtt                                                                                                                                                                                                   <  1:313369/251‑1 (MQ=255)
                                                  cccGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGCGCCACAGATTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTACACAAACCAGATGCGATGGTAATTGGGTAACTAAGTTCCACGAGAGTAAAGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAACGCCAGAATTAGTTACTTTcc                                                                                                                                                                                                >  1:270443/1‑251 (MQ=255)
                                                            cAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGaa                                                                                                                                                                                      <  1:294095/251‑1 (MQ=255)
                                                              aGTACGCCGCGGACCTTAGCGAGATGCAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACGCTCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGaata                                                                                                                                                                                    <  1:552082/251‑1 (MQ=255)
                                                                  cgccgcGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGATTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCt                                                                                                                                                                                <  1:510754/251‑1 (MQ=255)
                                                                  cgccgcGGACATTGTCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCCACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCt                                                                                                                                                                                >  2:83471/1‑251 (MQ=255)
                                                                              tATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATaaa                                                                                                                                                                                                                    <  2:508958/203‑1 (MQ=255)
                                                                              tATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATaaa                                                                                                                                                                                                                    >  1:508958/1‑203 (MQ=255)
                                                                                          gaggagCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGAGTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGGGTATATAAAGCCAGAATTAGTGGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTGTAGCGCttt                                                                                                                                                        <  2:25823/251‑1 (MQ=255)
                                                                                             gagCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGat                                                                                                                                                     <  1:231747/251‑1 (MQ=255)
                                                                                                                gACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCTAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCTCTTTGTTCATCAGGACGAATGGTCa                                                                                                                                   <  2:90194/250‑1 (MQ=255)
                                                                                                                   cAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTc                                                                                                                               <  1:340036/251‑1 (MQ=255)
                                                                                                                    aaCATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGCATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCg                                                                                                                              >  1:396778/1‑251 (MQ=255)
                                                                                                                     aCATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATGCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGAGTATCAGTTTGATTGGTCACGTCgg                                                                                                                             <  2:491271/251‑1 (MQ=255)
                                                                                                                       aTGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATTGTCACGTCGGCa                                                                                                                           <  1:179765/251‑1 (MQ=255)
                                                                                                                                tCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCAGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGACAATCTCTTcc                                                                                                                  >  2:570087/1‑250 (MQ=255)
                                                                                                                                  gCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGAAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCAt                                                                                                                 >  1:335966/1‑250 (MQ=255)
                                                                                                                                           cAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAAATGCGATGGTAATTGGGTAACTACGTTCCACGAAAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCCTGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTAGGCAATCTCTTCATACAGCGGAtt                                                                                                       >  2:429482/1‑251 (MQ=255)
                                                                                                                                             gCGGTAAGAATGAAGCTGGTTCCTTACACAAACCAGATGAGCTGGTAATTGGGAAACTACGTTCCCCGAGAGTCAAGACAATACTCTCCATAACGCGACATCAACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGaataa                                                                                                                                                                                   >  2:519725/1‑173 (MQ=255)
                                                                                                                                             gCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGaataa                                                                                                                                                                                   <  1:519725/173‑1 (MQ=255)
                                                                                                                                               ggTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCTAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGt                                                                                                   <  1:381598/251‑1 (MQ=255)
                                                                                                                                                gTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGtt                                                                                                  <  2:224702/251‑1 (MQ=255)
                                                                                                                                                                 tCATTAGACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCACGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGCGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCCTCAGTACGAATGGTCACGTCGGCAATCTCTGCATACAGCGGATTGCGTTCATTGGCCAACGCTTcc                                                                                 >  2:291630/1‑251 (MQ=255)
                                                                                                                                                                          aaacCCGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGACGGCACTCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTc                                                                         >  2:483676/1‑250 (MQ=255)
                                                                                                                                                                          aaacCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTc                                                                         >  1:523083/1‑250 (MQ=255)
                                                                                                                                                                           aacCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCAc                                                                       >  1:492010/1‑251 (MQ=255)
                                                                                                                                                                                     gATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCAc                                                                                                                                  <  2:526641/182‑1 (MQ=255)
                                                                                                                                                                                     gATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCAc                                                                                                                                  >  1:526641/1‑182 (MQ=255)
                                                                                                                                                                                     gATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGgtgt                                                             >  1:422136/1‑251 (MQ=255)
                                                                                                                                                                                           aaTTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAAGTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGGTCATTGGCCAACGCTGCCAGAACTTCACGCTGCGTTGGTTCAAc                                                       <  2:38929/251‑1 (MQ=255)
                                                                                                                                                                                              tGGGTAACTACGTTCCCCGAGAGTATCGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGCGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATCATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACATCGGCAATCTCTTCATACAGCGGTTTGCGTTCATTGGCCAACGCTTCCAGTACTTCACGCGGCGGTGTTTAAACGTg                                                    >  2:268369/1‑251 (MQ=255)
                                                                                                                                                                                                   aaCTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCCCTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTgcagc                                                <  1:257842/250‑1 (MQ=255)
                                                                                                                                                                                                   aaCTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGGGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCAGCAGTACGAATGGTCACGTCGGCAATCGCTTCATACAGCGGATTGCGGTCAATGGCCAACGCTGCCAGAACTACACGCGGCGGTGTTTCAGCGTgcagc                                                <  2:302112/250‑1 (MQ=255)
                                                                                                                                                                                                       acggtCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCTGTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGGCAGAATTAGTTGCTTTCCAGCATGTGAATAAGCTGGTTTGCAACCACTTTAGCGCTTTGAGCATCAGTTCGAATGGTCACGTCTGCAATCTCTTCATACAGCGGATGGCGTTCATTGGCCAACGCTTACAGAACTTCACGCGGCGGTGTTTCAACGTGCAGTAAc                                             <  2:179220/245‑1 (MQ=255)
                                                                                                                                                                                                       aCGTACCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACGTTAATTTCTGTACCTGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACGGCGGATTGCGTTCATTTGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAAc                                             <  1:268369/249‑1 (MQ=255)
                                                                                                                                                                                                         gTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGc                                         <  2:450158/251‑1 (MQ=255)
                                                                                                                                                                                                            ccccGAGAGTAACGACAATCCTCTCCATAACGAGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGACTGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGtt                                      >  1:340931/1‑251 (MQ=255)
                                                                                                                                                                                                                gagaGTAACGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCGGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTAATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGAGCAGCAACGGGCGtttttt                                  >  1:376747/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaCGACAATCCTCTCAATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATACAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCCCGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACg                            >  2:41393/1‑251 (MQ=255)
                                                                                                                                                                                                                        cGACAATCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGAAGCAACGGGCGTTTTTTATCACGCt                          >  1:39762/1‑251 (MQ=255)
                                                                                                                                                                                                                             aTCCTCTCCATAACGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCtgcgtg                     >  1:314217/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaCGCGTCATCCACTTTAAGTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACTAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGTCGGTGTTTCAACGTTCAGCAATGGGCGTTTGTTATCACTCTGCGTGTGTGCAAGtt           <  2:563416/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aaCGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATCAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGGTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCCATCTCTTCATACAGCGGATTGCGTTCATTGGCCCACGCTTCCAGAACTTCACGCTGCGGTGTTTCAACGTGCATCAACGGGCGTTTTTTATCCCGCTGCGTGCGTGCAAGttt          >  2:97073/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aaCGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTg          <  1:388261/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aaCGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTg          >  1:123417/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaCGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTg          >  2:434577/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aCGCGACATCCGCCTTAATTACTGTCCCCTCTGACGAGTGTAGATATAGTCAGAATTAGTTGCTTTCCTGCCCGGCAATAATCTGGTTTGCATCCACCTTAGCGCTTTGAGCATCAGTACGAATGTTCACGTCTGCAATCTCTTCATACATCAGATTGCGGTCATTGTCCAACGCTTGCAGCACTTCTCCCGGCCGTGTTTCAGCGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGc         <  1:377304/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aCGCGACATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGc         >  1:180660/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cATCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGGCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCg   <  1:383153/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aTCCACTTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGa  <  2:461717/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GGCTTTTATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGA  >  NC_000913/3518255‑3518745

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: