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)
     ctcCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGc                                                                                                                                                                                                                                                                                             <  1:736146/300‑1 (MQ=255)
          gCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTCTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACTCTCACCg                                                                                                                                                                                                                                                                                       >  2:158535/1‑301 (MQ=255)
                      aaaTGCGAATCCCACGCCAACCCAGGCATGGAGTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACGTTGCAGCATGAATGGGAGTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTgc                                                                                                                                                                                                                                                                           <  1:158534/300‑1 (MQ=255)
                          ccGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCgg                                                                                                                                                                                                                                                                       <  1:221351/300‑1 (MQ=255)
                            aaaTCCCACGCCAACCCAGGAATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGTGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGt                                                                                                                                                                                                                                                                     <  1:677291/300‑1 (MQ=255)
                              aTCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTaa                                                                                                                                                                                                                                                                   <  2:227420/301‑1 (MQ=255)
                                cccACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAAt                                                                                                                                                                                                                                                                 >  2:178646/1‑301 (MQ=255)
                                                         tctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTTACGTCATCTTCCGCCAGACgg                                                                                                                                                                                                                                          >  1:424205/1‑299 (MQ=255)
                                                                 gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATG‑CCTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCaa                                                                                                                                                                                                                                >  1:685051/1‑301 (MQ=255)
                                                                       cATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATctct                                                                                                                                                                                                                           <  2:685056/300‑1 (MQ=255)
                                                                           tAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACt                                                                                                                                                                                                                      <  2:612711/301‑1 (MQ=255)
                                                                                                 cGACATCCAGCTTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATTAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTgc                                                                                                                                                                                                <  1:730267/301‑1 (MQ=255)
                                                                                                             tACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCGGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGc                                                                                                                                                                                    <  2:424209/301‑1 (MQ=255)
                                                                                                              aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGc                                                                                                                                                                                    >  2:615671/1‑300 (MQ=255)
                                                                                                               cGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGTATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCt                                                                                                                                                                                   >  2:148447/1‑300 (MQ=255)
                                                                                                                    ggTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGAt                                                                                                                                                                             >  2:106499/1‑301 (MQ=255)
                                                                                                                      tGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATc                                                                                                                                                                            <  1:178645/300‑1 (MQ=255)
                                                                                                                                  aaTTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTGGCAGCATGAATGGGATTTCAGTGCGCTAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACAGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCa                                                                                                                                                               <  1:148447/300‑1 (MQ=255)
                                                                                                                                            tgtgctcttccgatcttACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAg                                                                                                                                                      <  2:805175/284‑1 (MQ=255)
                                                                                                                                             tGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGc                                                                                                                                                     >  2:466517/1‑300 (MQ=255)
                                                                                                                                                      cccTGATACTGCGCCACCTGTTCTTCTTCCGACGAGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGAt                                                                                                                                           >  2:432599/1‑301 (MQ=255)
                                                                                                                                                       ccTGATACTGCGCCACCTGTTCTTCTTCCGAAGTGAAACCACTTTGCAGCATGAATGGGAGTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCGTTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCTGATCGACCAGCAATTCACCGCGATAGCCATCAACGATc                                                                                                                                          <  1:615666/301‑1 (MQ=255)
                                                                                                                                                            tACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAgagatcggaagagcgtcg                                                                                                                                     >  1:805167/1‑284 (MQ=255)
                                                                                                                                                              cTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTcc                                                                                                                                   <  1:466513/301‑1 (MQ=255)
                                                                                                                                                                     ccgatctcttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATg                                                                                                                            <  1:18759/295‑1 (MQ=255)
                                                                                                                                                                         tttctttttcCGACTAGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCCTCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATgcagc                                                                                                                        <  1:432595/300‑1 (MQ=255)
                                                                                                                                                                          ttcttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATgcagca                                                                                                                       >  2:198166/1‑301 (MQ=255)
                                                                                                                                                                            cttcttcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATgagatcgg                                                                                                                     >  2:18759/1‑294 (MQ=255)
                                                                                                                                                                                 tcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAgt                                                                                                                <  1:548339/301‑2 (MQ=255)
                                                                                                                                                                                 tcCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAgg                                                                                                                >  2:548344/1‑301 (MQ=255)
                                                                                                                                                                                                    gcagcaTGAATGGGATTTCAGTGCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCAt                                                                                              >  2:395516/1‑300 (MQ=255)
                                                                                                                                                                                                            aaTGGGATTTCAGTGCGATAAATACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTTCGCACGACGGCTAAGCGCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCTGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCGACGATGCAGCACCGAAGGCTGAATATCCGGGCCCATCACGGGAgg                                                                                     >  1:647072/1‑301 (MQ=255)
                                                                                                                                                                                                                          gCGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTTACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAATCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCCATATCCATCACCGATCAGCGTCCTACGATGCAGCACCTAAGGCTGAATATCCGAGCCCCTCACGGTTGGTCTCCCCAgcgcc                                                                        >  1:576237/1‑299 (MQ=255)
                                                                                                                                                                                                                           cGATAAAGACCGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATAGCCGCGCCCATCACGGTAGGGATCCCCAGCGAAc                                                                       <  1:198165/300‑1 (MQ=255)
                                                                                                                                                                                                                                     cGATGCCATCAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGAt                                                             >  2:734972/1‑300 (MQ=255)
                                                                                                                                                                                                                                              cAATACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAGTCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCCTgg                                                   >  1:405832/1‑301 (MQ=255)
                                                                                                                                                                                                                                              aaaTACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATgg                                                   <  1:395512/300‑1 (MQ=255)
                                                                                                                                                                                                                                               cggtgactggagttcagacgtgtgctcttCCGATCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATggg                                                  <  2:770445/274‑1 (MQ=255)
                                                                                                                                                                                                                                               aaTACGGCTGCCCAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACAATGATCGCAGCATggg                                                  >  2:232754/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         cccAGTTTTTCTTCATGTTCCGGGCTTAAACCAGCATTGAGCATG‑CCTTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTAACGTCTCCGGATCGACCAGCAATTCACCGCGATAGCTATCAACGATCATCGTCCTACGATGCAGAACCGAAGGCAGAATATCTGCACCCAACACGGTAGGGATCCCCAGCGCACGTAACATGATCGCAGCATGCGCGGTGGcgt                                         >  1:788785/1‑298 (MQ=255)
                                                                                                                                                                                                                                                              tttttCTTCATGTTCCGGTCTTAAACCAGCATTGAGCATGAC‑GTTATACGCTCACCGCTTTGTAACTGCGCGGGTAAATTGACGTCATCGTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGTACGTACCATGATCGCAGCATGGGAGTTGGCTGCGcc                                     <  2:647077/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tcATGTTCCGGGCTTAAACCAGCATTGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACAATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCAc                            <  1:232750/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tcATGTTCCGGGCTTAAACCAGCATTGAGAATGAC‑CTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCCCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCAATCAACGAACAGCGTCCTACGATGCAGCACCGAAGGCTGAATAGCCCCGCCCATCACGCTAGGGATCCCCAGCGCACGTACCACGATCGCAGCATGGGAGTAGACTGCGCCATCTCGCa                             >  2:258257/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                gCTTAAACCGGCATTTAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCg                  <  2:788792/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             ttGAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCCTg    <  1:734966/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               gAGCATGAC‑TTTATACGCTCACCGCTTTTTAACTGCGCGGGTAAATTGACGTCATCTTCCGCCAGACGGCTAAGCTCAATCTCTTCACTAATTAGCCGCTGATATTCTTGCAGCAGTACCGGCTCCGGATCGACCAGCAATTCACCGCGATAGCCATCAACGATCAGCGTCCTACGATGCAGCACCGAAGGCTGAATATCCGCGCCCATCACGGTAGGGATCCCCAGCGCACGTACCATGATCGCAGCATGGGAGTTGGCTGCGCCATCTCGCACGACAACACCGACTAAGCGATCTTggg  >  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: ‑