Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,702,835 G→A 7.6% L12L (CTG→TTG)  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,8350GA 7.6% 158.5 / 2.8 56L12L (CTG→TTG) dppFdipeptide/heme ABC transporter ATPas
Reads supporting (aligned to +/- strand):  ref base G (28/22);  new base A (2/2);  total (31/25)
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.29e-01

AACCGTAAGGGTTCTGGAAGACGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGTATACGGGTGACGCGGCGCATGGAAGATGGCGTGC  >  NC_000913/3702538‑3703135
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                            
aaCCGTAAGGGTTCTGGAAGACGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCCGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCAATCCCGCGCTTTAACCAGACATTCCGGCGCCAACCTGCCTTTCCTCACCGGAAAATGTTTTTTCAGGTCGATAGCCTTCAACACCgg                                                                                                                                                                                                                                                                                                           >  1:334504/1‑301 (MQ=255)
    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)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                            
AACCGTAAGGGTTCTGGAAGACGATCTGGATTTTCTGCCGACGCAGCTTCTGCGCCTGCGGATCGTGCTTAAGCAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACGACACCTGGCAACGACGCCAGACGTTCTTTGTCCTGAGCAAATTCTGGCAGCGCACGCAGCAATGCCTGAGTATACGGGTGACGCGGCGCATGGAAGATGGCGTGC  >  NC_000913/3702538‑3703135

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