Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,346 T→C A217A (GCT→GCC creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,3460TC98.6% 213.8 / ‑5.2 70A217A (GCT→GCCcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (33/35);  total (34/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.09e-01

GCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACAT  >  NC_000913/4637099‑4637585
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
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                                                                                              >  2:558367/1‑221 (MQ=255)
                                                                                                                                                                              ggATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAg                                                                                              <  1:558367/221‑1 (MQ=255)
                                                                                                                                                                                 ttGCACTGGTGAGTGGCGCGGGCGTGGTTTGTTGGAGCACACTCGCTATGGCCCGTCTCCCGCGCGATTGTGATTCCGTCACTGACAAGAAGCCAGTTCCTCTGCCAGATCTCGGTAGTAGCGGGTTGCGTAAAAGTGCGCAGGCGCTGGAAAGTATGCGCGGGAAGCTGGAAGGGAAAACCTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGTCCACTGTCGGCGATTCGTGGAgc                                                              <  2:511138/250‑1 (MQ=255)
                                                                                                                                                                                                   cAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGc                                                                                                                                                                                              <  2:117255/104‑1 (MQ=255)
                                                                                                                                                                                                   cAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGc                                                                                                                                                                                              >  1:117255/1‑104 (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)
                                                                                                                                                                                                                               caTTGCCCGGTTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGtt               <  2:429031/250‑1 (MQ=255)
                                                                                                                                                                                                                                aTTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGttt              >  1:119145/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTAccg         <  1:196803/251‑3 (MQ=255)
                                                                                                                                                                                                                                     ccGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGAATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTAccg         >  1:496484/1‑249 (MQ=255)
                                                                                                                                                                                                                                      cGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTAccg         >  1:647997/1‑248 (MQ=255)
                                                                                                                                                                                                                                            aCTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACAt  <  1:89234/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
GCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACAT  >  NC_000913/4637099‑4637585

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: