Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,507,805 T→C E83G (GAA→GGA)  fecI ← RNA polymerase sigma factor FecI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,507,8050TC100.0% 119.1 / NA 32E83G (GAA→GGA) fecIRNA polymerase sigma factor FecI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

CGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCA  >  NZ_CP009273/4507686‑4507930
                                                                                                                       |                                                                                                                             
cGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGaaa                                                                                                            <  1:379769/139‑1 (MQ=255)
 gTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAAc                                                                                                           <  2:292130/139‑1 (MQ=255)
          gTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTccac                                                                                                   >  2:345787/1‑138 (MQ=255)
                cATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcg                                                                                            <  2:634991/139‑1 (MQ=255)
                cATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcg                                                                                            <  1:345787/139‑1 (MQ=255)
                cATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCCCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATCCGCTTTTCCCAGCGCTTTTCGGCGAAACAGGTCCACCATCAcg                                                                                            <  1:550304/139‑1 (MQ=255)
                                 gTAGGGTTTCGAGTTGGCTTTCGCGTTCCCCAGGTGAAGGCGTTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGc                                                                            <  1:479794/138‑1 (MQ=255)
                                 gTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTg                                                                             >  1:519413/1‑137 (MQ=255)
                                       tctcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGagg                                                                       >  1:257896/1‑137 (MQ=255)
                                         tcGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGaggaag                                                                    >  2:281829/1‑138 (MQ=255)
                                                       gcgTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCgg                                                     <  2:257896/139‑1 (MQ=255)
                                                                    gAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa                                                                                  >  1:407240/1‑97 (MQ=255)
                                                                    gAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa                                                                                  <  2:407240/97‑1 (MQ=255)
                                                                       ggCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGttt                                     >  1:504417/1‑139 (MQ=255)
                                                                       ggCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGttt                                     >  2:417459/1‑139 (MQ=255)
                                                                           cTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCt                                 <  2:504417/139‑1 (MQ=255)
                                                                            tCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa                                                                                  <  1:448711/89‑1 (MQ=255)
                                                                            tCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa                                                                                  >  2:448711/1‑89 (MQ=255)
                                                                            tCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGTAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCt                                 >  1:203629/1‑138 (MQ=255)
                                                                                   tCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTa                         >  1:477531/1‑139 (MQ=255)
                                                                                        cATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCg                     >  1:297326/1‑138 (MQ=255)
                                                                                        cATAAGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCg                     >  1:289024/1‑138 (MQ=255)
                                                                                            aGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCgg                                                                                                                 >  1:577707/1‑42 (MQ=255)
                                                                                            aGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCgg                                                                                                                 <  1:553244/42‑1 (MQ=255)
                                                                                            aGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCgg                                                                                                                 >  2:553244/1‑42 (MQ=255)
                                                                                            aGCGCCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCgg                                                                                                                 <  2:577707/42‑1 (MQ=255)
                                                                                              cgcCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCaaaaaa              <  2:303763/139‑1 (MQ=255)
                                                                                              cgcCAGCATCTCCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCaaaaaa              <  1:476050/139‑1 (MQ=255)
                                                                                                        tcCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCg                                   >  2:414925/1‑108 (MQ=255)
                                                                                                        tcCAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCg                                   <  1:414925/108‑1 (MQ=255)
                                                                                                          cAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCa  <  2:289024/139‑1 (MQ=255)
                                                                                                          cAGATACGCTTTTCCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCa  <  2:558064/139‑1 (MQ=255)
                                                                                                                       |                                                                                                                             
CGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCCCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGTTTCGCTGACCATTACCCGCAAAAAAGTGTCCTGGGCA  >  NZ_CP009273/4507686‑4507930

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: