Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,421,192 A→G 16.4% S405G (AGT→GGT)  ulaA → L‑ascorbate‑specific enzyme IIC permease component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,421,1920AG16.4% 142.3 / 2.5 65S405G (AGT→GGT) ulaAL‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (39/14);  new base G (2/10);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 4.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCA  >  NC_000913/4420954‑4421423
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
catggttgggggCTGTATGTGGGTGACCATCGGGTGACTGATTGCGGTTGGCTGCCTGGTCGCCTGCGTCTCCTCGATCCTGATTATTCCTGTCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                           <  2:616711/231‑1 (MQ=255)
cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                           <  1:790346/239‑1 (MQ=255)
cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                           <  1:516092/239‑1 (MQ=255)
   ggTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGGTGCGTCTGGGCGGTGAAACTCACCGGGATGAggg                                                                                                                                                                                                                                        >  1:1099749/1‑237 (MQ=255)
        ggggCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTgg                                                                                                                                                                                                                                   <  2:809147/239‑1 (MQ=255)
          ggCTTTATGTGGGGCACCATCGGGCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTGCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGa                                                                                                                                                                                                                                  <  1:746507/238‑1 (MQ=255)
            cTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATgg                                                                                                                                                                                                                               >  1:908639/1‑239 (MQ=255)
            cTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGGGCCTGGGTgg                                                                                                                                                                                                                               >  1:549657/1‑239 (MQ=255)
               tATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGGAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGGAACTCACCGGTATGAGTGCCTGGGTGGGCa                                                                                                                                                                                                                            >  1:446796/1‑239 (MQ=255)
                     gggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCTCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCg                                                                                                                                                                                                                       <  1:947706/238‑1 (MQ=255)
                      ggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTTGGCGGCTGGCGTGCCGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGGAAACTTTGGGTGTGTCTTGGCGGTGAAACTCCCCGGTATTGGTGTCTGGATGGGCAGGGCgg                                                                                                                                                                                                                      >  2:1057385/1‑238 (MQ=255)
                           cATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGt                                                                                                                                                                                                                <  1:449752/239‑1 (MQ=255)
                            aTCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTc                                                                                                                                                                                                               <  1:473043/239‑1 (MQ=255)
                                gTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGAt                                                                                                                                                                                                            >  1:165089/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTgg                                                                                                                                                                                                      >  1:61975/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTgg                                                                                                                                                                                                      >  1:318308/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGGGATGGGGATGATCGGAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGCTTCTgg                                                                                                                                                                                                      >  1:655707/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCCTGCGGCGCTGAAGATTTGTCTGGGGATGGGGGTGATCGGAATCTTTTGTTGCGTCTGGGCGGGGAAACTCACCGGTATGAGTGCCTGGGTGGGCATGGCGGCCTGGTCGATTCTgg                                                                                                                                                                                                      >  1:424453/1‑238 (MQ=255)
                                        aTTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATTAGTg                                                                                                                                                                                                                                        >  2:179126/1‑202 (MQ=255)
                                        aTTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTg                                                                                                                                                                                                                                        <  1:179126/202‑1 (MQ=255)
                                                  gCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCAACGGTATGAGTGCCTGGATGGGCATGGGGGACTGGTCGATTCTGGAACCGCCggtggt                                                                                                                                                                                        >  2:1068823/1‑236 (MQ=255)
                                                         ggTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCCCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTTTCTGGTGATGGGGGTGATCGGAATCTTTGGTTGCGTCTGGGCGGTGAAAATACCCGGGATGAGTGCCTGGCTGGGCATGGCGGGCTGGTCGATTCTGGCCCCGGCGCTGATGGACGGct                                                                                                                                                                                >  2:378756/1‑239 (MQ=255)
                                                         ggTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGGATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGct                                                                                                                                                                                >  1:1009716/1‑239 (MQ=255)
                                                             gCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGGCTGGTCGATTCTGGCACCGCCGATGATGCAAGGcttct                                                                                                                                                                             >  1:1149221/1‑238 (MQ=255)
                                                                                  attattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGAATTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTAt                                                                                                                                                       >  2:742832/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGCTGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGGAACTCACCGGTATGAGTGCCTGGATGGGCCTGGCGGGCTGGGCGATTCTGGCACCGCCGGTTGTGCAAGGCTTCTTTTCCATCGGGATTGCCTTTTTGGCCg                                                                                                                                                  >  2:362345/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGGAATTTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTAGGGGTGCCTGGATGGGCATGGCGGGCTGGTCGGTTCCGGCACCGCCGCTTATGCAAGGCTTTGTCTCCATCGGGATTGCCTTTTTGGCCg                                                                                                                                                  >  2:1191813/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGCTGATGCAAGGCTTCTTCTCCATCGGGATTGCCTTTATGGCCg                                                                                                                                                  >  1:1191251/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTACTGCCg                                                                                                                                                  >  1:452263/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATAGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGTCTGGTCGCTTCTGGCCCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATTGCCTTTATGGcc                                                                                                                                                   >  1:254567/1‑238 (MQ=255)
                                                                                                 aTCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGCTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATAATTGt                                                                                                                                        >  1:378575/1‑239 (MQ=255)
                                                                                                 aTCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTTTGGCCGTCATCATTGt                                                                                                                                        >  1:746308/1‑239 (MQ=255)
                                                                                                  tCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGCATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCCCCTTTATGGCCGGCATCATTTTa                                                                                                                                       >  1:69244/1‑239 (MQ=255)
                                                                                                       ttGTTCTTCTCTAACGCCACCATCGGCGTTTTCGCTAACCACTTCGGCGGCTGGCGTGCGGGGCTGGAGATTTGACTGGTGGTGGGGATCAGCGGAATGTTTGGTTGCGTCTGGGCGGTGGATCTCACCGGTATGAGTGGCCGGATGGGCATGGCGGCCTGGTGGTTTCTGGCACCGCCGATGACGTGAGGCAGCTTCTTCATCAGCATCGCCTTTATGCGCGTCATCATAGTaaatc                                                                                                                                   >  2:227532/2‑235 (MQ=255)
                                                                                                       aTGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGc                                                                                                                                  >  1:793225/1‑239 (MQ=255)
                                                                                                        ttttcttcGCTAACGCCCCCACAGGCGTGTTCCCTAACCACTTCGGCGGCTGGCTTGCGGCGCAGAGGTTTTGTCTGGTGTTGGGGTGGATCGAAAGCTTTGGTTGCGTCGGGGCGGTGAACCTCCCGGGTAGGAGTGCCTGGATGGGCTGGGCGGACTGGTCGATTCTGGCGCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                                 <  2:746308/237‑1 (MQ=255)
                                                                                                         gTTCTCCTCTAACCCCCCCATGGGGGGTTCGGCTCACCGTTGGGGCGGCGGGGGGGGGGCGCGAATGTTTGGGCTGGTGATGGTGATTAGCGAACGCTGTGGGGGGGTGGGGGCGTGTAAACACTCCGGTGGGGGTGCCGGGGGGGGTGTGGCGGGGTGTCCTATTCTGTCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCAc                                                                                                                                <  2:1099749/239‑1 (MQ=255)
                                                                                                                    aCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGTTGATCGAAATCTTTGGTTGCGTCTGGGCGGGGAAACTCACCGGGATGAGGGCCTGGATGGGCATGGCGGACTGGTCGGTTCTGGGCACGCCGGTGATGAAAGGCTTCTTCTCCATCGGGATCGCCTTTATGGCCGTCATCATTGGAATTTCACTGGCCTATATg                                                                                                                     >  2:845292/1‑239 (MQ=255)
                                                                                                                         accaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGttct                                                                                                                 >  1:1029899/1‑238 (MQ=255)
                                                                                                                            aTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGGGATGGGGCTGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGGCGATTCTGGCACCGCCGATGGTGCAAGGCTTCTTTTCCATCGGGATCGCCTTTATGGCCGTCATCATTGTAATTGCCCTGGCTtatc                                                                                                                       >  1:626850/1‑228 (MQ=255)
                                                                                                                              cGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTTGCTg                                                                                                           >  1:431068/1‑239 (MQ=255)
                                                                                                                                             ccctctggggcgggTGGGGTGCGGCGGTGAAGTTTTTTGTGGGGATGGGGTTGGTCAAATTCTTTGGTGGCGCGGGGGCGGTAAATCTCACGTGTATGGGTGCCTGGAGGGGCTGGGCGGACTGGGCTTTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGCTTATATGTTCTTCGCTGGCCGCGCGCTgcgcg                                                                                            <  2:793225/231‑1 (MQ=255)
                                                                                                                                                aCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGTCTTTATGGCCGTCATCATTGTAATTGCAATGGCCTATTTGTTGTTCCCTGGCCGCGGGCTGCGCGCa                                                                                          >  2:1084938/1‑239 (MQ=255)
                                                                                                                                                  ttCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGAGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCaga                                                                                        >  2:135950/1‑239 (MQ=255)
                                                                                                                                                         gcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGAt                                                                                 >  1:934348/1‑239 (MQ=255)
                                                                                                                                                           gggCGTGGGGCGGTGAAGTTTTGGCGGGTGGTGGGTATGTGCAAATTCTTTGGTTGCGTGTGGGGGGGGAAACCCCCCGGTATGGGTGCCTGGATGGGCTTGGCGGACTGGTCGTTTCTGCCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATg                                                                                <  2:921117/237‑1 (MQ=255)
                                                                                                                                                            ggCGTGCGGCGCTGAAGTTTTGTCTGGTGCTGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGTCGGGCCGATTCTGGCCACGCCGATGATGCAAGGCTTCTTCTCCATCGCGATAGCCTTTCTGGCCGTCATCATGGTAATTGCACTGGCTTATATGGTCTTCGCTGCCCGCGCCCTGCACGCTGAAGAAGCTGcg                                                                              >  2:589577/1‑238 (MQ=255)
                                                                                                                                                                       cTGAAGATTTTTCGGGTGATGGGGATGACAGAAATCTTTGGTGGCGTCTGGGCGGTCAAACTCACCGGTATGAGTGCCGAGAGGGGAATGGCGGACTGGTCTATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaacaac                                                                    <  1:1057385/238‑1 (MQ=255)
                                                                                                                                                                              tttGTCGGGTGATGGGTATGATCGAATTTTTTGGTGGCGTCTGGGCGGTGAACCTCACCGGTATGGGTGCCGGGAGGGGCATGGCGGCCTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGa                                                             <  1:292938/238‑1 (MQ=255)
                                                                                                                                                                                  gcgggTGATGGGGATGATAGAATTCTTGGGTGGCGTCTGGGGGGTAAAACGCCCCGGTATGAGTGCCGGGATGGGCATGGCGGACTGGTCGATTCTGGCGCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                        <  1:1068823/236‑1 (MQ=255)
                                                                                                                                                                                  ctgggTGATGGGGTTTAGGAAAATTTTTGGTTGCGTTGGGGCGGTGAACCTCCCCGGTATGGGTGCTTGGATGGGCTTGGCGGACTGGTCTTTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                        <  1:1191813/236‑1 (MQ=255)
                                                                                                                                                                                   ttGGTGATGGGGATAAGCGAAATCTGTGGGGGGGTGGGGGCGGTAAAACTCCCGGGTAGGGGTGCCGGGTGGGGCATGGCGGCCTGTGCTTTTCTGGCTCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                        <  2:1149221/237‑1 (MQ=255)
                                                                                                                                                                                    tGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAg                                                         >  2:917006/1‑236 (MQ=255)
                                                                                                                                                                                         aTGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTta                                                 >  2:822326/1‑239 (MQ=255)
                                                                                                                                                                                             ggATGATAGAAGTGTTTGGTTGCGTCTGGGCGGTGAAGCTCACCGGTAGGAGGGCGGGGAGGGGGAGGGGGGACTGGTCGATTGTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTtaata                                              <  2:828418/238‑1 (MQ=255)
                                                                                                                                                                                                 gatCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAg                                         >  1:460879/1‑239 (MQ=255)
                                                                                                                                                                                                   cctaaaTGTTGGGTTGCGTGGGGGGGGTGAAATCCCCGGGTATGGTGGTGTGGATGGTCGTGTGGGCCTGGTGTATCGTGACACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                                       <  2:460879/236‑1 (MQ=255)
                                                                                                                                                                                                      gaaTCTTTGGGTGGGGTCGGGGCGGTCAACCTCACCGGTTGGGGTGCCGGGAGGGGAATGGCGTACGGGGGGTTGTGGGACCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                                       <  2:378575/236‑1 (MQ=255)
                                                                                                                                                                                                          cTTTGGTTGCGTCGGGGGGGTGAAGCTCGCCGTTTGGAGGGTGGGGAGGGGCGTGGCGTCCGGTGCGTTTCTGTCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTa                                <  2:1202604/239‑1 (MQ=255)
                                                                                                                                                                                                           tttgggTTGCGCCGGGGAGGTGAAAGTCACCGGTAGGAGGGCCGGGAGGGGCATGGCGGGCTGGTCTATTCGGGCACCGCCGATGATGCAAGGCTTCGTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTa                                <  1:378756/235‑1 (MQ=255)
                                                                                                                                                                                                             tGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATg                              >  2:830634/1‑238 (MQ=255)
                                                                                                                                                                                                                gTTGCGTCTGGCCGGGAAAACTCCCGGGTATGAGTGCTGGGATGGGCATGGCGGACGGGTCAATGCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGAc                            <  1:44204/238‑1 (MQ=255)
                                                                                                                                                                                                                    gtctggggCGGGGAAACCCGCCGTGATGGTTGGCTGGAGGGGTGGGGCGTGCTGCTCATTCGTGGGACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACcgta                        <  2:337543/233‑1 (MQ=255)
                                                                                                                                                                                                                    agTCTGGGCGGCAAAACTCAGCGTATTGGGTGCCGGGATGGAGTTGGGGGCTTGTTCGATTTTGTCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACcgta                        <  1:195935/237‑1 (MQ=255)
                                                                                                                                                                                                                               tGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTAtt                  >  1:520209/1‑233 (MQ=255)
                                                                                                                                                                                                                               tGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTg               >  1:540393/1‑236 (MQ=255)
                                                                                                                                                                                                                                         gggggggAGGGCTTGAATGGGCGGGGGGGACGGGCCGATGGTGGAACCGCCGATGATGCAAGGTTTCATCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGGTATGATTATGACCGTACGTATTCTGGCTGTGTGTGGCa  <  1:227532/233‑1 (MQ=255)
                                                                                                                                                                                                                                              tGAGTGCCTGGATGGGCATGGCTGACTGGTCGATTCTGGCACCGCCGATGGTGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAACACAACTGGCAGAACAGTCTGCTtaa                                                >  1:649853/1‑188 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
CGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCA  >  NC_000913/4420954‑4421423

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