Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,967,130 Δ1 bp coding (1305/2247 nt) ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,1270T.95.7% 189.9 / 2.4 46coding (1308/2247 nt)ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (20/24);  minor base C (2/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGG  >  NC_000913/2966831‑2967419
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                    
aagaGCTCCGGCTGATCGAGAGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGATGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACttta                                                                                                                                                                                                                                                                                                   <  1:800132/300‑2 (MQ=255)
     ctcCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGc                                                                                                                                                                                                                                                                                             <  1:736146/300‑1 (MQ=255)
          gCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTCTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACTCTCACCg                                                                                                                                                                                                                                                                                       >  2:158535/1‑301 (MQ=255)
                      aaaTGCGAATCCCACGCCAACCCAGGCATGGAGTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACGTTGCAGCATGAATGGGAGTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTgc                                                                                                                                                                                                                                                                           <  1:158534/300‑1 (MQ=255)
                          ccGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCgg                                                                                                                                                                                                                                                                       <  1:221351/300‑1 (MQ=255)
                            aaaTCCCACGCCAACCCAGGAATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGTGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGt                                                                                                                                                                                                                                                                     <  1:677291/300‑1 (MQ=255)
                              aTCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTaa                                                                                                                                                                                                                                                                   <  2:227420/301‑1 (MQ=255)
                                cccACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAAt                                                                                                                                                                                                                                                                 >  2:178646/1‑301 (MQ=255)
                                                         tctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTTACGTCATCTTCCGCCAGACgg                                                                                                                                                                                                                                          >  1:424205/1‑299 (MQ=255)
                                                                 gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGCCTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCaa                                                                                                                                                                                                                                >  1:685051/1‑301 (MQ=255)
                                                                       cATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATctct                                                                                                                                                                                                                           <  2:685056/300‑1 (MQ=255)
                                                                           tAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACt                                                                                                                                                                                                                      <  2:612711/301‑1 (MQ=255)
                                                                                                 cGACATCCAGCTTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATTAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTgc                                                                                                                                                                                                <  1:730267/301‑1 (MQ=255)
                                                                                                             tACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCGGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGc                                                                                                                                                                                    <  2:424209/301‑1 (MQ=255)
                                                                                                              aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGc                                                                                                                                                                                    >  2:615671/1‑300 (MQ=255)
                                                                                                               cGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGTATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCt                                                                                                                                                                                   >  2:148447/1‑300 (MQ=255)
                                                                                                                    ggTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGAt                                                                                                                                                                             >  2:106499/1‑301 (MQ=255)
                                                                                                                      tGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATc                                                                                                                                                                            <  1:178645/300‑1 (MQ=255)
                                                                                                                                  aaTTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTGGCAGCATGAATGGGATTTCAGTGCGCTAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACAGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCa                                                                                                                                                               <  1:148447/300‑1 (MQ=255)
                                                                                                                                            tgtgctcttccgatcttACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAg                                                                                                                                                      <  2:805175/284‑1 (MQ=255)
                                                                                                                                             tGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGc                                                                                                                                                     >  2:466517/1‑300 (MQ=255)
                                                                                                                                                      cccTGATACTGCGCCACCTGTTCTTCTTCCGACGAGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGAt                                                                                                                                           >  2:432599/1‑301 (MQ=255)
                                                                                                                                                       ccTGATACTGCGCCACCTGTTCTTCTTCCGAAGTGAAACCACTTTGCAGCATGAATGGGAGTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCGTTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCTGATCGACCAGCAATTCACCGCGATAGCCATCAACGATc                                                                                                                                          <  1:615666/301‑1 (MQ=255)
                                                                                                                                                            tACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAgagatcggaagagcgtcg                                                                                                                                     >  1:805167/1‑284 (MQ=255)
                                                                                                                                                              cTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTcc                                                                                                                                   <  1:466513/301‑1 (MQ=255)
                                                                                                                                                                     ccgatctcttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATg                                                                                                                            <  1:18759/295‑1 (MQ=255)
                                                                                                                                                                         tttctttttcCGACTAGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCCTCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATgcagc                                                                                                                        <  1:432595/300‑1 (MQ=255)
                                                                                                                                                                          ttcttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATgcagca                                                                                                                       >  2:198166/1‑301 (MQ=255)
                                                                                                                                                                            cttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATgagatcgg                                                                                                                     >  2:18759/1‑294 (MQ=255)
                                                                                                                                                                                 tcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAgt                                                                                                                <  1:548339/301‑2 (MQ=255)
                                                                                                                                                                                 tcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAgg                                                                                                                >  2:548344/1‑301 (MQ=255)
                                                                                                                                                                                                    gcagcaTGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCAt                                                                                              >  2:395516/1‑300 (MQ=255)
                                                                                                                                                                                                            aaTGGGATTTCAGTGCGATAAATACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTTCGCACGACGGCTAAGCGCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCTGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCGACGATGCAGCACCGAAGGCTGAATATCCGGGCCCATCACGGGAgg                                                                                     >  1:647072/1‑301 (MQ=255)
                                                                                                                                                                                                                          gCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTTACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAATCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCCATATCCATCACCGATCAGCGTCCTACGATGCAGCACCTAAGGCTGAATATCCGAGCCCCTCACGGTTGGTCTCCCCAgcgcc                                                                        >  1:576237/1‑299 (MQ=255)
                                                                                                                                                                                                                           cGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATAGCCGCGCCCATCACGGTAGGGATCCCCAGCGAAc                                                                       <  1:198165/300‑1 (MQ=255)
                                                                                                                                                                                                                                     cGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGAt                                                             >  2:734972/1‑300 (MQ=255)
                                                                                                                                                                                                                                              cAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAGTCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCCTgg                                                   >  1:405832/1‑301 (MQ=255)
                                                                                                                                                                                                                                              aaaTACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATgg                                                   <  1:395512/300‑1 (MQ=255)
                                                                                                                                                                                                                                               cggtgactggagttcagacgtgtgctcttCCGATCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATggg                                                  <  2:770445/274‑1 (MQ=255)
                                                                                                                                                                                                                                               aaTACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACAATGATCGCAGCATggg                                                  >  2:232754/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         cccAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGCCTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTAACGTCTCCGGATCGACCAGCAATTCACCGCGATAGCTATCAACGATCATCGTCCTACGATGCAGAACCGAAGGCAGAATATCTGCACCCAACACGGTAGGGATCCCCAGCGCACGTAACATGATCGCAGCATGCGCGGTGGcgt                                         >  1:788785/1‑298 (MQ=255)
                                                                                                                                                                                                                                                              tttttCTTCATGTTCCGGTCTTAAACCAGCATTGAGCATGACGTTATACGCTCACCGCTTTGTAACTGCGCGGGTAAATTGACGTCATCGTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGTACGTACCATGATCGCAGCATGGGAGTTGGCTGCGcc                                     <  2:647077/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tcATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACAATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCAc                            <  1:232750/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tcATGTTCCGGGCTTAAACCAGCATTGAGAATGACCTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCCCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCAATCAACGAACAGCGTCCTACGATGCAGCACCGAAGGCTGAATAGCCCCGCCCATCACGCTAGGGATCCCCAGCGCACGTACCACGATCGCAGCATGGGAGTAGACTGCGCCATCTCGCa                             >  2:258257/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                gCTTAAACCGGCATTTAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCg                  <  2:788792/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             ttGAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTg    <  1:734966/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               gAGCATGACTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCTTggg  >  1:588522/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                    
AAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGACTTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTGGG  >  NC_000913/2966831‑2967419

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