Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,421,193 A→G 15.7% 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,1930AG15.7% 145.9 / 2.0 69S405G (AGT→GGT) ulaAL‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (39/16);  new base G (2/10);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 7.81e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCA  >  NC_000913/4420955‑4421424
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
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)
                                                                    cgtcccccgttctgtgttttCCTGCGTTTTTCCCGGTTGTTCTGCTCCACCGCCGCCAGGGGCGTTTTCGCTAACCTCTCCGGCGGCGGGCGTGGGGGTGTGAAGATTTGTCGGGTGGGGGGGATGATGAAAATCTTGGTTGGCGTCGGGGGGGTAAAACTCGCCGGTAGGAGTGCCTGGAGGGGCATGGGGGACTGGTCTTTTCTGCCACCGCCGATGATGCAAGGCTTCTTCTCCAt                                                                                                                                                                      <  2:318308/221‑1 (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)
                                                                                                                                                                 tcggcggctgagattttGTCTGGTGGGGGGGGTGTTGGAATTGTTGGGGTGGGTGTGGGGGGGAAACCGCCCGGGTATGTGTGCCTGGTGGGGCGGGGCGGCCTGGTGTTTTGGGCCCCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGaa                                                                           <  2:1009716/224‑1 (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)
                                                                                                                                                                                        cgcacaaaaagatgaaactttttGGTTGTGTGGGGGGGGGGAACCTCCGGGGTATGAGGGGTTGGTGGGGCAGGGCGGTCTGTGCTTTTCTGGCCCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGc                                                    <  2:763545/218‑1 (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)
                                                                                                                                                                                                                           gggggtgaaccccggcGGTATGCGTGCGTGGAGGGGCGTGGCCGACTGGTCTATTCTGGCTCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTAtt                  <  2:540393/223‑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/4420955‑4421424

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: