Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,016,354 C→A 16.1% A330E (GCA→GAA)  fliF → flagellar basal body M‑ring protein FliF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077792,016,3540CA16.1% 193.0 / 4.9 69A330E (GCA→GAA) fliFflagellar basal body M‑ring protein FliF
Reads supporting (aligned to +/- strand):  ref base C (44/12);  new base A (2/10);  total (46/23)
Fisher's exact test for biased strand distribution p-value = 9.65e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTG  >  NC_007779/2016117‑2016586
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
cGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGc                                                                                                                                                                                                                                          <  1:1394911/238‑1 (MQ=255)
   cTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCATCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCaaa                                                                                                                                                                                                                                       <  2:1848092/238‑1 (MQ=255)
   cTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAAt                                                                                                                                                                                                                                      <  1:1096203/239‑1 (MQ=255)
         ttgGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCACTCAAGAACAAACGGAAGCACAGTATCGCCCTAACGGTGCTGAATCGCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGCGCGAGCACCGTGGTTCCGGTTATCCGGGCGGCGCACCGGCGGCGTTGTCGAATCAACCGGCACCTGCGACTCCAGCGTCACGCAGCACGCACCCGGCACTTCaaacc                                                                                                                                                                                                                                >  2:160706/1‑237 (MQ=255)
           ggTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATa                                                                                                                                                                                                                               >  1:1702085/1‑238 (MQ=255)
                           aCGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCg                                                                                                                                                                                                              >  1:967278/1‑239 (MQ=255)
                             gcgcAGGTTATGGCGCAGCTGGACATCGCGAGTAAAGAACAAACGGAAGAACAGTAACGCCCTAACGGTGACGAATCACATGCGGCGCTACGTTCACGCCAGCTTAATGAGAGCGAGCACAGCGGTTCCGGTTTTCCGAGCTGCGTACCGGGGGCGTTGTCGAGTCAACAGGCACCTGCGACTAACGCGCCAACCAGCACGCCGCCGGAGTATCAAAATAAACGCacgcaacagtagc                                                                                                                                                                                                             >  2:33886/4‑234 (MQ=255)
                             gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCaa                                                                                                                                                                                                                                        <  1:1653058/211‑1 (MQ=255)
                             gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCaa                                                                                                                                                                                                                                        >  2:1653058/1‑211 (MQ=255)
                             gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGa                                                                                                                                                                                                             >  1:716375/1‑238 (MQ=255)
                             gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGa                                                                                                                                                                                                             >  1:1690009/1‑238 (MQ=255)
                             gCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGa                                                                                                                                                                                                             >  1:1529269/1‑238 (MQ=255)
                               caggtttACGGCGCCGCTGAACTGCACCAGTAAATAACAAACGGAAGAACGGTACCGCCCTAGCGGTGATGAATCTCATGCGGCGCTTGTGTACCGCCAGCTTAAGGGGGGCGAGCAAAGCGTTTCCGGTTATCCGGGCGGCGTACCGGGGGGTTTGCCAAACCACCCGCCCCTTGAGAATACCGCGCAATCCACCCCCCCCCCGGCAAATAAAAAAAACCGCCAGCAGCAGGCGAGca                                                                                                                                                                                                          <  2:700575/234‑1 (MQ=255)
                                  ggTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGcacc                                                                                                                                                                                                        >  1:1173820/1‑238 (MQ=255)
                                       cGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGc                                                                                                                                                                                                  >  1:1991903/1‑239 (MQ=255)
                                            cAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAg                                                                                                                                                                                             >  1:8601/1‑239 (MQ=255)
                                                tGGATTTCGCAAGTAAAGAACAAACGGAAAAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGTTTAATGAGAGCGAGCAAGGCGGTTCCGTTTACCCGGGCGGCGTCCCGGGGGCGTTGTCGAATCACCCGGCCCCTGCGAATAACGCGCAAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTggg                                                                                                                                                                                         <  1:1517463/239‑1 (MQ=255)
                                                tGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTACCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTCCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCCCCTGCGAATAACGCGCCAATCAGCACGCCTCCGGAAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTggg                                                                                                                                                                                         <  1:1942548/239‑1 (MQ=255)
                                                 ggACTTCGCAAGTAAAGAACAAACGGAAGAACAGTACCGCCCTAACGGTGATGAATCTCAGGCGGCGCTTCGTTCACGCCGGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTACCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGAAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTggg                                                                                                                                                                                         <  1:2012127/238‑1 (MQ=255)
                                                         ccAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGc                                                                                                                                                                                                                                          >  1:1926898/1‑181 (MQ=255)
                                                              aaaGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGcaca                                                                                                                                                                             <  1:163268/237‑1 (MQ=255)
                                                                       aCGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGa                                                                                                                                                                  >  1:195623/1‑239 (MQ=255)
                                                                       aCGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAAGGAGAGCGAGCAAAGCGGTTCCGGTTACCCGGGCGGCTTACCGGGGGGGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGAAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGa                                                                                                                                                                  <  1:1812246/239‑1 (MQ=255)
                                                                        cGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCAGGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTCCCGGGGGCGTTGTCGAATCAACCGGCCCCTGCGAATAACGCGCCAATCAGCACGCCTCCGGAAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGaa                                                                                                                                                                 <  1:1341174/239‑1 (MQ=255)
                                                                            agaaCAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCa                                                                                                                                                             >  2:1838466/1‑239 (MQ=255)
                                                                                   tATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGCAGCACACAGCGGAATGAAACCAGTAACt                                                                                                                                                       >  1:494384/1‑238 (MQ=255)
                                                                                             aCGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCCCCCAGCGGAATGAAACCAGTAAATACGAAGTCGa                                                                                                                                             >  2:710516/1‑238 (MQ=255)
                                                                                             aCGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGa                                                                                                                                             >  2:1175895/1‑238 (MQ=255)
                                                                                                   atgaATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGcacca                                                                                                                                                                                                       >  1:1735099/1‑174 (MQ=255)
                                                                                                   atgaATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATa                                                                                                                                                                                              >  1:1318632/1‑183 (MQ=255)
                                                                                                   atgaATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCAcc                                                                                                                                      >  1:1206501/1‑239 (MQ=255)
                                                                                                   atgaATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCAcc                                                                                                                                      >  1:1090476/1‑239 (MQ=255)
                                                                                                   atgaATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCAcc                                                                                                                                      >  1:344483/1‑239 (MQ=255)
                                                                                                   atgaATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCAcc                                                                                                                                      >  1:659549/1‑239 (MQ=255)
                                                                                                         ctcATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCg                                                                                                                                 >  1:295309/1‑238 (MQ=255)
                                                                                                           cATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATa                                                                                                                                                                                              >  1:869978/1‑175 (MQ=255)
                                                                                                             tGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCAt                                                                                                                             >  1:1754535/1‑238 (MQ=255)
                                                                                                                 gcgcTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAg                                                                                                                                                                          >  1:34019/1‑189 (MQ=255)
                                                                                                                   gctctgTTCACGCCGGCTAGAGGGGGGGGAGCAAGGCGGTCCCGGTTATCCGGGCGGCGTCCCGGGGGGGTGGTCGAACAACCGGCCCCCTGAAAATACCGCGCCAACCACCCCCCCGCCGGAAAATAAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCaaaa                                                                                                                      <  2:295309/234‑1 (MQ=255)
                                                                                                                       cGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCCGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCGTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATACTGGGCCGCGTAGCACAAATCGGAATGCAACCACAAACTACAACGTCGATCGCACCATCCGACACACGAAAATGa                                                                                                                   >  2:1982385/1‑238 (MQ=255)
                                                                                                                                  gttgaaTGGGAGCGAGCAAGGGGGTCCCGGTTAGCCGGGCGGCGTACCGGGGGCGTTGTCGAACCACCCGCCCCCTCCGAAAAACGCGCAAATAACCCCCCCCCCGGAAAATAAAAATAACCGCCAGCAGCAGGCGAg                                                                                                                                                                                                            <  2:1861183/134‑1 (MQ=255)
                                                                                                                                   cTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTGCCGGGGGCGTTGTAGAATCAACGGGCAACGGCGAGTAACGCGCCAATCAGCACGCCGCCGGCAAATAAAAATAACCGCAAGCAGCAGGCGAGCACCACCAGCAGTAGGGGGCCGTGTAGCAAACAGCGGAATGAAAACAGTACCTACGCAGCAGAGCGCAGCATCCGTCACACCACAACGAACGTGCGCGAtgt                                                                                                      >  2:1439066/1‑239 (MQ=255)
                                                                                                                                   cTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGAtgt                                                                                                      >  1:92804/1‑239 (MQ=255)
                                                                                                                                          gaggagagcAAATCGTTTCCGGTCAACGGGTGACGGTACAGGGGGCGTATACGAACCACCCGCACCCTGCAATACCGGACCAATCCGCGCCACCTCCAGAAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACg                                                                                                <  2:8601/233‑1 (MQ=255)
                                                                                                                                            agcgagcAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGcacc                                                                                                                                                                                                        >  2:1434304/1‑132 (MQ=255)
                                                                                                                                            agcgagcAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGcacc                                                                                                                                                                                                        <  1:1434304/132‑1 (MQ=255)
                                                                                                                                            agcgagcAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGc                                                                                                                                                                                                  >  1:1861183/1‑138 (MQ=255)
                                                                                                                                            agcgagcAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGt                                                                                                                                                                                            >  1:10620/1‑144 (MQ=255)
                                                                                                                                            agcgagcAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACg                                                                                                >  1:1641796/1‑236 (MQ=255)
                                                                                                                                                agcAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTgg                                                                                                             >  1:1887252/1‑219 (MQ=255)
                                                                                                                                                           tCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATAc                                                                                                                           >  1:676087/1‑194 (MQ=255)
                                                                                                                                                                ttAGCCGGGCGGCGCCCCGGGGGCGTGGGCAAAACAACCGGCCCCTGCGAATACCGCGCCAATCAGCACGCCTCCGGAAAAGCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCgtggtg                                                                          <  1:105509/238‑1 (MQ=255)
                                                                                                                                                                taaTCGGGGCGGCGTACCGGGGGCGTTGCCAAAACACCCGGCACCTGCAAATAACCCGCCAATCACCCCCCCCCCGGCAAATCAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCgtggtg                                                                          <  1:1547214/236‑1 (MQ=255)
                                                                                                                                                                    ccGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAAtt                                                                      >  1:1077722/1‑238 (MQ=255)
                                                                                                                                                                        gcggcgTCCGGGGGGCGTTGTCGAAACACCCGGCCCCTGCGAATAACGCGCAAACCACCACGCCTCCGGCAAACCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACaa                                                                  <  1:710516/238‑1 (MQ=255)
                                                                                                                                                                            cgTCCCGGGGGCGTTGGCGAATCACCCGGCACCTGCGAATAACGCGCCAATCAGCACCCCCCGGGAAAACAAAAAAAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAg                                                                                                                                                                          <  2:1281549/130‑1 (MQ=255)
                                                                                                                                                                                       gCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGaaa                                                                                                                                                                >  1:1281549/1‑130 (MQ=255)
                                                                                                                                                                                            gTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCaa                                                                                                  >  2:1108314/1‑187 (MQ=255)
                                                                                                                                                                                            gTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGAAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCaa                                                                                                  <  1:1108314/187‑1 (MQ=255)
                                                                                                                                                                                                  tCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGAAAATCAAACTAACCGACAGCAGCAGGCGAGCACCACCACCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTCCGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTAGCGGTCGTAGTGACTTACAGAACCGTCACAGATGGCAAACCGtcg                                         >  2:1248204/1‑236 (MQ=255)
                                                                                                                                                                                                     aCCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGcc                                       >  1:211686/1‑237 (MQ=255)
                                                                                                                                                                                                                gCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTCGCCTCgcagcacgga                           >  2:777653/1‑235 (MQ=255)
                                                                                                                                                                                                                        cgcgCCAATCAGCACGCCTCCGGCAAATCAGAATAACCGCCAGCAGCAGGAGAGCACCACCAGCAACAGTGGGCCGCGTAGCACACAGCGGAATGAAACGAGTAACTACGAAGTCTACCGCACCATTCGTCATACCAAAATCAACGTGGGCTATGTGCACAGTCTGTCAGTCGCGATCGCGGAGAATTACGAAACCGTGTCAGATAGCAAACCGTGGCCTCCCACCAACGACCAGATgg                   >  2:1500904/1‑238 (MQ=255)
                                                                                                                                                                                                                        cgcgCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGa                   >  1:1488570/1‑239 (MQ=255)
                                                                                                                                                                                                                         gcgcCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGaa                  <  1:1838466/239‑1 (MQ=255)
                                                                                                                                                                                                                              aaccAGCACGCCTCCGGCAAACAAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCaaa             <  1:1796976/236‑1 (MQ=255)
                                                                                                                                                                                                                                  ggctcccccTCCGATAAATAAAAAGAGCCGTCAGCAGCAGGCGATCACCCCCAGCAATAGTGGGCCGTGTAGCACACAGCGTTATGTAACCAGTAACTACGAAGTCGATCGCTCCATTCGTCATTCCAAAATGAACGTGGGCTATGTGCAACGTCTGTCATTCGCGGTCGTGTTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAACAACGAACAGATGAAGCAAATTg          <  2:2019865/232‑1 (MQ=255)
                                                                                                                                                                                                                                          cTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCa                                                        >  1:1567541/1‑185 (MQ=255)
                                                                                                                                                                                                                                          cTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTg  >  2:964668/1‑239 (MQ=255)
                                                                                                                                                                                                                                             cGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCa                                                        >  2:2008399/1‑182 (MQ=255)
                                                                                                                                                                                                                                             cGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCa                                                        <  1:2008399/182‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
CGCCTATTGTTGGTAACGGTAATATTCACGCCCAGGTTACGGCGCAGCTGGACTTCGCCAGTAAAGAACAAACGGAAGAACAGTATCGCCCTAACGGTGATGAATCTCATGCGGCGCTTCGTTCACGCCAGCTTAATGAGAGCGAGCAAAGCGGTTCCGGTTATCCGGGCGGCGTACCGGGGGCGTTGTCGAATCAACCGGCACCTGCGAATAACGCGCCAATCAGCACGCCTCCGGCAAATCAAAATAACCGCCAGCAGCAGGCGAGCACCACCAGCAATAGTGGGCCGCGTAGCACACAGCGGAATGAAACCAGTAACTACGAAGTCGATCGCACCATTCGTCATACCAAAATGAACGTGGGCGATGTGCAACGTCTGTCAGTCGCGGTCGTGGTGAATTACAAAACCTTGCCAGATGGCAAACCGTTGCCTCTCAGCAACGAACAGATGAAGCAAATTGAAGATCTG  >  NC_007779/2016117‑2016586

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: