Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,925 G→C R77P (CGC→CCC)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,9250GC94.9% 184.5 / 3.2 60R77P (CGC→CCC) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base C (26/30);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

GCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCG  >  NC_000913/4636680‑4637171
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGt                                                                                                                                                                                                                                                   <  1:65134/251‑1 (MQ=255)
  ctgcggggCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGATGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCCTGGGCAACTGGCGCCGGCGTTTAATCAGCTCCAACATCCCCCggtt                                                                                                                                                                                                                                                 >  1:332292/5‑248 (MQ=255)
        ggCCTGTAATGCGTATCGGCATGCGGTTGTTGGTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGGCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGCCGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGcc                                                                                                                                                                                                                                           <  2:406125/251‑1 (MQ=255)
         gtctgtaatttggaTCGGCAGGCGGTTGTTTCTGGGCTATTTTTGCCTGGTGGTGGTGGCAGCCTGGTGCGTACTGGCCATTTTTGGCAAAGAAGTTAAACTGGGCGTGCGAAGAGCAACGGAGGGGCTGTTGATCGTCTCCGCAACGTTGCTGACGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGCCCCAACGCATGGGCAACTGTCGCAGGCGTTTAATCAGCTTCATCATCCCCCGTTTcgcgtca                                                                                                                                                                                                                                          <  2:262240/238‑4 (MQ=255)
             gTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGCAGGGCCACGGAGGGGCCGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGCTTTGCTCTCTGGGGACCCAACGCCTGGGCAACTGGCGCAGGCGTTTAATCAACTACAACAACCCCCGTTTCGCGCCAatat                                                                                                                                                                                                                                      >  1:436514/1‑251 (MQ=255)
             gTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAatat                                                                                                                                                                                                                                      >  1:399158/1‑251 (MQ=255)
              tAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATc                                                                                                                                                                                                                                     >  1:221315/1‑251 (MQ=255)
               aaTGCGTATAGGCAGGCGGTTGTTGCTGGGCTATTTTTTACTTGTGGCGGTGGCAGTCTGGTTCTTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCg                                                                                                                                                                                                                                    <  2:259695/251‑1 (MQ=255)
                      aTCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGGACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCAt                                                                                                                                                                                                                             <  1:104322/251‑1 (MQ=255)
                          gCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAAc                                                                                                                                                                                                                         <  1:394330/251‑1 (MQ=255)
                                            gCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTTATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAatat                                                                                                                                                                                                       <  1:1692/251‑1 (MQ=255)
                                             cTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATc                                                                                                                                                                                                      >  1:173566/1‑251 (MQ=255)
                                               aTTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAt                                                                                                                                                                                                    <  1:151620/251‑1 (MQ=255)
                                                  ttttACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGt                                                                                                                                                                                                  <  1:363252/250‑1 (MQ=255)
                                                      aCTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATGTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCtat                                                                                                                                                                                              <  2:41162/250‑1 (MQ=255)
                                                           tggcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAc                                                                                                                                                                                         <  1:426394/250‑1 (MQ=255)
                                                            ggcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACCCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCg                                                                                                                                                                                       <  1:455230/251‑1 (MQ=255)
                                                                 tgacAGCCTGCTTCGTACTGGCCCTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGAAACCGCAACGTTGCTGGCGGAGCTGGCGCGCCCCGATTTGCTCTCTGGGGACCCAACGCATGGGACACTGGCGCAGGCGTTTAATCAGCTACAACATCACCCGTTTCGCGCCAATATAGGTGGCATTAACAAAGTGCCGAATGAATATCATGTcaattagaccgagtcg                                                                                                                                                                                  >  2:43865/4‑236 (MQ=255)
                                                                    cAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGGGCGCAATGAATATCATGTCTATATGACCGATGCGCAg                                                                                                                                                                               <  1:371166/251‑1 (MQ=255)
                                                                             tCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGt                                                                                                                                                                       <  1:38035/250‑1 (MQ=255)
                                                                              cGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTAt                                                                                                                                                                     <  1:325787/251‑1 (MQ=255)
                                                                                           tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGCTTTGATCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAAACCATCCACCGTTTCCCGCCAATATCGGTGGAATTAACAAAGTGCGCAATGACTATCATGACTATATGACCGATGCGCCGGGCAAAGTATTGTTCGATTCGGc                                                                                                                                                        >  2:295256/1‑251 (MQ=255)
                                                                                           tttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCta                                                                                                                                                                                               >  1:373734/1‑212 (MQ=255)
                                                                                           tttGTCAAAGAAGTTAAACCGGGCGGGCGAAGAGCAACGGAGGGGACGTTGATCGACTCCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCta                                                                                                                                                                                               <  2:373734/212‑1 (MQ=255)
                                                                                             tGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCa                                                                                                                                                       <  1:156095/250‑1 (MQ=255)
                                                                                              gTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCaa                                                                                                                                                      >  1:444128/1‑250 (MQ=255)
                                                                                                           aaCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGAAAATAAAGCCGTTggg                                                                                                                                        >  2:55037/1‑250 (MQ=255)
                                                                                                              cGGGCGTGCGAAGAGCAACGGAGGGGCCGTTGATCGACAACGCAACGTTTCTGGCGGAGCTGGCGCGTCCCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATAGTTCGATTCGGCAAATAAAGCCGTTGGActg                                                                                                                                     <  2:286864/251‑3 (MQ=255)
                                                                                                                   gTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGatt                                                                                                                                 >  1:442787/1‑250 (MQ=255)
                                                                                                                   gTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAAAATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCACAGTATTGTTCGCTTCGGCAAATAAAGCCGTTGGACAGGattt                                                                                                                                >  1:51882/1‑250 (MQ=255)
                                                                                                                   gTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGATGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCCTTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGGATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGatt                                                                                                                                 >  1:394840/1‑250 (MQ=255)
                                                                                                                           agCAACGGAGGGGACGTTAATCGATACCGCAACGTTACTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCt                                                                                                                        <  1:244566/251‑1 (MQ=255)
                                                                                                                                       gACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCt                                                                                                            <  1:456900/251‑1 (MQ=255)
                                                                                                                                         cggtGATCGACACTGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCGCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGTCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATCAAGCCTTTGGACAGGATTATTCGCGCTGGAATGACGTCTgg                                                                                                          <  2:13648/248‑1 (MQ=255)
                                                                                                                                           ttGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGAAAGGCTTATTCGCGCTGGAATGACGTCTggat                                                                                                        >  2:315265/1‑249 (MQ=255)
                                                                                                                                             gATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACAGCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTa                                                                                                       <  1:213877/250‑1 (MQ=255)
                                                                                                                                                  acacCGCAACGTTGCTGACGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCCAAGCATGGGCGACTGGCACAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCAGCGGCCTTCACAAAGGGCGCAATGAATATCATGTCTATAAGACCGCTGCGCAGGGCCACGTATTGTTCGCTTAAGCAAATAAAGCCGTTGGACTGGTTTATTCGCGCTGGAATGACGTCTGACTAACATTg                                                                                                 >  2:277732/1‑251 (MQ=255)
                                                                                                                                                      cGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTg                                                                                             >  1:80322/1‑251 (MQ=255)
                                                                                                                                                                   gCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTgc                                                                                <  1:315265/251‑1 (MQ=255)
                                                                                                                                                                     ggAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTgcgc                                                                              >  1:451899/1‑251 (MQ=255)
                                                                                                                                                                                    cGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCACCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCaa                                                                >  1:337320/1‑250 (MQ=255)
                                                                                                                                                                                     aaTTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCACTGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCaaa                                                               >  2:294024/2‑250 (MQ=255)
                                                                                                                                                                                        ttGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATc                                                            <  2:77073/250‑1 (MQ=255)
                                                                                                                                                                                            cctcgtgggACCCAACGCATGGGCAACTGGCTCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGcc                                                       <  2:442787/245‑1 (MQ=255)
                                                                                                                                                                                              tctGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCg                                                      >  1:197549/1‑250 (MQ=255)
                                                                                                                                                                                               ctGGGGACCCAACGCCTGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGCCCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGATTGACGTCTGGCTAACGTTACGGGGTCAGTCTGGGGCGCGCAGCACGTTGCAAAATCCTGCCgct                                                    >  2:120840/1‑249 (MQ=255)
                                                                                                                                                                                                tGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATc                                                   >  1:250938/1‑251 (MQ=255)
                                                                                                                                                                                                 ggggACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAACGATGCGCAGGGCAAAGTATTGTACGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATcc                                                  >  2:266087/1‑251 (MQ=255)
                                                                                                                                                                                                  gggACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACAACCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCTTGCCGATccc                                                 <  2:446241/251‑1 (MQ=255)
                                                                                                                                                                                                   ggACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCg                                                >  1:174672/1‑251 (MQ=255)
                                                                                                                                                                                                       ccAACGCGTGGGCACCGGGAGCAGGCTTTTAATCAGCTTCAACATCCCCCTTTTCGCTCCAATATCGGTGGTATTAACAAAGTGAGCATTTAATACAATGTTTATATGACCGATGCGCGGGGCAAAGTATTGTTCGATTCTGCAAATAAAGCCGTTGGACAGTATTATTCTCGCTGGAATGACTTCTGGCTAACGTTGCTTGGTCAGTAGGGTGCGCGCAGCACGTTGCAAAATCCTGCAGTTCCCGAAAg                                            <  2:352143/251‑1 (MQ=255)
                                                                                                                                                                                                                gggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGAGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGat                                   >  1:469197/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGACATTAACAAAGTGAGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGATGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCGGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCAGAAAGTTCTGTGatg                                  >  1:420718/1‑250 (MQ=255)
                                                                                                                                                                                                                   cAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGatgta                                <  1:187187/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggggtccggtttttcatcatcTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATAAACAAAGTGCGCAATTAATATCATGTCTAAATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATGCGCGCTGGAATTCCGTCTGGCTAACGTGGCGTGGTCAGTTTGTTGCGCGCCGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGtt                            <  2:310801/236‑1 (MQ=255)
                                                                                                                                                                                                                        ggCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTg                           <  1:55037/251‑1 (MQ=255)
                                                                                                                                                                                                                           gcAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGCTGCGCAGGGAAAAGTATTGTTCGATTCGGCAAATAAAGCAGTTGGCCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGCAAGTTCTGTGATGTATGTTGCCg                        >  2:416133/1‑251 (MQ=255)
                                                                                                                                                                                                                               gCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCAcc                    >  1:114378/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aGCTACAACTTCCCCCGTTTCGCGCCAATATCGGGGGCATTAACAAAGTGCGCAATGAATATCATTTCTATATGACCGATGCGCAGGGCAAAGTATGGTCCGATTCTGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGAc          <  2:339236/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 caTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGCAAAGTATTGTTCGATCCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTAGGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTGCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTcgcg  <  2:451899/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCG  >  NC_000913/4636680‑4637171

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: