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→GGTflgKflagellar 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

GTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAC  >  NC_000913/1138287‑1138755
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
gTATGAACATGGAAAATCTGTTCTGAATAACTCAAGTCGGGGGGGTCGCGGCGGATAAGACTCTGTAATTGAAGGCTTATAAGGATCCTCCATGCCCAGCTTGATTAAAAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTggg                                                                                                                                                                                                                             <  2:313207/251‑1 (MQ=255)
                   gTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCt                                                                                                                                                                                                              <  1:519613/247‑1 (MQ=255)
                    gtcggAATAACCAAAGTCCGGCGGGTGTCTGCCGATAATCCTTTGTAATTGAAGGCTTTTAAGGAACCTCCATGTCCAGCTTGATTAAGAACGCCAGGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGt                                                                                                                                                                                                          <  1:471338/246‑1 (MQ=255)
                        gaaaaacccAAGTCCGGGGGGTCGCTGCGGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCAGGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGTTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTc                                                                                                                                                                                                       <  2:271793/241‑1 (MQ=255)
                        gAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTc                                                                                                                                                                                                       <  1:291908/249‑1 (MQ=255)
                           tAACTCAAGTCCGGCGGGTCGCTGCCGCTAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGAGCTGGCGGGTGGGGTGGGCATGGTGTCTACGTTTCTtgtg                                                                                                                                                                                                  >  1:353946/1‑251 (MQ=255)
                             aCTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCCATGGGGTCTACGTTTCTGGtgtg                                                                                                                                                                                                >  1:208345/1‑251 (MQ=255)
                                     ccGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTTGGGTGGCAATGGGGTTTACGTTTCTGGTGGGCAgagggg                                                                                                                                                                                        >  1:495497/1‑246 (MQ=255)
                                        gcggGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTa                                                                                                                                                                                     <  1:82851/251‑1 (MQ=255)
                                                         atactcctaaatggAAGGCTTTTAGGGAACCGCCACGTCCAGCTTGATTAAGAACGCCAGGAGCGGACTGACCGCGGCCGGGGCGGCGTAAAATACGGCAAGTATTAATATCTCCAGCTAGAACGTGTCTGTATAGACCCGCCAAAGCACTATTATGGCGCAGGCCAAGAGCTCGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTacca                                                                                                                                                                    <  2:133728/238‑1 (MQ=255)
                                                          taCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCCATGGGGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTacca                                                                                                                                                                    >  1:439524/1‑250 (MQ=255)
                                                           aCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTaccaac                                                                                                                                                                  <  1:180865/251‑1 (MQ=255)
                                                                gTAAGTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGtt                                                                                                                                                            <  1:626952/251‑1 (MQ=255)
                                                                  aaTTGAAGGCTTATAAGGACCCCCCATGCCCAGCTTGATTAATACCGCCATGAGCGGACTGAACCCGGCCCAGGCGGCTTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGTATATCCCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTTAGTATGATGCGTTTATTACCAACCAGTTAc                                                                                                                                                          <  1:160878/251‑1 (MQ=255)
                                                                           cTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCCATGGTGTCTTCGTTTCTTGTGTGCAGCGTGGGTATGATGCGTTTTTTACCAACCAGTTACGTTCGgcgc                                                                                                                                                 >  1:321349/1‑251 (MQ=255)
                                                                             tataAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGCAAGGAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTTGGTGGGCAATGGGGTCTACGGTTCTGGGGTGCAGGGGGGGTATGAGGCGTTTTTTTCCCACCAGTGTCGTGTGgcgc                                                                                                                                                 >  1:141790/1‑248 (MQ=255)
                                                                                aaGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTGGGGAATGGTGTCTACGTTTTTGGTGTGCAGCGTGAGTGTGGTGCGGTTATTACAAACCCATTACCTGGGGCGCATACg                                                                                                                                            >  2:602431/1‑251 (MQ=255)
                                                                                             gTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCCGTTACGTGCGGCGCAGACGCAAAGTAGCGGTc                                                                                                                               >  1:28125/1‑251 (MQ=255)
                                                                                                           ataaCGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCCCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTAt                                                                                                                 <  1:189999/251‑1 (MQ=255)
                                                                                                             aaCGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCCCTGGCGGATGGGTTGGGAATGGGGTCTACGTTTTTGGTGGGCAGCGTGAGGATGAGGCGTTTATTAACCACCAGATAAGTGCGGCGCAGACGCAAAGGAGCggcacgactgcacgcaatt                                                                                                                >  1:540782/1‑233 (MQ=255)
                                                                                                                     gAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCACCAGCTATAACGTTGCCGGATATACCCGCCAAACAACTATTACGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGATGGGTTGGCAATGGAGGCAACGATTCTGGCGTGCAACGAGAGTATGATGAGCTTATTACCAACCAGTTTCGGACGGCGCCGACGCCAAGTAGCGGTCTGACAGCCCGCtctgtgccgaggg                                                                                                       >  2:347438/1‑239 (MQ=255)
                                                                                                                               aaCGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTTGGTTGGCCATGGGGGCTACGTTTCTGGTGTGCAGCGTGAGGCTGATGGGTTTATTACCCACCAGTTCCTTGCGGCGCAGACGCAAAAGAGCGTTCTTACTGCCCGCTATTAAAAGATGTCGAAAAt                                                                                                >  2:531339/1‑248 (MQ=255)
                                                                                                                                 cgcgGCCCAGGCGGCGTTAAATCCGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACa                                                                                           <  1:571744/251‑1 (MQ=255)
                                                                                                                                 cgcgGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGGGGCTAAGTTTCTGGTGTGCAGCGGGAGGATGGTGCGTTTAATACCAAACAGTTACGTGGGGCGCAGACGCAAAGTAGCGGGCTGAATGCCCGCCATGAGCAGGAGTCGACAATCGAc                                                                                            >  1:79148/1‑250 (MQ=255)
                                                                                                                                      ccagggcggcGTTAAAACCGTCAAGTAATAATAGACCCAGCTAAAACGTGGCTGGATATACCCGCCAAACCACTAATATGGCGCAGGCGAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCGACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATg                                                                                      <  2:28125/247‑1 (MQ=255)
                                                                                                                                               cgTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTTCCGGATATACCCGACAAAACACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGGGAGTATAGTGCGTTTATTACCAACCCGTTTAGTGCGGGGCAGACGCAAAGTAGCGGTCTGGCTGTCCGCTATGAGCAGCTGTCGAAAATCGACCATAAGCTCTccacc                                                                             >  1:590104/1‑251 (MQ=255)
                                                                                                                                                gTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAAACGACAATATGGTCTccacac                                                                            >  1:88621/1‑249 (MQ=255)
                                                                                                                                                  taaatctggctaGTAAAAATATCCCCAGCTATATCGGTGGCGGTAATCCCCGCCAAACCTCTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAg                                                                           <  2:226917/239‑1 (MQ=255)
                                                                                                                                                        cGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCCAAGTAGCGGTCTGACTGCCCGCTATGAGccg                                                                                                           >  2:285587/1‑210 (MQ=255)
                                                                                                                                                        cGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAg                                                                                                           <  1:285587/212‑1 (MQ=255)
                                                                                                                                                         ggCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCCGTACCTCtt                                                                   >  1:438303/1‑251 (MQ=255)
                                                                                                                                                                     tatCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGGGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGCCGCCAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGGCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAacac                                                       >  1:621818/1‑251 (MQ=255)
                                                                                                                                                                      atCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTAt                                                                                                                 <  1:568868/192‑1 (MQ=255)
                                                                                                                                                                      atCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAacca                                                       >  1:536420/1‑248 (MQ=255)
                                                                                                                                                                      atCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCTACGTTTCTGGTGGGCAGCGGGAGGATGATGCGTTTATTAACAACCAGGTTAGGTCGGCGCAGACGCAAAGGAGCGGTCTTACTGCCCGCTAt                                                                                                                 >  2:568868/1‑192 (MQ=255)
                                                                                                                                                                                  aaCGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGGGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGGCAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGAttt                                          >  1:195518/1‑251 (MQ=255)
                                                                                                                                                                                               atccccGACAAACCATTACTTTGGCGCAGGCCAAGAGCACGGTGGGCGCTGGCGGTGTGGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCACGCTATGAGCCGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAAAACAGATGCAGTATTTCTTCTCCAGcat                              <  2:353946/247‑3 (MQ=255)
                                                                                                                                                                                                       ccaaCCACTATTATGGCGCAGGCCAATAGCACGTTGGGTGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTTTGGTGTGCAGCGTGAGTATGATTCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGCCTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCCCCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCt                     <  2:13508/249‑1 (MQ=255)
                                                                                                                                                                                                       caacccACTATTATGGCGCAGGCCAATAGCACGTGGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCt                     <  2:438303/247‑1 (MQ=255)
                                                                                                                                                                                                       cAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGAGGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCt                     <  1:50964/251‑1 (MQ=255)
                                                                                                                                                                                                          aCCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGGGTCTACGTTTCTGGTGTGCAGCGGGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGGACCTCTTCGCTGGCCACACAGATGCAGGATTTCTTCACCCGCCTGCCAACGCTggg                  >  1:603310/1‑250 (MQ=255)
                                                                                                                                                                                                                attatGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGGCATGGTGTCCACGTTTCTTGTGTGCAGCGTGAGGATGATGCGTTTATTACCAACCAGTTTAGTGTGGCCCAGACGCAAAATAGCGGTCTGACTGCCCGATATATGCAGATGGCGAAAATCGACAATATGCTCTTCACAAGTACCTCTTCGCTGGCACCACAGATGCAGGGTTTCTTCACCCGCCTGCAAACGCTGGTGAGTa             >  1:587391/1‑250 (MQ=255)
                                                                                                                                                                                                                attatGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTTTACGATTCTGGGGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACTCAAAGTAGAGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTaa            <  2:185552/251‑1 (MQ=255)
                                                                                                                                                                                                                   atGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAAcg          <  1:628151/250‑1 (MQ=255)
                                                                                                                                                                                                                         cAGGCCAATAGCACGTTGGGCGCTGGCGGGTGGGGTGGCAATGGTGGCCACCTTTCTGTGGGGCCGGGTGAGGATGATTCGTTTTTTACCAACCGGTTACGTTCGGGGCCGACGCACAATAGCGGGCTTACTGCCCACTATGAACCGATATCGAAAAACGACAACATGCGCTCCACCAGTACCTCTTCACTGGCAAACCATCTACAGGGTTTCTTCACCCGCCTGCAAACGccagggaataaagcaggagg   >  2:525660/1‑232 (MQ=255)
                                                                                                                                                                                                                          aGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAc  <  1:478712/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGccc                                                                                                                      <  1:171250/122‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
GTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAC  >  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: