Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,141,053 G→A S390F (TCC→TTC)  ptsA ← putative PTS multiphosphoryl transfer protein PtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,141,0530GA93.0% 121.7 / 1.9 44S390F (TCC→TTC) ptsAputative PTS multiphosphoryl transfer protein PtsA
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (18/22);  total (20/24)
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.46e-02

CGATAACCGAGGAACGGGTTTGCCTCTGCGGGAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAAC  >  NC_000913/4140813‑4141298
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
cgATAACCGAGGAACGGGTTTGCCTCTGCGGGAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGAt                                                                                                                                                                                                                                             <  1:211674/251‑1 (MQ=255)
      ccGAGGAACGGTTTTGCCTCTGCGGGAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGc                                                                                                                                                                                                                                       <  1:166035/251‑1 (MQ=255)
       cGAGGAACGGGTTTGCCTCTGCGGGAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCg                                                                                                                                                                                                                                      >  1:73489/1‑251 (MQ=255)
                      ccTCTGCGGGAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAGTAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAACTGTTGTACAACTCGCTTTAGCCCGGTGCGCTGGTGCGATCCATCTAGAGCATTTCAGTGCGGAACACACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCACTTTCAATGCgg                                                                                                                                                                                                                       >  1:290600/1‑251 (MQ=255)
                             gggaaagatcagatcctCAACGGGTTTGTCGCCGCCACTGCCCCTAGAGCGCACACTAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAACTGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGCGCATTTCAGGGCGGAACAAACCAACGCCTTCCGCACCATTGCCGACTGAGGCCTGGGCTTCCACGCAGTGAACGATGTTAGCGGCACTTTCAGTTCGGATACCGt                                                                                                                                                                                                                >  2:344869/17‑251 (MQ=255)
                               gaaggTTCAGAGAATAAACTGGTTGGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCCGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCCTTCGCCCGGAGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTTCGCTTCCACGAACTGAGCGATGTTAGCGGCAATTTCaacgcttataccttca                                                                                                                                                                                                              <  2:226822/247‑15 (MQ=255)
                               gAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAATTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCa                                                                                                                                                                                                              <  1:119800/251‑1 (MQ=255)
                               gAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCATTGCGGATACCGTc                                                                                                                                                                                                               >  2:344410/1‑250 (MQ=255)
                                   gTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCgg                                                                                                                                                                                                                       >  1:79425/1‑238 (MQ=255)
                                   gTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCgg                                                                                                                                                                                                                       <  2:79425/238‑1 (MQ=255)
                                      cataTAATCAACTGGTTGGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTTCAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCTCCATTGCCGTATGCGGCCTGCGCTTTCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTAc                                                                                                                                                                                                       <  1:150294/248‑1 (MQ=255)
                                      cAGATAATCAACGGTGTTGTCGCCGCCAATTTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTa                                                                                                                                                                                                        <  1:13789/250‑1 (MQ=255)
                                          tAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGACTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGAAATTTCAATGCGGATACCGTCAGCGGTACGGGc                                                                                                                                                                                                   >  1:289699/1‑251 (MQ=255)
                                                       tGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGCACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCACTTTCCATGCGGTTACCGTCAGCGGTACGGGCTTGTTGGGTCAGc                                                                                                                                                                                      >  2:343008/1‑251 (MQ=255)
                                                       tGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGGTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACTCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGAACGGGCtttgtgatacacc                                                                                                                                                                                      >  2:181610/1‑240 (MQ=255)
                                                             cgccAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGaca                                                                                                                                                                                <  1:344410/251‑1 (MQ=255)
                                                                       cATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCt                                                                                                                                                                      >  1:58646/1‑251 (MQ=255)
                                                                           gTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACg                                                                                                                                                                  >  1:304551/1‑251 (MQ=255)
                                                                                      aatGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCt                                                                                                                                                       >  1:28243/1‑251 (MQ=255)
                                                                                                         gATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCt                                                                                                                                     <  2:84026/250‑1 (MQ=255)
                                                                                                                 ggtctgaccaaaaatgtcgttgcaCTCGCTTTCGCCCGCGTCGTTTGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCATCGACATACTGACCACTGCCGAATGCGTCCTGCGCTTTCACGAAGTGAGCGATGTTAGCGTCAATTTCAATGCGGATTCCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGGTGTTGCGCCCGTTGGGCGTCCTGTACGCGGGCTTCGTTTTGATAATAAc                                                                                                                            <  2:152403/228‑1 (MQ=255)
                                                                                                                  gcCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACg                                                                                                                           >  1:159121/1‑251 (MQ=255)
                                                                                                                          aaaaaTGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACgg                                                                                                                   <  1:106421/251‑1 (MQ=255)
                                                                                                                                                   cccGGTGCGCTGGTGCGATCCATATAGAGAATTTCAGTGCGGAACAAACCAAAGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGGTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGAGAATAACGGGCTAAGGCGTCCCCCGGCTAAAcacaaatcgc                                                                                          >  2:257683/1‑242 (MQ=255)
                                                                                                                                                       ggggcgcTGGTGCGATCCATATAGAGCATTGCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGcccc                                                                                       <  2:304551/248‑1 (MQ=255)
                                                                                                                                                        gTGCGCGGGTGCGAGCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCg                                                                                      <  1:157085/251‑1 (MQ=255)
                                                                                                                                                          gcgcTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGc                                                                                    <  1:241728/251‑1 (MQ=255)
                                                                                                                                                           cgctgtggCGATCCCTATATCGCAATTCAGTTCGGAACAAGCCAACGCCTTCCGCACCATTGCCGAAGGCGGCCTGCGCTTCCACGAAGGGCGCGATGTTAGCGGCAATTGCAATGCGGATCCCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCTGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCGCAACCACAATCGCCCCGGct                                                                                   <  2:242388/244‑2 (MQ=255)
                                                                                                                                                                      tCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGt                                                                            >  1:11564/1‑247 (MQ=255)
                                                                                                                                                                            taGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGacac                                                                                                                                                                               <  1:65269/142‑1 (MQ=255)
                                                                                                                                                                            taGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGacac                                                                                                                                                                               >  2:65269/1‑142 (MQ=255)
                                                                                                                                                                                          gCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGcc                                                    >  1:135744/1‑251 (MQ=255)
                                                                                                                                                                                           cGGAACAAACCAACGCCTTCCGCTCCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTGAGCGGCAATTGCAATGCGGATACCGTCAGCGGTACGGGCGTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTTATAATAACGAGCTACGGCTTCCCCTGGCTCATCCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCa                                                   <  2:42379/251‑1 (MQ=255)
                                                                                                                                                                                                ttgctgtttggcTTCCGCACTATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACACAGGGCGTCCTGTACGCGGGcttctt                                                                                                                                        <  2:172773/150‑1 (MQ=255)
                                                                                                                                                                                                ttgctgtttggcTTCCGCACTATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACACAGGGCGTCCTGTACGCGGGcttctt                                                                                                                                        >  1:172773/12‑161 (MQ=255)
                                                                                                                                                                                                cAAACCAACGCCTTCCGCACCATTGCCTAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAg                                              <  2:238696/251‑1 (MQ=255)
                                                                                                                                                                                                 aaaCCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCt                                                                                                       <  1:103451/193‑1 (MQ=255)
                                                                                                                                                                                                 aaaCCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCt                                                                                                       >  2:103451/1‑193 (MQ=255)
                                                                                                                                                                                                                  aCCATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCca                            <  1:217523/251‑1 (MQ=255)
                                                                                                                                                                                                                    cATTGCCGAATGCGGCCTGCGCTTCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCCACCACAATCGCCCCGGCGTTGCCGTCGGTATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCCTAAATATacacc                          >  2:264980/1‑250 (MQ=255)
                                                                                                                                                                                                                           gAATGCGGCCTGCGCTTCCCCGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCGGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACaccaacca                   <  1:310411/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggCGGCCTGCGCTTCCACGAAGTTAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTTCGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                <  2:199293/250‑1 (MQ=255)
                                                                                                                                                                                                                                           tCCACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACTCGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTga    <  1:238655/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ccACGAAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCCGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTgaac  >  2:36424/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
CGATAACCGAGGAACGGGTTTGCCTCTGCGGGAATGTTCAGATAATCAACGGGTTTGTCGCCGCCAATGTCCATAGTGCGCACAATAATGCTGCGTCCGTTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAAC  >  NC_000913/4140813‑4141298

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