Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,322 T→C S209S (TCT→TCC creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,3220TC98.6% 223.5 / ‑5.5 71S209S (TCT→TCCcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (35/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCT  >  NC_000913/4637078‑4637568
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
gTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGcc                                                                                                                                                                                                                                                  <  1:121264/251‑1 (MQ=255)
 tgggcAGTATGGTGCGCGCAGCACGTTGCCAAATCCTGCAGCTCCCGCACGTTCGGGGATGTATGTTGCCGCGCCGATTCTGGACGGCGCGCGGCTTCTTGGCGTTTTGAGCGGCGGCCAACCGAACGCGGCGATGGCTCCGGTCCTTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCAATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCGCCATTGccc                                                                                                                                                                                                                                                 >  2:130575/5‑251 (MQ=255)
     cAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCt                                                                                                                                                                                                                                             <  1:278620/251‑1 (MQ=255)
         aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                                                                                                                                                                                                                         >  1:151128/1‑251 (MQ=255)
         aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                                                                                                                                                                                                                         >  1:426605/1‑251 (MQ=255)
         aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                                                                                                                                                                                                                         >  1:226275/1‑251 (MQ=255)
               cgcgacgcaagTTGAAAAAGCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCTGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTAt                                                                                                                                                                                                                                   <  2:233709/241‑1 (MQ=255)
                     gcaCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGAt                                                                                                                                                                                                                             >  1:429031/1‑251 (MQ=255)
                      caCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGCCGTTTTGATCGGAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCTGCGCCGTTTTGCTGGGGATTGCACTGGTGATTGGCGCAGTCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGAtt                                                                                                                                                                                                                            <  1:295852/251‑1 (MQ=255)
                            gCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTc                                                                                                                                                                                                                       <  1:30617/250‑1 (MQ=255)
                                      gCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCCGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATa                                                                                                                                                                                                             >  1:77052/1‑250 (MQ=255)
                                        cGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTAGCTATGCCGATTCCGTCACTGACACTAAGc                                                                                                                                                                                                          >  2:196803/1‑251 (MQ=255)
                                          aTCCCGAAAGTTCGGTGCTGTATGTTGCCTCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGGGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGGGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGc                                                                                                                                                                                                          <  2:116988/249‑1 (MQ=255)
                                          aTCCCGAAAGTTCGGTGATGTCTGTTGCCGAGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTACGTCACTGACAATAAGccc                                                                                                                                                                                                        >  1:67139/1‑251 (MQ=255)
                                                aaaGTTCGGTGTTGTATGTTGCCGCGCCGATTATGGACGGCTCGAGGCTTATTGGCGTTTTGAGCGGAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCct                                                                                                                                                                                                  <  2:529102/251‑1 (MQ=255)
                                                 aaGTTCGGAGATGTATGTTGCCGCGCCGATTATGGACGACTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctc                                                                                                                                                                                                 >  1:174373/1‑251 (MQ=255)
                                                       ggtgATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcttcctg                                                                                                                                                                                           >  1:130186/2‑246 (MQ=255)
                                                                      agcgccGATTATGGACGGCTCGCGGCTTATTGACGTTTTGAGCGTAGGCAAAACGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGAGATTGGCGCAGGCATGGTTTGGAGGATCAACCGATCCATTGCCCGGCTCACTCGCTTAGCCGCTTCCGTCACTGGCAATAAGCCCGTTCctcttcctgaaatcggccgtagcg                                                                                                                                                                            >  2:546005/5‑231 (MQ=255)
                                                                            gATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGCTTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGc                                                                                                                                                                      >  2:287923/1‑251 (MQ=255)
                                                                             aTTATGGACGGCTCGCGGCTTATTGGCGATTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGACATTTTGCTGGGGATTGCACTGGGGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCTCCATTGACCGGCTCACTCGCTATACCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGACCTCGGCAGTAGCGAGTTGCg                                                                                                                                                                     >  1:441493/1‑251 (MQ=255)
                                                                                aTGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCCAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGCTTGGTGGATCACCCGCTCCCTTGCCCGGCTCACACGCTATGCCGATTCCGTCACTGACAATAAGCCCGCTCCTCTTCCTGAACTCGGCAGTAGCGTGTTGCGTa                                                                                                                                                                   >  2:173072/1‑250 (MQ=255)
                                                                                     cGGCGCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCACCGGTCATTAAGCGTAGCGCGCGGCGATTTTTTTGGGCCCGCGCCATTTTGCTGGGGATTGAACTGGTGGTTGGCGCAGGCCTGGTTTGTTGGATCAACCGCTCCATTGCCCGGCTCACTCTCTATGTCGGTTCCGTCACCGACAATCAGCCCGTTCCTCTTCCTGAACTCCGCAGTATCGGGTTGAGTAAACttt                                                                                                                                                             >  2:491606/1‑249 (MQ=255)
                                                                                      ggCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGGGATTGGCGCGGGCATGGTTTGGTGGATCAACCGCTCCATTGCACGGCTCCCTCGCAATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGACCTCGGCAGTAGCGGGTTGCGTaagctggc                                                                                                                                                            >  2:178703/1‑245 (MQ=255)
                                                                                        cTCGCGGCTTATTGGCGTTTTTAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGAGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTACCTTGCCCGTCTCACTCGCTATGCCGGTTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGACACCGGCAGTAGCGAGTTgagatagcttggtc                                                                                                                                                          >  2:366674/1‑238 (MQ=255)
                                                                                         gcgcgTCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCa                                                                                                                                                         <  2:229176/250‑1 (MQ=255)
                                                                                          cgcgGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGAGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGCACTCGGCAGTAGCGAGTTGCGTaagcttgagc                                                                                                                                                          >  2:264149/1‑241 (MQ=255)
                                                                                            cgGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGgc                                                                                                                                                      <  1:4276/251‑1 (MQ=255)
                                                                                               cTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCt                                                                                                                                                   >  1:350888/1‑251 (MQ=255)
                                                                                                 tATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTgg                                                                                                                                                 >  1:526041/1‑251 (MQ=255)
                                                                                                     ggCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAg                                                                                                                                             <  1:585194/251‑1 (MQ=255)
                                                                                                       cGTTTTGAGTGTAGGCAAACCGACCGCGGCGGTGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTa                                                                                                                                           <  1:434781/251‑1 (MQ=255)
                                                                                                         ttttGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcttcctgaactc                                                                                                                                                                                      >  1:149106/1‑196 (MQ=255)
                                                                                                         ttttGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcgtcctgaactc                                                                                                                                                                                      <  2:149106/206‑12 (MQ=255)
                                                                                                          tttGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATg                                                                                                                                         <  2:217120/250‑1 (MQ=255)
                                                                                                             gAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATgcgc                                                                                                                                      <  1:337048/250‑1 (MQ=255)
                                                                                                               gCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTg                                                                                                                                   <  1:287923/251‑1 (MQ=255)
                                                                                                                         aCCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGa                                                                                                                           >  1:200249/1‑249 (MQ=255)
                                                                                                                          ccGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAgg                                                                                                                        <  2:600317/251‑1 (MQ=255)
                                                                                                                          ccGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCATGAAGCTGGAAgg                                                                                                                        <  1:340328/251‑1 (MQ=255)
                                                                                                                               cgcgGCGATGGCACCGGTCATTAAACGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCCTTTTGCTGGGGATTGCACTGGCGATTGGCGCAGGCATGGTTTGGTGGCTCCACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCCCTGCCAATAAGCCCGTTCCTCTTCATGAACTCGGCAGTCGCGATTTGCATAAGCTTGCGCAGGCGCTGGCAAGTATGCGCGTGAAGCTGGAAgggcaa                                                                                                                    >  2:599127/1‑247 (MQ=255)
                                                                                                                                 gggcgATGGCGCCGGGCATTAAGCGTAGCGAGCGGCGAAGTTTATGGGCCAGCGCCATTTTGCTGGGGAGTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGTCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAAc                                                                                                                 <  2:283024/250‑1 (MQ=255)
                                                                                                                                  ggcgATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGATGGAAGGGAAAAACt                                                                                                                >  2:555904/1‑251 (MQ=255)
                                                                                                                                    cgATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTGGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACtat                                                                                                              <  1:267653/251‑1 (MQ=255)
                                                                                                                                        gtgtCCGGTCATTAAGCGTAGCGAGCGGCGATTTTTATGGGCCAGCGCCATTTGGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTg                                                                                                          <  2:174373/248‑1 (MQ=255)
                                                                                                                                        ggCTCCGGTCATTAATCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTg                                                                                                          <  2:335950/251‑1 (MQ=255)
                                                                                                                                          cTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTCTGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCCTGGTTTGGTGGCTCAACCGCTCCATTGCCCGGCTCACTCGCTCTGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGGGTTGCGTCAGCTTGCGCCAGCGCTGGACAGTATTCGCGTGAAGCTGGGAGGGAAAAACTATATTg                                                                                                          >  2:114347/1‑249 (MQ=255)
                                                                                                                                             cGGTCATTACGCGTAGCGAGCGGGGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACCGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGACCGTTCCTCTTCCTGAACTCGGCAGTAGCGCGTTGCGTAAGCTTGAGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGATACACTATATTGAGc                                                                                                       >  2:628074/1‑249 (MQ=255)
                                                                                                                                                    gaaGCGGAGCGAGGGTCGAATTGTAGGGGCCAGCGCTAGTTTGCTGGGGATGGCACTGGTGAGTGGCGCAGGCATGGGTTGGTGGAGCATCCGCGCCAGTGCCCGGCTCACTCGCTATGCCGAGGCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACACGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGCAAATTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGATCAGTATGtt                                                                                               <  2:67139/249‑1 (MQ=255)
                                                                                                                                                           agcgagcgGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTatgcgt                                                                                        <  1:555904/250‑6 (MQ=255)
                                                                                                                                                             cgagcgGCGAATTTTATGGGCCATCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTAGGCCGATTTCGTCACGGACAATAAGCCTGTTTCTCTTCCTGAACTCGGCAGTAGCGAGGTGCGTAAGCTTGCGCAGGAGCTGGAAAGGATGCGCGTGAAGATGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTaa                                                                                     <  2:77052/251‑1 (MQ=255)
                                                                                                                                                               agcgGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGAGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGCAAGAATGCGCGTGAAGCTGGAAGGGAAAAACTAAATTGAGCAGTATGTTTAGGCGTTAAcc                                                                                   >  2:604786/1‑250 (MQ=255)
                                                                                                                                                                  ggcgAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCAt                                                                                >  1:411026/1‑251 (MQ=255)
                                                                                                                                                                     gAATTTTATGGGCCAGCGCCAGTTTGCTGGGGATTGCACTGGGGGTTGGCGCAGGCATGGTTTGGTTGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCGCTGACAGTAAGCCCGTTCCTCGTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAg                                                                             <  2:183838/251‑1 (MQ=255)
                                                                                                                                                                       aTTTTATGGGCCATCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCt                                                                           <  2:335265/251‑1 (MQ=255)
                                                                                                                                                                          ttATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTaa                                                                         <  1:270227/250‑1 (MQ=255)
                                                                                                                                                                              gggCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAg                                                                     <  1:604786/250‑1 (MQ=255)
                                                                                                                                                                                     gtccATTTTTGTGGGGATTGCACTGGTGATGTGCGCAGGCATGGTGTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCGATGCCGATTCCGTCACTGTCAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTgg                                                             <  2:196753/249‑1 (MQ=255)
                                                                                                                                                                                         aTTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGAGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTggaggc                                                         >  1:483267/1‑251 (MQ=255)
                                                                                                                                                                                              gctgGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTc                                                    <  1:491606/249‑1 (MQ=255)
                                                                                                                                                                                              gctagcgagtGCACTGGTGATTGGGGCAGTTATGTTTTGGTGGATCAACCGCTCCATTGCCCGGCTCTCTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTGCCTCTTCCTGAACTCTGCAGTAGCGAGTTGGGTAAGCTTGCGCATGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGCGCTAAAAAGCCCGCTGGCGGCGATTc                                                    <  2:426605/242‑1 (MQ=255)
                                                                                                                                                                                                   ggATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAg                                                                             <  1:558367/221‑1 (MQ=255)
                                                                                                                                                                                                   ggATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAg                                                                             >  2:558367/1‑221 (MQ=255)
                                                                                                                                                                                                      ttGCACTGGTGAGTGGCGCGGGCGTGGTTTGTTGGAGCACACTCGCTATGGCCCGTCTCCCGCGCGATTGTGATTCCGTCACTGACAAGAAGCCAGTTCCTCTGCCAGATCTCGGTAGTAGCGGGTTGCGTAAAAGTGCGCAGGCGCTGGAAAGTATGCGCGGGAAGCTGGAAGGGAAAACCTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGTCCACTGTCGGCGATTCGTGGAgc                                             <  2:511138/250‑1 (MQ=255)
                                                                                                                                                                                                                        cAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGc                                                                                                                                                                             >  1:117255/1‑104 (MQ=255)
                                                                                                                                                                                                                        cAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGc                                                                                                                                                                             <  2:117255/104‑1 (MQ=255)
                                                                                                                                                                                                                                gTTTGGTGGATCAACCTCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGGAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTccgcc                  <  1:126654/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggATCAACAGCACCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGAGAAGCTGGACGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGCAAttttcagggacggaccgcccccgga            >  2:89234/1‑229 (MQ=255)
                                                                                                                                                                                                                                           cAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACt                                                               >  2:505209/1‑195 (MQ=255)
                                                                                                                                                                                                                                           cAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACt                                                               <  1:505209/195‑1 (MQ=255)
                                                                                                                                                                                                                                              ccGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCAATGGCGGCGATTCGTGGCGCGGCGGACATTTTACGCGGCAGTCCGCCACCGGAAGgggtgg    >  1:403852/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gCTCCATTGCCCGGCTCACTCGGTATGCAGATTCTGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGAGTGAATTTTACGCGAAGGAccggcccagggatgggtgggg  >  2:146076/1‑233 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
GTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCT  >  NC_000913/4637078‑4637568

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