Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,151 A→G A152A (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,1510AG98.6% 216.3 / ‑4.8 70A152A (GCA→GCG) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (36/32);  total (37/33)
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.66e-01

TAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCG  >  NC_000913/4636908‑4637383
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
tAATCCGCTACAACATATCCCGTTTCGCGCCAATATCGGTGGTATTAACAAAGTTCGCATTGA‑TATTCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATCACGTCTGCCTAACGTTGCGTGGTCAGTATGTTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGCCCGGTGATGTATGTTGCCGCTCCGAtt                                                                                                                                                                                                                                     <  2:97660/250‑1 (MQ=255)
tAATCAGCTACAACATCGCCCGTTTCGCGCCACTATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAACGCAGTTGGGCAGGATTATTCGCGCTGGAACGACGTCTGGCTAACGTTGCGTGGTCAGTTTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCTGAAAGTTCGGTGATGTATGGTGCCGCGCCggtt                                                                                                                                                                                                                                     >  2:270227/1‑247 (MQ=255)
    cAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGa                                                                                                                                                                                                                                <  1:134192/251‑1 (MQ=255)
          caacaTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTc                                                                                                                                                                                                                          <  1:304694/251‑1 (MQ=255)
             caTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTTCCGCGCCGATTATGGACGGCTcgcg                                                                                                                                                                                                                       >  2:260222/1‑251 (MQ=255)
             caTCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTcgcg                                                                                                                                                                                                                       >  1:508434/1‑251 (MQ=255)
                  cgcttttCGCGCCATTTGCGGTTGCATTAACAAAGTGCGCAATGAATAT‑CATGTCAATATGACCGATTCGCAGGTTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGGGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTa                                                                                                                                                                                                                  <  2:656271/247‑1 (MQ=255)
                  ccctgttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAAGAAAGACGTTGGGCAGGATTTTTCGCGCTGGAATGTCGGCTGGTTAACGTTGCGTGGTCCGTATGGTGCGCGCCGCACGTTGGAAAATCCTGCCGATCCCGAAAGTGCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTa                                                                                                                                                                                                                  <  1:594315/246‑1 (MQ=255)
                      tttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTgg                                                                                                                                                                                                              >  1:335265/1‑251 (MQ=255)
                           cgcCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGttt                                                                                                                                                                                                         >  1:14871/1‑251 (MQ=255)
                              cAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTg                                                                                                                                                                                                       >  1:229176/1‑250 (MQ=255)
                                                aaaaGTGCGAAAGGAATGT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCGGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTGGCAAACCCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAAcg                                                                                                                                                                                    <  2:590219/250‑1 (MQ=255)
                                                 aaaGTGCGCAATGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCACACCGAAcgc                                                                                                                                                                                   >  2:121264/1‑251 (MQ=255)
                                                        gctaTGAATAT‑CATGTCTATATGACCGATTCGGAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTGGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATg                                                                                                                                                                            <  2:255465/248‑1 (MQ=255)
                                                          aaTGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑c                                                                                                                                                                         >  1:196753/1‑251 (MQ=255)
                                                          aaTGAATAT‑CATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑c                                                                                                                                                                         >  1:600317/1‑251 (MQ=255)
                                                                 aatcctgtctatctgctcgtttcgcAGGGTAAAGTATTGCTCGATTCGGAAAAGAATGCCGTTGGGAAGGTGTTTTCGCGCTGGAATGATGACTGTCTAACGGTGCGTGGTGAGTAGGGTGCTCGCAGCATGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGGAGGTTGCCGCGCCGATTATGGACGGGGCGCGGCTGATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGT‑CTCCGGt                                                                                                                                                                   <  2:602867/229‑1 (MQ=255)
                                                                    cATGTCTATATGACTGATTCGCAGGTTAAAGTTTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGTCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCa                                                                                                                                                                 <  2:14871/250‑1 (MQ=255)
                                                                                  cGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTagcgagc                                                                                                                                                  >  1:183838/1‑251 (MQ=255)
                                                                                    attcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGCTGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGc                                                                                                                                                                                                             >  2:82391/4‑191 (MQ=255)
                                                                                    attcgcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGc                                                                                                                                                                                                             <  1:82391/188‑1 (MQ=255)
                                                                                    agtcgcAGGGTAACGTATTGTGCGATGCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCTCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAAGCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGGAGCGAgcgg                                                                                                                                                <  2:126663/248‑1 (MQ=255)
                                                                                        gcAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAAGCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGaa                                                                                                                                            <  1:620334/251‑1 (MQ=255)
                                                                                               aaaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATg                                                                                                                                     <  1:580329/251‑1 (MQ=255)
                                                                                                aaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATgg                                                                                                                                    <  1:260222/251‑1 (MQ=255)
                                                                                                  gTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTATGGCTAACGTTGAGTGGACAGTATGGTGCGAGCAGCACGTAGCAAAATCCTGCCGATCCCGAAAGTGCGGAGATGTATGTTGACGCGCAGATTATGGACGGCTCGCGGCTTATTGGAGTTTTGAGCGTAGGCAAACCGAAAGCGGCG‑TGGACTCCGGGCCTTAAGCGTCGCGAGAGGCGGATTTTATgggg                                                                                                                                  >  2:203695/1‑250 (MQ=255)
                                                                                                   tATTGTTCGATGCGGCAAATAAAACCGTTGGGCAGGATTATTCGCGCTGGACTGACGTCTGGATAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGACAGTTCGGTGATGTATGTTGCCGAGACGATTATGGACGGCCCGCGGCTTAATGGCGTTTTGATCGTAGGCTACACGCACGCGGAGCAGG‑CTACGGTCATGACGCGTAGCGAGCCGCGAATTTTATTGGcc                                                                                                                                 >  2:126654/1‑251 (MQ=255)
                                                                                                    aTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGGTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCTTTAAGCGTAGCGAGCGGCGAATTTTAGGGGcc                                                                                                                                 >  2:337048/1‑250 (MQ=255)
                                                                                                      tgttcctttCGGAAAAAAAAGCCGTGGGGCAGGATTATTCGCGCTGGAATAACTTCTGTCTAACGTTGCTTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAg                                                                                                                               <  2:508434/243‑1 (MQ=255)
                                                                                                          cGATTCGGCAAATAAAGACGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGAAAACTCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCACGGCTTATTGGCGTTTTTAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAgcaccc                                                                                                                          >  2:585194/1‑247 (MQ=255)
                                                                                                              tCGGCAAATAAAGCCGTTGGGCAGGATTACTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTGGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAttt                                                                                                                       <  1:120465/250‑1 (MQ=255)
                                                                                                               cGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttg                                                                                                                     <  1:229325/251‑1 (MQ=255)
                                                                                                               cGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttg                                                                                                                     >  1:335950/1‑251 (MQ=255)
                                                                                                                 gCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTttgct                                                                                                                   >  1:310793/1‑249 (MQ=255)
                                                                                                                      tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTgggg                                                                                                               >  1:673389/1‑250 (MQ=255)
                                                                                                                      tAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGAATATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGAAGGCCAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTAGTGGCAAGCGCCATTTTGCTggggg                                                                                                              >  2:68293/1‑250 (MQ=255)
                                                                                                                         ggCCGTTGGGCAGGATTATTCGGGTTGGAATGAGGTCTGGCTAGCGTTGCGTGGTTAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTg                                                                                                           <  2:183773/250‑1 (MQ=255)
                                                                                                                         aGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGAATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGAGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTg                                                                                                           >  1:246817/1‑251 (MQ=255)
                                                                                                                           ccGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGc                                                                                                          >  1:116988/1‑250 (MQ=255)
                                                                                                                                gggcAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGt                                                                                                    >  1:676809/4‑251 (MQ=255)
                                                                                                                                     ggATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTg                                                                                               >  1:283024/1‑251 (MQ=255)
                                                                                                                                        ttattCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTTTGGGCCAGCGCCATTTTGCTGGGGTTTGCACTGGTGATTgtcg                                                                                            >  2:340328/1‑248 (MQ=255)
                                                                                                                                            tCGCGCTGGAATGACGTCTGGCTCACGTTGCGTGGTCAGTATGGTGCGCGAAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCACGCGGATTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCCTTAAGCGTAGCGAGCGTCGATTTTTATGGGCCAGCGCCCTTTTGCAGGGGATTGACCGGGTGATTGGCGCAgg                                                                                        >  2:631409/1‑251 (MQ=255)
                                                                                                                                                  tGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGAACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTAGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGt                                                                                  >  1:76117/1‑251 (MQ=255)
                                                                                                                                                   ggAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATgg                                                                                   <  1:181025/249‑1 (MQ=255)
                                                                                                                                                    gaatgactttgtgcataCGTTGCGGGGGCAGTTTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGttt                                                                                <  2:489152/235‑1 (MQ=255)
                                                                                                                                                    gAATGACGTCTGTCTAACGTTGCGTGGGCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGGTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGttt                                                                                <  1:631409/251‑1 (MQ=255)
                                                                                                                                                     aaTGACGGCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTtg                                                                               <  1:572861/251‑1 (MQ=255)
                                                                                                                                                        gACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTAATTGGCCCAGGCATGGTTtggtg                                                                            >  1:237566/1‑251 (MQ=255)
                                                                                                                                                                gttAACGTGGAGTGGTCAGTATGGTGCGCGAGGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAAcc                                                                    <  2:673389/249‑1 (MQ=255)
                                                                                                                                                                     cgtggcgtgtcaaGTATGGGGCGCGCAGCCCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGGTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATGTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGctcc                                                               <  2:310793/239‑2 (MQ=255)
                                                                                                                                                                           gTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGcc                                                         <  1:121264/251‑1 (MQ=255)
                                                                                                                                                                            tgggcAGTATGGTGCGCGCAGCACGTTGCCAAATCCTGCAGCTCCCGCACGTTCGGGGATGTATGTTGCCGCGCCGATTCTGGACGGCGCGCGGCTTCTTGGCGTTTTGAGCGGCGGCCAACCGAACGCGGCGATGG‑CTCCGGTCCTTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCAATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCGCCATTGccc                                                        >  2:130575/5‑251 (MQ=255)
                                                                                                                                                                                cAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGTTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGAt                                                                                                                                                                             <  1:399393/131‑1 (MQ=255)
                                                                                                                                                                                cAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCt                                                    <  1:278620/251‑1 (MQ=255)
                                                                                                                                                                                    aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                                >  1:226275/1‑251 (MQ=255)
                                                                                                                                                                                    aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                                >  1:151128/1‑251 (MQ=255)
                                                                                                                                                                                    aTGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACt                                                >  1:426605/1‑251 (MQ=255)
                                                                                                                                                                                          cgcgacgcaagTTGAAAAAGCCTGCCGATCCCGAAAGTTCGGGGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCTGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTAt                                          <  2:233709/241‑1 (MQ=255)
                                                                                                                                                                                                gcaCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGAt                                    >  1:429031/1‑251 (MQ=255)
                                                                                                                                                                                                 caCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGCCGTTTTGATCGGAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCTGCGCCGTTTTGCTGGGGATTGCACTGGTGATTGGCGCAGTCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGAtt                                   <  1:295852/251‑1 (MQ=255)
                                                                                                                                                                                                       gCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTc                              <  1:30617/250‑1 (MQ=255)
                                                                                                                                                                                                                 gCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCCGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATa                    >  1:77052/1‑250 (MQ=255)
                                                                                                                                                                                                                   cGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTAGCTATGCCGATTCCGTCACTGACACTAAGc                 >  2:196803/1‑251 (MQ=255)
                                                                                                                                                                                                                     aTCCCGAAAGTTCGGTGCTGTATGTTGCCTCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGGGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGGGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGc                 <  2:116988/249‑1 (MQ=255)
                                                                                                                                                                                                                     aTCCCGAAAGTTCGGTGATGTCTGTTGCCGAGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTACGTCACTGACAATAAGccc               >  1:67139/1‑251 (MQ=255)
                                                                                                                                                                                                                     aTCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAAc                                                                     <  1:261896/197‑1 (MQ=255)
                                                                                                                                                                                                                           aaaGTTCGGTGTTGTATGTTGCCGCGCCGATTATGGACGGCTCGAGGCTTATTGGCGTTTTGAGCGGAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCct         <  2:529102/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaGTTCGGAGATGTATGTTGCCGCGCCGATTATGGACGACTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctc        >  1:174373/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggtgATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcttcctg  >  1:130186/2‑246 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         
TAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAT‑CATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGG‑CTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCG  >  NC_000913/4636908‑4637383

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: ‑