Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,533 C→T G52G (GGC→GGT)  flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,5330CT95.5% 121.7 / NA 47G52G (GGC→GGT) flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/23);  minor base G (2/0);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAC  >  NC_000913/1138287‑1138755
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
gTATGAACATGGAAAATCTGTTCTGAATAACTCAAGTCGGGGGGGTCGCGGCGGATAAGACTCTGTAATTGAAGGCTTATAAGGATCCTCCATGCCCAGCTTGATTAAAAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTg‑gg                                                                                                                                                                                                                             <  2:313207/251‑1 (MQ=255)
                   gTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCt                                                                                                                                                                                                              <  1:519613/247‑1 (MQ=255)
                    gtcggAATAACCAAAGTCCGGCGGGTGTCTGCCGATAATCCTTTGTAATTGAAGGCTTTTAAGGAACCTCCATGTCCAGCTTGATTAAGAACGCCAGGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGt                                                                                                                                                                                                          <  1:471338/246‑1 (MQ=255)
                        gaaaaacccAAGTCCGGGGGGTCGCTGCGGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCAGGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGTTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTc                                                                                                                                                                                                       <  2:271793/241‑1 (MQ=255)
                        gAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTc                                                                                                                                                                                                       <  1:291908/249‑1 (MQ=255)
                           tAACTCAAGTCCGGCGGGTCGCTGCCGCTAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGAGCTGGCGGGTG‑GGGTGGGCATGGTGTCTACGTTTCTtgtg                                                                                                                                                                                                  >  1:353946/1‑251 (MQ=255)
                             aCTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCCATGGGGTCTACGTTTCTGGtgtg                                                                                                                                                                                                >  1:208345/1‑251 (MQ=255)
                                     ccGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTT‑GGGTGGCAATGGGGTTTACGTTTCTGGTGGGCAgagggg                                                                                                                                                                                        >  1:495497/1‑246 (MQ=255)
                                        gcggGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTa                                                                                                                                                                                     <  1:82851/251‑1 (MQ=255)
                                                         atactcctaaatggAAGGCTTTTAGGGAACCGCCACGTCCAGCTTGATTAAGAACGCCAGGAGCGGACTGACCGCGGCCGGGGCGGCGTAAAATACGGCAAGTATTAATATCTCCAGCTAGAACGTGTCTGTATAGACCCGCCAAAGCACTATTATGGCGCAGGCCAAGAGCTCGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTacca                                                                                                                                                                    <  2:133728/238‑1 (MQ=255)
                                                          taCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCCATGGGGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTacca                                                                                                                                                                    >  1:439524/1‑250 (MQ=255)
                                                           aCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTaccaac                                                                                                                                                                  <  1:180865/251‑1 (MQ=255)
                                                                gTAAGTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑tt                                                                                                                                                            <  1:626952/251‑1 (MQ=255)
                                                                  aaTTGAAGGCTTATAAGGACCCCCCATGCCCAGCTTGATTAATACCGCCATGAGCGGACTGAACCCGGCCCAGGCGGCTTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGTATATCCCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTTAGTATGATGCGTTTATTACCAACCAG‑TTAc                                                                                                                                                          <  1:160878/251‑1 (MQ=255)
                                                                           cTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCCATGGTGTCTTCGTTTCTTGTGTGCAGCGTGGGTATGATGCGTTTTTTACCAACCAG‑TTACGTTCGgcgc                                                                                                                                                 >  1:321349/1‑251 (MQ=255)
                                                                             tataAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACG‑CAAGGAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTT‑GGTGGGCAATGGGGTCTACGGTTCTGGGGTGCAGGGGGGGTATGAGGCGTTTTTTTCCCACCAG‑TGTCGTGTGgcgc                                                                                                                                                 >  1:141790/1‑248 (MQ=255)
                                                                                aaGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTGGGGAATGGTGTCTACGTTTTTGGTGTGCAGCGTGAGTGTGGTGCGGTTATTACAAACCCA‑TTACCTGGGGCGCATACg                                                                                                                                            >  2:602431/1‑251 (MQ=255)
                                                                                             gTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCCG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTc                                                                                                                               >  1:28125/1‑251 (MQ=255)
                                                                                                           ataaCGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTAt                                                                                                                 <  1:189999/251‑1 (MQ=255)
                                                                                                             aaCGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCCCTGGCGGATG‑GGTTGGGAATGGGGTCTACGTTTTTGGTGGGCAGCGTGAGGATGAGGCGTTTATTAACCACCAG‑ATAAGTGCGGCGCAGACGCAAAGGAGCggcacgactgcacgcaatt                                                                                                                >  1:540782/1‑233 (MQ=255)
                                                                                                                     gAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCACCAGCTATAACGTTGCCGGATATACCCGCCAAACAACTATTACGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGATG‑GGTTGGCAATGGAGGCAACGATTCTGGCGTGCAACGAGAGTATGATGAGCTTATTACCAACCAG‑TTTCGGACGGCGCCGACGCCAAGTAGCGGTCTGACAGCCCGCtctgtgccgaggg                                                                                                       >  2:347438/1‑239 (MQ=255)
                                                                                                                               aaCGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTT‑GGTTGGCCATGGGGGCTACGTTTCTGGTGTGCAGCGTGAGGCTGATGGGTTTATTACCCACCAG‑TTCCTTGCGGCGCAGACGCAAAAGAGCGTTCTTACTGCCCGCTATTAAAAGATGTCGAAAAt                                                                                                >  2:531339/1‑248 (MQ=255)
                                                                                                                                 cgcgGCCCAGGCGGCGTTAAATCCGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACa                                                                                           <  1:571744/251‑1 (MQ=255)
                                                                                                                                 cgcgGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATGGGGGCTAAGTTTCTGGTGTGCAGCGGGAGGATGGTGCGTTTAATACCAAACAG‑TTACGTGGGGCGCAGACGCAAAGTAGCGGGCTGAATGCCCGCCATGAGCAGGAGTCGACAATCGAc                                                                                            >  1:79148/1‑250 (MQ=255)
                                                                                                                                      ccagggcggcGTTAAAACCGTCAAGTAATAATAGACCCAGCTAAAACGTGGCTGGATATACCCGCCAAACCACTAATATGGCGCAGGCGAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCGACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATg                                                                                      <  2:28125/247‑1 (MQ=255)
                                                                                                                                               cgTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTTCCGGATATACCCGACAAAACACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGGGAGTATAGTGCGTTTATTACCAACCCG‑TTTAGTGCGGGGCAGACGCAAAGTAGCGGTCTGGCTGTCCGCTATGAGCAGCTGTCGAAAATCGACCATAAGCTCTccacc                                                                             >  1:590104/1‑251 (MQ=255)
                                                                                                                                                gTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGGGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAAACGACAATATGGTCTccacac                                                                            >  1:88621/1‑249 (MQ=255)
                                                                                                                                                  taaatctggctaGTAAAAATATCCCCAGCTATATCGGTGGCGGTAATCCCCGCCAAACCTCTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAg                                                                           <  2:226917/239‑1 (MQ=255)
                                                                                                                                                        cGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCCAAGTAGCGGTCTGACTGCCCGCTATGAGccg                                                                                                           >  2:285587/1‑210 (MQ=255)
                                                                                                                                                        cGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAg                                                                                                           <  1:285587/212‑1 (MQ=255)
                                                                                                                                                         ggCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCCGTACCTCtt                                                                   >  1:438303/1‑251 (MQ=255)
                                                                                                                                                                     tatCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGGGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGCCGCCAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGGCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAacac                                                       >  1:621818/1‑251 (MQ=255)
                                                                                                                                                                      atCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTAt                                                                                                                 <  1:568868/192‑1 (MQ=255)
                                                                                                                                                                      atCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAacca                                                       >  1:536420/1‑248 (MQ=255)
                                                                                                                                                                      atCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGGGTCTACGTTTCTGGTGGGCAGCGGGAGGATGATGCGTTTATTAACAACCAGGTTAGGT‑CGGCGCAGACGCAAAGGAGCGGTCTTACTGCCCGCTAt                                                                                                                 >  2:568868/1‑192 (MQ=255)
                                                                                                                                                                                  aaCGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGGGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGGCAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGAttt                                          >  1:195518/1‑251 (MQ=255)
                                                                                                                                                                                               atccccGACAAACCATTACTTTGGCGCAGGCCAAGAGCACGGTGGGCGCTGGCGG‑TGTGGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCACGCTATGAGCCGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAAAACAGATGCAGTATTTCTTCTCCAGcat                              <  2:353946/247‑3 (MQ=255)
                                                                                                                                                                                                       ccaaCCACTATTATGGCGCAGGCCAATAGCACGTTGGGTGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTTTGGTGTGCAGCGTGAGTATGATTCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGCCTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCCCCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCt                     <  2:13508/249‑1 (MQ=255)
                                                                                                                                                                                                       caacccACTATTATGGCGCAGGCCAATAGCACGTGGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCt                     <  2:438303/247‑1 (MQ=255)
                                                                                                                                                                                                       cAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGAGGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCt                     <  1:50964/251‑1 (MQ=255)
                                                                                                                                                                                                          aCCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGGTGGCAATGGGGTCTACGTTTCTGGTGTGCAGCGGGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGGACCTCTTCGCTGGCCACACAGATGCAGGATTTCTTCACCCGCCTGCCAACGCTggg                  >  1:603310/1‑250 (MQ=255)
                                                                                                                                                                                                                attatGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGGCATGGTGTCCACGTTTCTTGTGTGCAGCGTGAGGATGATGCGTTTATTACCAACCAG‑TTTAGTGTGGCCCAGACGCAAAATAGCGGTCTGACTGCCCGATATATGCAGATGGCGAAAATCGACAATATGCTCTTCACAAGTACCTCTTCGCTGGCACCACAGATGCAGGGTTTCTTCACCCGCCTGCAAACGCTGGTGAGTa             >  1:587391/1‑250 (MQ=255)
                                                                                                                                                                                                                attatGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGGGTTTACGATTCTGGGGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACTCAAAGTAGAGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTaa            <  2:185552/251‑1 (MQ=255)
                                                                                                                                                                                                                   atGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAAcg          <  1:628151/250‑1 (MQ=255)
                                                                                                                                                                                                                         cAGGCCAATAGCACGTTGGGCGCTGGCGGGTG‑GGGTGGCAATGGTGGCCACCTTTCTGTGGGGCCGGGTGAGGATGATTCGTTTTTTACCAACCGG‑TTACGTTCGGGGCCGACGCACAATAGCGGGCTTACTGCCCACTATGAACCGATATCGAAAAACGACAACATGCGCTCCACCAGTACCTCTTCACTGGCAAACCATCTACAGGGTTTCTTCACCCGCCTGCAAACGccagggaataaagcaggagg   >  2:525660/1‑232 (MQ=255)
                                                                                                                                                                                                                          aGGCCAATAGCACGTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAc  <  1:478712/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gTTGGGCGCTGGCGGTTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGccc                                                                                                                      <  1:171250/122‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
GTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTG‑GGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAG‑TTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAC  >  NC_000913/1138287‑1138755

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