Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,702,839 T→C 8.6% Q10Q (CAA→CAG dppF ← dipeptide/heme ABC transporter ATPas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,702,8390TC 8.6% 143.9 / 2.5 53Q10Q (CAA→CAGdppFdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base T (27/21);  new base C (3/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTAAGGGTTCTGGAAGACGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGTATACGGGTGACGCGGCGCATGGAAGATGGCGTGCGCA  >  NC_000913/3702542‑3703138
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
ttAAGGGTTCTGGAAGACGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgt                                                                                                                                                                                                                                                                                                          <  2:212927/300‑1 (MQ=255)
tcagtcgtgtgctcttccGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgt                                                                                                                                                                                                                                                                                                          <  2:561136/284‑1 (MQ=255)
              aGACGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCctc                                                                                                                                                                                                                                                                                            <  2:408907/301‑1 (MQ=255)
                  gATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGGAGCCCCTTGc                                                                                                                                                                                                                                                                                        >  1:240299/1‑301 (MQ=255)
                       ggATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCCCGTTCCAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGgttgtagatcggaagaacgtcgggtagg                                                                                                                                                                                                                                                                                   >  1:560422/1‑278 (MQ=255)
                       cgATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCCGATCCTGCCCCTGGTAATACAGCTCGCCCCCGGTGGGCATTTCAATCCTCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATCAGCCCACTACTACCAGCGTTTACCACCGTACAAGAATAAACGAACGGCCATACAGCGCTTTGACACGACGTTCACGCGGGACCATGCCTTCTCTAACCGGAAAATGTTTTTCCAGGTCGATACCCTAAAAAAGCGGTTATTGCAGGTAGGCCTCTTGCGTACt                                                                                                                                                                                                                                                                                   >  1:463260/2‑301 (MQ=255)
                                   gACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGc                                                                                                                                                                                                                                                                       <  2:240566/301‑1 (MQ=255)
                                      gCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTc                                                                                                                                                                                                                                                                    <  2:629465/301‑1 (MQ=255)
                                          tttCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCgg                                                                                                                                                                                                                                                                <  2:334909/300‑1 (MQ=255)
                                               gcgcCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGATTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAATGTCGGCCTCCAGGaatca                                                                                                                                                                                                                                                           >  2:79954/1‑301 (MQ=255)
                                                  ccTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGcaacaacc                                                                                                                                                                                                                                                        >  1:423306/1‑298 (MQ=255)
                                                   gggcaggatctGCTTAAGCACGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCAtcatca                                                                                                                                                                                                                                                        >  2:85232/11‑301 (MQ=255)
                                                       ggATTGTGCTTAAGCAGCTTCTTCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCCTCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTGACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTGTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTgg                                                                                                                                                                                                                                                   <  2:463824/301‑1 (MQ=255)
                                                                   aGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCCATCATCGTCAGCACCCGACCGAGGGTCGCTTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTCGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCCGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAACCTGCAACAGCAGCTGTTGCAGGGCGGCCCCTTCCGTACAAATCGTTTCCACCTCCCCGCACCACCACGTGGGTAATGGCAttt                                                                                                                                                                                                                                       >  1:324621/1‑301 (MQ=255)
                                                                        gATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTgtaa                                                                                                                                                                                                                                   >  2:160772/1‑298 (MQ=255)
                                                                          tCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGTTTTACCGCAACCAGATTCGCCCACTACTGCCAGCGTTTTCCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCACCCGTTCCGGCGCGAACTTGCATTTCTTCACCGGATTATGTTTTTTCAGGTCGATAGCCTGCACCAGCTGTTGTTGCAGGGTCGCCTCTTGCATACTCATAGTCTCGGCCTCCCGCAATCAACCAGTGGGACATGGACTTTCGACTg                                                                                                                                                                                                                                 >  1:609554/1‑301 (MQ=255)
                                                                                 ccTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCg                                                                                                                                                                                                                          >  2:112688/1‑301 (MQ=255)
                                                                                   atgTAATACAGCGTGCCACCGGTGGGTATTTCAATCATCGTCAGCAACCTACCGAGGGTCGATGTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTGAAAAGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTATTGCGAACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTc                                                                                                                                                                                                                        <  1:160621/299‑1 (MQ=255)
                                                                                         tACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAg                                                                                                                                                                                                                  >  2:464641/1‑301 (MQ=255)
                                                                                         tACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAg                                                                                                                                                                                                                  >  1:234966/1‑301 (MQ=255)
                                                                                                           ggCATTTCAATCATCGTCAGCAACCGACCGAGGGCCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttc                                                                                                                                                                                                >  2:269156/1‑301 (MQ=255)
                                                                                                             aagttgattcatcGTCAGTAACCGGCCGGGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCGCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttctt                                                                                                                                                                                              <  1:254870/293‑1 (MQ=255)
                                                                                                              aTTTCAATCATCGTCAGCAACCGACCGAGTGGCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACGGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGttcttc                                                                                                                                                                                             <  1:477993/301‑1 (MQ=255)
                                                                                                                              taaatcgaccctcggTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGCTAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTTACGCCCGTCAGCGAGCATATTCACAGCCAGTTCTTCAGCGCCACAACAGTcc                                                                                                                                                                             >  1:474012/14‑300 (MQ=255)
                                                                                                                                   tgaccctcggTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGc                                                                                                                                                                        <  2:474613/293‑1 (MQ=255)
                                                                                                                                             tCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAg                                                                                                                                                              <  2:423869/301‑1 (MQ=255)
                                                                                                                                                   tACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGAGGCAGCGCggg                                                                                                                                                        >  2:225176/1‑301 (MQ=255)
                                                                                                                                                     ccGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGAAATCTGTCAGTGGAATAGGGGCAGCGCGGGtt                                                                                                                                                      >  2:284320/1‑301 (MQ=255)
                                                                                                                                                            cAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAgcagg                                                                                                                                               <  2:571650/301‑1 (MQ=255)
                                                                                                                                                                   gCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGt                                                                                                                                        >  1:135594/1‑301 (MQ=255)
                                                                                                                                                                                            ctcttccgttctaacgatacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACg                                                                                                              <  1:105109/289‑1 (MQ=255)
                                                                                                                                                                                                  gtgttaatcgtagtgctttcctgttcTTTAACTAGACTTTCCGGCGCGAACTTGCCTTTCTTCACCTGAGAATTTTTTTTCAGTTCGATAGCCTGCAACAGCGGTTGTTGCAGGTTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACAcc                                                                                                         <  1:250745/276‑1 (MQ=255)
                                                                                                                                                                                                  ccgttctacgaaacgCCATCCAGCTCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACAcc                                                                                                         <  1:54875/294‑1 (MQ=255)
                                                                                                                                                                                                        aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTAGGGCGGTCGTACTTGCCGGGAACGACCCCTGGCaa                                                                                                   >  1:578409/1‑301 (MQ=255)
                                                                                                                                                                                                        aacgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGagatcggaagaa                                                                                                  >  2:105233/1‑290 (MQ=255)
                                                                                                                                                                                                         acgaaacgCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACAccagatcgg                                                                                                  >  2:54943/1‑294 (MQ=255)
                                                                                                                                                                                                             aacgCCATCCAGCGCTTTAACCAGATGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAacgacg                                                                                              <  1:112564/301‑1 (MQ=255)
                                                                                                                                                                                                                gCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCa                                                                                           >  2:471659/1‑301 (MQ=255)
                                                                                                                                                                                                                     ccAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGt                                                                                      <  2:235232/301‑1 (MQ=255)
                                                                                                                                                                                                                            tttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGt                                                                               <  2:135719/301‑1 (MQ=255)
                                                                                                                                                                                                                              tAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCCGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGCCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGCCGCCAGACGTTCTTTCTcc                                                                             >  1:556952/1‑301 (MQ=255)
                                                                                                                                                                                                                              tAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTAGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTcc                                                                             >  2:149073/1‑301 (MQ=255)
                                                                                                                                                                                                                                  cgggcggtcgtgctcgAACATGCTTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAAGAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTGCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAg                                                                         <  1:148922/287‑1 (MQ=255)
                                                                                                                                                                                                                                      cGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGAAGACAGACGTTCTTTGTCCTGAGCaaa                                                                     >  1:564118/1‑301 (MQ=255)
                                                                                                                                                                                                                                           gggtgcgattgtGCCTTTCTGCACTGGATAATGTGTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTg                                                                <  1:268858/290‑1 (MQ=255)
                                                                                                                                                                                                                                                 gAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCACGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCCGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTAGCAACGACGCCAGAAGTTCTTTGTCCTGAGCAAATTCTGGCAgcg                                                          >  1:39547/1‑301 (MQ=255)
                                                                                                                                                                                                                                                   aCATGCCTTTTTTTCCCGGCTATTTTTTTTTCAGTTCTCTAGCGTTCAGTATCGGTGGTTGCATGGTGGCCTCTTGCGTATTCATTGTTGCGGCCTCCAGGCATCATCAAGTGGGTAATGGCATTTGTACTGACGCCGGTCAGCGGGCATATTCTGCGCCGGTTCTTCAGCGTGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGTCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGcgca                                                        <  1:224933/301‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGcgcacgca                                                    >  2:540996/1‑301 (MQ=255)
                                                                                                                                                                                                                                                               tctttccctacccgacgctcttccgatcgcccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCTTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGcc                                            <  1:476674/271‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tttcagacgtgtgctcttccgatctcccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAg                                        <  2:153776/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGt                                       >  2:94427/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             tttttCGGGTCGGTGGCTTGCAACGGCGGTTGTTGCAGGGTGGCCTCTTGCGTAGTCATAGGGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGTATACGGGTg                              <  1:471058/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     gTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACATTCTTTGTCCTGAGCCAATTCTGCCCGCCCCCGAAGCAATCCCCGAGTATACGGGtaccccggcc                      >  1:168319/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                                            cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCAGatcggaagagcacacgactgaactcca               >  2:477275/2‑275 (MQ=255)
                                                                                                                                                                                                                                                                                            cccTGCAACAACCGCTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAAAGCCTGAgagatcggaagagcgtcgtgtagggg               >  1:153625/2‑276 (MQ=255)
                                                                                                                                                                                                                                                                                                      gCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGTATACGGGTGACGCGGCGCATGGAAGATGGCGTgc     >  2:358926/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                         tttgttgCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGTATACGGGTGACGCGGCGCATGGAAGATGGCGTGCGCa  <  2:564832/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
GTAAGGGTTCTGGAAGACGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGTATACGGGTGACGCGGCGCATGGAAGATGGCGTGCGCA  >  NC_000913/3702542‑3703138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: