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

CGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCA  >  NC_000913/4420955‑4421424
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
catggttgggggCTGTATGTGGGTGACCATCGGGTGACTGATTGCGGTTGGCTGCCTGGTCGCCTGCGTCTCCTCGATCCTGATTATTCCTGTCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                           <  2:616711/231‑1 (MQ=255)
cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                           <  1:790346/239‑1 (MQ=255)
cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                           <  1:516092/239‑1 (MQ=255)
   ggTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGGTGCGTCTGGGCGGTGAAACTCACCGGGATGAggg                                                                                                                                                                                                                                        >  1:1099749/1‑237 (MQ=255)
        ggggCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTgg                                                                                                                                                                                                                                   <  2:809147/239‑1 (MQ=255)
          ggCTTTATGTGGGGCACCATCGGGCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTGCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGa                                                                                                                                                                                                                                  <  1:746507/238‑1 (MQ=255)
            cTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATgg                                                                                                                                                                                                                               >  1:908639/1‑239 (MQ=255)
            cTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGGGCCTGGGTgg                                                                                                                                                                                                                               >  1:549657/1‑239 (MQ=255)
               tATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGGAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGGAACTCACCGGTATGAGTGCCTGGGTGGGCa                                                                                                                                                                                                                            >  1:446796/1‑239 (MQ=255)
                     gggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCTCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCg                                                                                                                                                                                                                       <  1:947706/238‑1 (MQ=255)
                      ggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTTGGCGGCTGGCGTGCCGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGGAAACTTT‑GGGTGTGTCTTGGCGGTGAAACTCCCCGGTATTGGTGTCTGGATGGGCAGGGCgg                                                                                                                                                                                                                      >  2:1057385/1‑238 (MQ=255)
                           cATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGt                                                                                                                                                                                                                <  1:449752/239‑1 (MQ=255)
                            aTCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTc                                                                                                                                                                                                               <  1:473043/239‑1 (MQ=255)
                                gTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGAt                                                                                                                                                                                                            >  1:165089/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑Tgg                                                                                                                                                                                                      >  1:61975/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑Tgg                                                                                                                                                                                                      >  1:318308/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGGGATGGGGATGATCGGAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGCTTC‑Tgg                                                                                                                                                                                                      >  1:655707/1‑238 (MQ=255)
                                     cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCCTGCGGCGCTGAAGATTTGTCTGGGGATGGGGGTGATCGGAATCTTT‑TGTTGCGTCTGGGCGGGGAAACTCACCGGTATGAGTGCCTGGGTGGGCATGGCGGCCTGGTCGATTC‑Tgg                                                                                                                                                                                                      >  1:424453/1‑238 (MQ=255)
                                        aTTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATTAGTg                                                                                                                                                                                                                                        >  2:179126/1‑202 (MQ=255)
                                        aTTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTg                                                                                                                                                                                                                                        <  1:179126/202‑1 (MQ=255)
                                                  gCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCAACGGTATGAGTGCCTGGATGGGCATGGGGGACTGGTCGATTC‑TGGAA‑CCGCCggtggt                                                                                                                                                                                        >  2:1068823/1‑236 (MQ=255)
                                                         ggTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCCCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTTTCTGGTGATGGGGGTGATCGGAATCTTT‑GGTTGCGTCTGGGCGGTGAAAATACCCGGGATGAGTGCCTGGCTGGGCATGGCGGGCTGGTCGATTC‑TGGCC‑CCGGCGCTGATGGACGGct                                                                                                                                                                                >  2:378756/1‑239 (MQ=255)
                                                         ggTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGGATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGct                                                                                                                                                                                >  1:1009716/1‑239 (MQ=255)
                                                             gCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGGCTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGcttct                                                                                                                                                                             >  1:1149221/1‑238 (MQ=255)
                                                                    cgtcccccgttctgtgttttCCTGCGTTTTTCCCGGTTGTTCTGCTCCACCGCCGCCAGGGGCGTTTTC‑GCTAACCTCTCCGGCGGCGGGCGTGGGGGTGTGAAGATTTGTCGGGTGGGGGGGATGATGAAAATCTTG‑GTTGGCGTCGGGGGGGTAAAACTCGCCGGTAGGAGTGCCTGGAGGGGCATGGGGGACTGGTCTTTTC‑TGCCA‑CCGCCGATGATGCAAGGCTTCTTCTCCAt                                                                                                                                                                      <  2:318308/221‑1 (MQ=255)
                                                                                  attattCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGAATTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTAt                                                                                                                                                       >  2:742832/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGCTGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGGAACTCACCGGTATGAGTGCCTGGATGGGCCTGGCGGGCTGGGCGATTC‑TGGCA‑CCGCCGGTTGTGCAAGGCTTCTTTTCCATCGGGATTGCCTTTTTGGCCg                                                                                                                                                  >  2:362345/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGGAATTTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTAGGGGTGCCTGGATGGGCATGGCGGGCTGGTCGGTTC‑CGGCA‑CCGCCGCTTATGCAAGGCTTTGTCTCCATCGGGATTGCCTTTTTGGCCg                                                                                                                                                  >  2:1191813/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGCTGATGCAAGGCTTCTTCTCCATCGGGATTGCCTTTATGGCCg                                                                                                                                                  >  1:1191251/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTACTGCCg                                                                                                                                                  >  1:452263/1‑239 (MQ=255)
                                                                                       tCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATAGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGTCTGGTCGCTTC‑TGGCC‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATTGCCTTTATGGcc                                                                                                                                                   >  1:254567/1‑238 (MQ=255)
                                                                                                 aTCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGCTTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATAATTGt                                                                                                                                        >  1:378575/1‑239 (MQ=255)
                                                                                                 aTCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTTTGGCCGTCATCATTGt                                                                                                                                        >  1:746308/1‑239 (MQ=255)
                                                                                                  tCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGCATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCCCCTTTATGGCCGGCATCATTTTa                                                                                                                                       >  1:69244/1‑239 (MQ=255)
                                                                                                        ttGTTCTTCTCTAACGCCACCATCGGCGTTTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGGGCTGGAGATTTGACTGGTGGTGGGGATCAGCGGAATGTTT‑GGTTGCGTCTGGGCGGTGGATCTCACCGGTATGAGTGGCCGGATGGGCATGGCGGCCTGGTGGTTTC‑TGGCA‑CCGCCGATGACGTGAGGCAGCTTCTTCATCAGCATCGCCTTTATGCGCGTCATCATAGTaaatc                                                                                                                                   >  2:227532/2‑235 (MQ=255)
                                                                                                        aTGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGc                                                                                                                                  >  1:793225/1‑239 (MQ=255)
                                                                                                         ttttcttcGCTAACGCCCCCACAGGCGTGTTC‑CCTAACCACTTCGGCGGCTGGCTTGCGGCGCAGAGGTTTTGTCTGGTGTTGGGGTGGATCGAAAGCTTT‑GGTTGCGTCGGGGCGGTGAACCTCCCGGGTAGGAGTGCCTGGATGGGCTGGGCGGACTGGTCGATTC‑TGGCG‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                                 <  2:746308/237‑1 (MQ=255)
                                                                                                          gTTCTCCTCTAACCCCCCCAT‑GGGGGGTTCGGCTCACCGTTGGGGCGGCGGGGGGGGGGCGCGAATGTTTGGGCTGGTGATGGTGATTAGCGAACGCTGT‑GGGGGGGTGGGGGCGTGTAAACACTCCGGTGGGGGTGCCGGGGGGGGTGTGGCGGGGTGTCCTATTC‑TGTCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCAc                                                                                                                                <  2:1099749/239‑1 (MQ=255)
                                                                                                                     aCGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGTTGATCGAAATCTTT‑GGTTGCGTCTGGGCGGGGAAACTCACCGGGATGAGGGCCTGGATGGGCATGGCGGACTGGTCGGTTC‑TGGGC‑ACGCCGGTGATGAAAGGCTTCTTCTCCATCGGGATCGCCTTTATGGCCGTCATCATTGGAATTTCACTGGCCTATATg                                                                                                                     >  2:845292/1‑239 (MQ=255)
                                                                                                                          accaTCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGttct                                                                                                                 >  1:1029899/1‑238 (MQ=255)
                                                                                                                             aTCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGGGATGGGGCTGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGGCGATTC‑TGGCA‑CCGCCGATGGTGCAAGGCTTCTTTTCCATCGGGATCGCCTTTATGGCCGTCATCATTGTAATTGCCCTGGCTtatc                                                                                                                       >  1:626850/1‑228 (MQ=255)
                                                                                                                               cGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTTGCTg                                                                                                           >  1:431068/1‑239 (MQ=255)
                                                                                                                                              ccctctggggcgggTGGGGTGCGGCGGTGAAGTTTTTTGTGGGGATGGGGTTGGTCAAATTCTTT‑GGTGGCGCGGGGGCGGTAAATCTCACGTGTATGGGTGCCTGGAGGGGCTGGGCGGACTGGGCTTTTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACT‑GCTTATATGTTCTTCGCTGGCCGCGCGCTgcgcg                                                                                            <  2:793225/231‑1 (MQ=255)
                                                                                                                                                 aCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGTCTTTATGGCCGTCATCATTGTAATTGCAATGGCCTATTTGTTGTTCCCTGGCCGCGGGCTGCGCGCa                                                                                          >  2:1084938/1‑239 (MQ=255)
                                                                                                                                                   ttCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGAGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCaga                                                                                        >  2:135950/1‑239 (MQ=255)
                                                                                                                                                          gcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGAt                                                                                 >  1:934348/1‑239 (MQ=255)
                                                                                                                                                            gggCGTGGGGCGGTGAAGTTTTGGCGGGTGGTGGGTATGTGCAAATTCTTT‑GGTTGCGTGTGGGGGGGGAAACCCCCCGGTATGGGTGCCTGGATGGGCTTGGCGGACTGGTCGTTTC‑TGCCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATg                                                                                <  2:921117/237‑1 (MQ=255)
                                                                                                                                                             ggCGTGCGGCGCTGAAGTTTTGTCTGGTGCTGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGTCGGGCCGATTC‑TGGCC‑ACGCCGATGATGCAAGGCTTCTTCTCCATCGCGATAGCCTTTCTGGCCGTCATCATGGTAATTGCACTGGCTTATATGGTCTTCGCTGCCCGCGCCCTGCACGCTGAAGAAGCTGcg                                                                              >  2:589577/1‑238 (MQ=255)
                                                                                                                                                                 tcggcggctgagattttGTCTGGTGGGGGGGGTGTTGGAATTGTTG‑GGGTGGGTGTGGGGGGGAAACCGCCCGGGTATGTGTGCCTGGTGGGGCGGGGCGGCCTGGTGTTTTG‑GGCCC‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGaa                                                                           <  2:1009716/224‑1 (MQ=255)
                                                                                                                                                                        cTGAAGATTTTTCGGGTGATGGGGATGACAGAAATCTTT‑GGTGGCGTCTGGGCGGTCAAACTCACCGGTATGAGTGCCGAGAGGGGAATGGCGGACTGGTCTATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaacaac                                                                    <  1:1057385/238‑1 (MQ=255)
                                                                                                                                                                               tttGTCGGGTGATGGGTATGATCGAATTTTTT‑GGTGGCGTCTGGGCGGTGAACCTCACCGGTATGGGTGCCGGGAGGGGCATGGCGGCCTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGa                                                             <  1:292938/238‑1 (MQ=255)
                                                                                                                                                                                   gcgggTGATGGGGATGATAGAATTCTTG‑GGTGGCGTCTGGGGGGTAAAACGCCCCGGTATGAGTGCCGGGATGGGCATGGCGGACTGGTCGATTC‑TGGCG‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                        <  1:1068823/236‑1 (MQ=255)
                                                                                                                                                                                   ctgggTGATGGGGTTTAGGAAAATTTTT‑GGTTGCGTTGGGGCGGTGAACCTCCCCGGTATGGGTGCTTGGATGGGCTTGGCGGACTGGTCTTTTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                        <  1:1191813/236‑1 (MQ=255)
                                                                                                                                                                                    ttGGTGATGGGGATAAGCGAAATCTGT‑GGGGGGGTGGGGGCGGTAAAACTCCCGGGTAGGGGTGCCGGGTGGGGCATGGCGGCCTGTGCTTTTC‑TGGCT‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                        <  2:1149221/237‑1 (MQ=255)
                                                                                                                                                                                     tGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAg                                                         >  2:917006/1‑236 (MQ=255)
                                                                                                                                                                                        cgcacaaaaagatgaaacttttt‑GGTTGTGTGGGGGGGGGGAACCTCCGGGGTATGAGGGGTTGGTGGGGCAGGGCGGTCTGTGCTTTTC‑TGGCC‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGc                                                    <  2:763545/218‑1 (MQ=255)
                                                                                                                                                                                          aTGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTta                                                 >  2:822326/1‑239 (MQ=255)
                                                                                                                                                                                              ggATGATAGAAGTGTTT‑GGTTGCGTCTGGGCGGTGAAGCTCACCGGTAGGAGGGCGGGGAGGGGGAGGGGGGACTGGTCGATTG‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTtaata                                              <  2:828418/238‑1 (MQ=255)
                                                                                                                                                                                                  gatCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAg                                         >  1:460879/1‑239 (MQ=255)
                                                                                                                                                                                                    cctaaaTGTTG‑GGTTGCGTGGGGGGGGTGAAATCCCCGGGTATGGTGGTGTGGATGGTCGTGTGGGCCTGGTGTATCG‑TGACA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                                       <  2:460879/236‑1 (MQ=255)
                                                                                                                                                                                                       gaaTCTTTGGGTGGGGTCGGGGCGGTCAACCTCACCGGTTGGGGTGCCGGGAGGGGAATGGCGTAC‑GGGGGGTTG‑TGGGACCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                                       <  2:378575/236‑1 (MQ=255)
                                                                                                                                                                                                           cTTT‑GGTTGCGTCGGGGGGGTGAAGCTCGCCGTTTGGAGGGTGGGGAGGGGCGTGGCGTCCGGTGCGTTTC‑TGTCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTa                                <  2:1202604/239‑1 (MQ=255)
                                                                                                                                                                                                            tttgggTTGCGCCGGGGAGGTGAAAGTCACCGGTAGGAGGGCCGGGAGGGGCATGGCGGGCTGGTCTATTC‑GGGCA‑CCGCCGATGATGCAAGGCTTCGTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTa                                <  1:378756/235‑1 (MQ=255)
                                                                                                                                                                                                              t‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATg                              >  2:830634/1‑238 (MQ=255)
                                                                                                                                                                                                                 gTTGCGTCTGGCCGGGAAAACTCCCGGGTATGAGTGCTGGGATGGGCATGGCGGACGGGTCAATGC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGAc                            <  1:44204/238‑1 (MQ=255)
                                                                                                                                                                                                                     gtctggggCGGGGAAACCCGCCGTGATGGTTGGCTGGAGGGGTGGGGCGTGCTGCTC‑ATTCGTGGGA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACcgta                        <  2:337543/233‑1 (MQ=255)
                                                                                                                                                                                                                     agTCTGGGCGGCAAAACTCAGCGTATTGGGTGCCGGGATGGAGTTGGGGGCTTGTTCGATTT‑TGTCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACcgta                        <  1:195935/237‑1 (MQ=255)
                                                                                                                                                                                                                           gggggtgaaccccggcGGTATGCGTGCGTGGAGGGGCGTGGCCGACTGGTCTATTC‑TGGCT‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTAtt                  <  2:540393/223‑1 (MQ=255)
                                                                                                                                                                                                                                tGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTAtt                  >  1:520209/1‑233 (MQ=255)
                                                                                                                                                                                                                                tGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTg               >  1:540393/1‑236 (MQ=255)
                                                                                                                                                                                                                                          gggggggAGGGCTTGAATGGGCGGGGGGGACGGGCCGATGG‑TGGAA‑CCGCCGATGATGCAAGGTTTCATCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGGTATGATTATGACCGTACGTATTCTGGCTGTGTGTGGCa  <  1:227532/233‑1 (MQ=255)
                                                                                                                                                                                                                                               tGAGTGCCTGGATGGGCATGGCTGACTGGTCGATTC‑TGGCA‑CCGCCGATGGTGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAACACAACTGGCAGAACAGTCTGCTtaa                                                >  1:649853/1‑188 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         
CGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCC‑GATGTTCTTCTCTAACGCCACCATCGGCGTGTTC‑GCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTT‑GGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTC‑TGGCA‑CCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTGGCA  >  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: ‑