Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,079,070 A→C 13.3% G87G (GGT→GGG cmtB ← putative mannitol‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,079,0700AC13.3% 135.2 / 3.3 70G87G (GGT→GGGcmtBputative mannitol‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base A (35/21);  new base C (3/9);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 2.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAAC  >  NC_000913/3078843‑3079291
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
ggcggcTTGTTTTCAATAATTTTAGTCTGTAGGGGGTGAGATAATTAGCCGCGGCTGTTAATGACCGCTAAGTGTTTTTCTGGTGATGTTGTTAAGAGTTGTGCGAGTATTTCTTCTTCGCACAGTAATTCCCTTAACGCCTTAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAaca                                                                                                                                                                                                                                <  2:193553/223‑1 (MQ=255)
gCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAacac                                                                                                                                                                                                                               >  1:671594/1‑229 (MQ=255)
                      tatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCACTGTGCGAATCGGCATAGGCAGCAGAGAGTCCGTTGAGTAATTTAATTGGTTCGTCATTCCCCGGCAAATAACCACCTTGTTCAAGTAATGTCaaaaa                                                                                                                                                                                                        >  2:634979/1‑227 (MQ=255)
                      tatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGCTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCCATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATACACCCCTTGTTCAAGTAATGTCAAAgc                                                                                                                                                                                                        >  1:334934/1‑228 (MQ=255)
                        tCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCCGAGAGTCCGATGAGTAATTTAATTGGTTGTTCATTCCCCGGAAAATAAACACCTTGTTCAAGAAATGTCaaaaaaa                                                                                                                                                                                                      >  2:656433/1‑225 (MQ=255)
                        tCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGTTTCGTCATTCCCCGGAAAATAAAAACCTTGTTCAAGTAATGTCAAAcaca                                                                                                                                                                                                      >  1:634835/1‑229 (MQ=255)
                             tAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCCGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAAGGTCAAAGCCATccc                                                                                                                                                                                                  >  1:461324/1‑228 (MQ=255)
                             tAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATccc                                                                                                                                                                                                  >  1:724822/1‑228 (MQ=255)
                             tAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATccc                                                                                                                                                                                                  >  1:524249/1‑228 (MQ=255)
                             tAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGAAATGTCAAAGACATccc                                                                                                                                                                                                  >  1:124495/1‑228 (MQ=255)
                                  gagagaGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGACAATAAACACCTTGTTCACGTAATGTCAAAGACATCCCGGttt                                                                                                                                                                                             >  1:641411/1‑228 (MQ=255)
                                     agagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTa                                                                                                                                                                                                                  >  2:69520/1‑204 (MQ=255)
                                     agagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAg                                                                                                                                                                                         <  1:571911/229‑1 (MQ=255)
                                     agagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTGCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTa                                                                                                                                                                                                                  <  1:69520/204‑1 (MQ=255)
                                       agTAATTATCCGCGGCTGATAATTTCCGCTAATTGTTTTTCTGCTGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGGAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGCTGGATAATTTAATTGGGTCGTCATTCCCCGGCAAGTAAACCCCTTGTTCCAGAAATGGCCAAGACAACCCGGGTTTataag                                                                                                                                                                                       >  2:660830/1‑225 (MQ=255)
                                       agTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGAAATCCCGGTTTTAAgg                                                                                                                                                                                        >  1:555894/1‑227 (MQ=255)
                                       agTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCCGGAAATAAACACCTTGTTCAATAAATGTCACAGACAACCCGGTTTTaaagc                                                                                                                                                                                       >  2:391256/1‑226 (MQ=255)
                                         ggaaTTTTCCGTGGGTGATAATGTCTGCTAATTGGTTTTCTGATGAGTGTGGTAAGAGTGTATCGAGGATTTCTGCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTGGTTCAAGTAATGTCAAAGACATCCCGGTGTTAAgcg                                                                                                                                                                                       <  2:355718/226‑1 (MQ=255)
                                          tAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGCGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGCGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAACACAAACCCCGTGTTCAAGTAATGTCAAAGAAAACCCGGTTTTAAGCAcc                                                                                                                                                                                     >  2:291666/1‑229 (MQ=255)
                                                  cgcgGCTGTTAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATAAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGCAATCCCGGTTTTAACAGCCCCACAATc                                                                                                                                                                             >  2:399747/1‑229 (MQ=255)
                                                  cgcgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGTCATCGGCAGCAGAGAGTCCGATAAGTAATTTAATTGGTTCGCATTTCCCCGGAAAATAACACCTTTGTTAAAGTAATGTCAAGAACATCCCGGTTTAAAGCACCCAACATTc                                                                                                                                                                             >  2:244774/1‑229 (MQ=255)
                                                  cgcgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCCGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGACAATAAACACCTTGTTCAAGAAATGTCAAAAAAATCCCGGTTTTAAACCCCCCAAATTc                                                                                                                                                                             >  2:725343/1‑229 (MQ=255)
                                                  cgcgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGGGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGGCCCACATTc                                                                                                                                                                             >  2:359980/1‑229 (MQ=255)
                                                      gCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGt                                                                                                                                                                         <  1:649103/229‑1 (MQ=255)
                                                                  cgttatttGTTTTTCTGAGTATGCTGTTCAGAGTTTTTCGAGTTTTTCTTCTGCGGACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGc                                                                                                                                                             <  2:91480/223‑1 (MQ=255)
                                                                    cTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATCCCGGACGGCCATAAGGc                                                                                                                                                             >  1:91480/1‑227 (MQ=255)
                                                                    cTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGCATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTTCCGCATGAgg                                                                                                                                                              >  1:713750/1‑226 (MQ=255)
                                                                                  gatgatTCTTTTAAGAGGTGGTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTg                                                                                                                                             <  1:391256/229‑1 (MQ=255)
                                                                                   atgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGACAATAAACACCTTGTTAAAGTAATGTCCAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGCAGCGCATGAGGACTTGCCAccaccgtgac                                                                                                                                            >  2:523561/1‑220 (MQ=255)
                                                                                   atgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAATAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGc                                                                                                                                            <  2:124495/229‑1 (MQ=255)
                                                                                   aggttgCTGTGAAGAGTTGGTCGGGCATTTCTTCTTCGCACAGTAATTCACTGAACGCCTGAATAGCGCCAATGGGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACCCCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGc                                                                                                                                            <  2:713749/225‑1 (MQ=255)
                                                                                     gatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCga                                                                                                                                          <  1:634980/229‑1 (MQ=255)
                                                                                        tcttttAAGAGTTGTGCGATTTTTTCTTTTCCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGaa                                                                                                                                       <  2:634834/225‑1 (MQ=255)
                                                                                                gagTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATAtaataa                                                                                                                               <  1:659322/229‑1 (MQ=255)
                                                                                                     ttttgtgTATTTCGTCTTCCCTCAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTAGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACAGCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGcc                                                                                                                          <  2:562476/223‑1 (MQ=255)
                                                                                                     gttCGAGTATTTTTTCTTCGCACAGTAATTCCCTTAACGCCGGAATAGCGCCAATGTGAGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGcc                                                                                                                          <  2:461324/229‑1 (MQ=255)
                                                                                                      ttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGGGAATATAATAAGGGCCa                                                                                                                         >  1:25284/1‑229 (MQ=255)
                                                                                                      ttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATCAGGCATTGCCACGCCTGGTGCGCGAATATCATACGGGCCa                                                                                                                         >  1:189886/1‑229 (MQ=255)
                                                                                                                    cttcGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGACAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCCTTGTTAAtggtgg                                                                                                            >  1:398001/1‑228 (MQ=255)
                                                                                                                      tcGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGATATTATTAAGGGCCATTGTTAATGGTGGaa                                                                                                          >  1:430794/1‑228 (MQ=255)
                                                                                                                                aTTCAATTAACGCCTGAATAGCGACAATGTGCCAATCGGACACGGCAGCAGAGAGTACCATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTGAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCACACATTCCGGTCGCGCATGAAGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGc                                                                                               <  1:607212/229‑1 (MQ=255)
                                                                                                                                  tcccTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCt                                                                                             <  1:725344/226‑1 (MQ=255)
                                                                                                                                      ttAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGGGGAATCTTTAATTGCCTGAAt                                                                                         >  1:506194/1‑229 (MQ=255)
                                                                                                                                      ttAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGGGCGGGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGaa                                                                                          >  1:455378/1‑228 (MQ=255)
                                                                                                                                      ttAACGCCTGAAGAGTGGCAATGTGCGAATCGGCATCGGCAGCAGAGTGTCCGATGAGTAAATTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACCTTCAGGTCGCGCCTGACGCATTGCCACGCCTGGTGCGGGACTTAAAAAAGGGCAATTGTTAATGTGGAATGCTTAAATTCCCGaa                                                                                          >  2:269102/1‑227 (MQ=255)
                                                                                                                                        aaCGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTACTTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTCCCACGCCTGGGGCGAGAATATATTAAGGGCCATTGTTAATGGGGGAATCTTTAATTGCCTGAATGt                                                                                       >  2:490836/1‑229 (MQ=255)
                                                                                                                                        aaCGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCCCGCCTGGGGCGACAATATAATAAAGGCCATTGTTATTGGGGGAACCTTTAATTGCTTGAAtta                                                                                       >  2:84662/1‑227 (MQ=255)
                                                                                                                                        aaCGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGaa                                                                                          >  2:443606/1‑226 (MQ=255)
                                                                                                                                        aaCGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTAGAGCATGAGGCATTGCCAAGCCTGGTGCGCGAATATAATAAGGGCCATTGTTATGGAGGGACTATTTATTTGCCTGaagta                                                                                       >  1:756590/1‑226 (MQ=255)
                                                                                                                                        aaCGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGGCGCGCATGAGGCATTGCCACGCCTGGTGCGCGAATCTAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGt                                                                                       >  1:550758/1‑229 (MQ=255)
                                                                                                                                                     cgcCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAAttatagcgggt                                                                          >  1:480805/1‑220 (MQ=255)
                                                                                                                                                       cccatggtggAATCGGCATCGGCAGCAAAGAGGCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACCCCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGatat                                                                        <  2:376353/220‑1 (MQ=255)
                                                                                                                                                          aTGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGt                                                                     >  1:244203/1‑229 (MQ=255)
                                                                                                                                                          aTGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGGGGAATCTTTAATTGCCGGAATGTAATTTTCGCTGATATAGt                                                                     >  2:637682/1‑229 (MQ=255)
                                                                                                                                                                  atcggcatcggcAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTATTGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATcc                                                             <  1:656434/229‑1 (MQ=255)
                                                                                                                                                                  atcggcaatcgcAGAAAAAAGTCCGAGTAGTAATTTAATTGGCTCGTCATTCCCCGGAAAATAAACCCCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATcc                                                             <  2:524249/229‑1 (MQ=255)
                                                                                                                                                                     ggcctcggagGCAGAGGGTCCCATTGTTAAATTAATTGGTTCGTCATTCCCCGGAAAATAAACCCCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGc                                                          <  2:642403/229‑1 (MQ=255)
                                                                                                                                                                          cggcaGAAGAGGTTCCATTGAGGAATTTAATTGGTTCGGCAATCCCCGTAAAATAACCACCTGGTTCAAGTAATGTAAAAGACTTGCCGGTTTCAAGCGCCCCACATTCCGGTCGCGCATGAGTCATTGCCACGCCGGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTTGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATg                                                      <  2:550758/228‑1 (MQ=255)
                                                                                                                                                                                            tgatttatttaagtGCTTCGCACTACCACGGAAAATAAACCCCTTGTACAAGGAAGGTCAAAGGCATCCCGGTTTTAAGCGCCCCACATTCGTGCGGCCCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCGGAATGTAATTCTCGATTATATAGTTTTTATCCAGAAATGATACCATCGAGAAATCGAt                                   <  2:205917/223‑1 (MQ=255)
                                                                                                                                                                                              tgaaatttaattGGTTCGTCTTTCCCCGGAAAATAAACACCTTGGTCAAGTAATGCCAAAGTCTTCCCGGGTTTAAGCGCCCCATATTCCGGTCGCGCATGAGGCATTGCCCCGCCTGGTGCGAGAATATAATAAGGTCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGAt                                   <  2:189886/224‑1 (MQ=255)
                                                                                                                                                                                                taatttaattGGTTCTTCATTCCCGGAAAAATGACCTCCTTTTTCAAATAATGTCAAAGACATCCCGTTTTTAAGCGCCCCACATGCCGGTCGCTCATGATGCATTGCCACGCCTTGTGCGGGAATATAATAAGGGCCATTGTTAATGTTGGAATCTTTAATTGCCTGAATGTAATTCTCGTTGCTATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCt                               <  2:724821/229‑1 (MQ=255)
                                                                                                                                                                                                      aattGGTTCGTTATTCCCCGGAAAATAAACCCCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGc                         <  1:637683/229‑1 (MQ=255)
                                                                                                                                                                                                        ttGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATc                                              >  1:423695/1‑206 (MQ=255)
                                                                                                                                                                                                        ttGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCAGTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATc                                              <  2:423695/206‑1 (MQ=255)
                                                                                                                                                                                                          ggTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCAAGCCTGGTGCGCGAATATAATAAGGGCCATTGTTAATGGTGGCATCTTTAATTGCCTGAATGTAATTTTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAAt                     >  2:541597/1‑229 (MQ=255)
                                                                                                                                                                                                                  aTTCCCCGGAAAATAACCCCCTTGTTCAAGTAATTGCAAAGACCTCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCt             <  2:244203/229‑1 (MQ=255)
                                                                                                                                                                                                                  aTTCCCCGGAAAATAAACTCCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTGAAGCGCCCCCCATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCt             <  1:523561/229‑1 (MQ=255)
                                                                                                                                                                                                                    ttcccGGAAAATTAACCCCTTTTGCAAGTAAGTTCAAAGAACTCCCGGTTTTAAGCTCCCCACATCTCGGGTGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTGTAATTGCCGGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGa           <  2:506194/227‑1 (MQ=255)
                                                                                                                                                                                                                    tCCCCTGAAAATAACCACCTTTTTCAATGAATGACAAAGACATCCCGTTTTGAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCAGTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGa           <  2:555893/229‑1 (MQ=255)
                                                                                                                                                                                                                    tCCCCGGAAAATAACCACCTTGTTCAAGTAATGTCAAAGCCATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGTTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGa           <  1:660831/229‑1 (MQ=255)
                                                                                                                                                                                                                    tCCCCGGAAAAAAAACCCCTTGGTCAAGGAATGTCAAAAGCATCCCGGTGTTAAGCGCCCCAAATTCCGGTCGCGCATGAGGCATTGCCACGCCTGTTGGGAGAATATAATAAGGACCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGa           <  2:398001/229‑1 (MQ=255)
                                                                                                                                                                                                                             aaTAACTACCTTGTTCAAGTACTGTCAAACGCATCCCGGTTTTAAGCGCCCCACATTCCAGTCGCGCATGAGGTATTGCCACGTCTGGTGCGAGAATATAATAAGTGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATTATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAAc  <  1:359980/229‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
GCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTTTAAGCGCCCCACATTCCGGTCGCGCATGAGGCATTGCCACGCCTGGTGCGAGAATATAATAAGGGCCATTGTTAATGGTGGAATCTTTAATTGCCTGAATGTAATTCTCGCTGATATAGTTTTTATCCAGCAATGATACCATCGAGAAATCGATAGCTTCCTGCCAATCTTTTGCTGAATGTATAAC  >  NC_000913/3078843‑3079291

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