Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,481 A→G A262A (GCA→GCG creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,4810AG98.3% 190.3 / ‑4.5 61A262A (GCA→GCGcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (25/34);  total (25/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 = 7.42e-01

CGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGT  >  NC_000913/4637235‑4637723
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
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)
                                                                   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)
                                                                                                             gccgATTCCGTCACTGACAATAAGCCCTGTCCTCTTCCTGAACTCGGCATTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATAGTGAGCAGTATGTTTATGCGTTACCCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCa                                                                                                                                   <  1:565372/248‑1 (MQ=255)
                                                                                                                 aTTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAAt                                                                                                                               >  1:133329/1‑251 (MQ=255)
                                                                                                                 aTTCCGTCACCGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCTGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAAt                                                                                                                               <  1:398753/251‑1 (MQ=255)
                                                                                                                       tCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGCCGCAAAATGCGCGt                                                                                                                         >  1:318241/1‑251 (MQ=255)
                                                                                                                         ccTGACAATAAGCCCTTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGATGCAAAATGCGCGTAt                                                                                                                       <  2:350888/250‑1 (MQ=255)
                                                                                                                             tcAATAAGCCCGTTCCTATTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAACACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGCAGGTCCGCCACCGGAAGTAGGGGGTCGGTTTCCCGACCACTTTCTGACGcctattgccagcttgt                                                                                                                     >  2:149511/2‑234 (MQ=255)
                                                                                                                              cAATAAGCCCGTTCCTCTTCCTGAACGCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAgg                                                                                                                  <  1:173072/251‑1 (MQ=255)
                                                                                                                                  aaGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAggg                                                                                                                                                                   >  1:347867/1‑196 (MQ=255)
                                                                                                                                  aaGCCCGTTCCTCTTCCTGAACACGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAgtg                                                                                                                                                                   <  2:347867/198‑1 (MQ=255)
                                                                                                                                        gTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGa                                                                                                        <  2:133329/251‑1 (MQ=255)
                                                                                                                                           cctcttccggaacgtggcagtagCGAGTGGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGGAGGTTTATGCGGTAACCCATGATATAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTGTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGTCATTGGTAGAAAc                                                                                                     <  2:647997/233‑1 (MQ=255)
                                                                                                                                                     aactcGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGAAGAAACGTTACTACGc                                                                                           >  1:545939/3‑251 (MQ=255)
                                                                                                                                                      actcGGCAGTAGCGAGTTGCGTAATCTTGCGCAGGCTCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCTTGGCGCGGCGGAAATTTTACGCGAATGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTTCTACGcc                                                                                          <  1:203695/250‑1 (MQ=255)
                                                                                                                                                                  cGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGCGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTAAGCGAAGGTCCGCCAACGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTTGTAGCAAAGTTCCTACACCAGGCAAGAcctg                                                                              >  2:306957/1‑248 (MQ=255)
                                                                                                                                                                            aagcttgcgcAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCa                                                                    <  1:130575/245‑1 (MQ=255)
                                                                                                                                                                             agcttgcgcAGGCGCTGGAAAGTATGCGCGGGAAGCGGGAAGGGAAAAACTATATTGAGCAGGATGTTTATGCGTTAACCCATGATCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGAGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACGGTAGAATCTTCAg                                                                   <  2:237566/246‑1 (MQ=255)
                                                                                                                                                                                gtgagcAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATGGAGCAGTATGTTTATGCGTTAACCCATGAGCGAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGAAAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGaa                                                                <  1:353859/247‑1 (MQ=255)
                                                                                                                                                                                      aGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGtt                                                          <  1:14477/251‑1 (MQ=255)
                                                                                                                                                                                       ggCGCTGGAAAGTATGCGGGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTCTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGTGATTCGTGGCGAGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCGCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGTCAAGCCTGGAGAATCGTCAGGAAGTCGTTc                                                         <  1:672853/251‑1 (MQ=255)
                                                                                                                                                                                            gggAAAGTGTGCGAGTTAAGCGGGAAGGGAAAAATCATATTGAGCAGTATGTTGATTCGTTAACCCATGGGTTAAAACGCCCACTGGCGGGGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAATGGGTGGTTCGTTTGATCGACAACAGTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAGACTTTACTACGGCAGGCAAGACTGGTGAATCGTCAGGAAGTCGTTCGGAct                                                    <  2:596264/250‑1 (MQ=255)
                                                                                                                                                                                               aaaGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGActgtt                                                 <  2:676809/251‑3 (MQ=255)
                                                                                                                                                                                                aaGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGActgttg                                                <  1:114347/251‑4 (MQ=255)
                                                                                                                                                                                                       cgcgTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGGGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCGGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGAttg                                         >  2:331274/1‑249 (MQ=255)
                                                                                                                                                                                                         cgTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGAATGGAGAATCGTCAGGAAGTCGTTCTGACTgttgtt                                              >  1:383989/1‑244 (MQ=255)
                                                                                                                                                                                                         cgTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGGCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGGTCTGACTgttgtt                                              <  2:383989/244‑1 (MQ=255)
                                                                                                                                                                                                            gAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGAggcgg                                    >  1:88346/1‑251 (MQ=255)
                                                                                                                                                                                                             aaGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTAAGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGAATGTTGTGTATGTggcggg                                   >  2:10495/1‑250 (MQ=255)
                                                                                                                                                                                                              aGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCa                                  <  2:226275/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCattat                              <  1:149511/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTtcgc                         >  1:273418/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ttGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCgg              >  1:410181/1‑249 (MQ=255)
                                                                                                                                                                                                                                     gAGCAGTATGTTTATGCCTTAACCCATGAGCTAAAAAGCCCCCTGGCGGCGATTCGTGGCGCTGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAgc           <  2:151128/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCa      <  2:410181/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCAc     >  1:119926/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGGGGCTCGTTTTACCGCCAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGGCGTTCTGACTGTTGTTGTTGTGGCGGCCTTTTTTCGCCGCGTcctgcaacgcggcaccgg  >  2:590671/1‑234 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
CGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGT  >  NC_000913/4637235‑4637723

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