Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 701,378 +T coding (217/1221 nt) nagC ← N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913701,3781.T83.8% 84.3 / 12.7 38coding (217/1221 nt)nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base T (14/17);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGG  >  NC_000913/701243‑701501
                                                                                                                                             |                                                                                                                           
gCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGTGCCCCCCGt                                                                                                                       >  2:665661/1‑143 (MQ=255)
        tAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCct                                                                                                               >  1:281746/1‑143 (MQ=255)
             ttgttctGCCACCCCTTTGCTGCGCTGATCACACAGCGTGATGGTGGCGTCATGACGCCCGCTCCGTACGCCGCTTGCGTCGTCATTGCAGGTTTCGGTGACGATGGATATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGa                                                                                                          <  2:714238/140‑1 (MQ=255)
                 tctGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCaa                                                                                                      <  1:372895/143‑1 (MQ=255)
                     ccAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAActtc                                                                                                  >  1:258898/1‑143 (MQ=255)
                           ttttGCTGCTGAGATCAAACAGAGCGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGAt                                                                                            <  1:3929/142‑1 (MQ=255)
                           cTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGAt                                                                                            <  2:304622/143‑1 (MQ=255)
                            tttGCTGCTGAGATCAAACAGAGTGAAGCTAGCGTCATGTCGACCAAGCCGTACGCCGATTGCGAGGAAAGTTGCGGATCGGAGGACGAGGGAGATAGCGCGACGGCCTCCCGGTGGAGGCCTGCTGAGCAACGTCTGTGATc                                                                                           >  2:398454/1‑143 (MQ=255)
                                 tgctgAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGccc                                                                                      >  1:213733/1‑143 (MQ=255)
                                       gaTCAAACTGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTc                                                                                <  1:708555/143‑1 (MQ=255)
                                        aTCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGACATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCTCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTCCg                                                                               >  1:322690/1‑143 (MQ=255)
                                        aTCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCg                                                                               >  1:107721/1‑143 (MQ=255)
                                         tCAAACAGAGTGATGGTGGCGTCATGACGACCAAGACGTAAGCCGATTGGGTGGAAATTGCGGGTTTCGGTGAAGATGGAGATAGCGCGGCGGCATCCCGGTGGAGGCCTGCTGATCAACTTCTTTGCTCAGCGCGCGCTCGa                                                                              >  2:557497/1‑143 (MQ=255)
                                             aCAGAGGGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTGTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAg                                                                          <  1:35673/143‑1 (MQ=255)
                                                  gtggtggtgGCGTCGTGAGGGCCAAGCTGTGGGGCGATTGGGTGGAGATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGAc                                                                     <  2:458038/143‑1 (MQ=255)
                                                   tggtggtggGGTCAGGCCGGCGCAGGCGTACGCGGAGTGCGTGGAAAGTGCGGGTTTCGGTGACGAGGGAGATAGCGGGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGatg                                                                    <  2:650250/140‑3 (MQ=255)
                                                   tGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATACGCTGAcg                                                                    >  1:583283/1‑143 (MQ=255)
                                                     aTGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGGTCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGAcgcg                                                                  >  2:407501/1‑143 (MQ=255)
                                                              tCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTgg                                                         >  2:94261/1‑143 (MQ=255)
                                                               cATGACGGCCAAGCTGTACGCCGATGGTGTTGAGAGTGCGGGTTGCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGCTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGt                                                        <  2:41970/143‑1 (MQ=255)
                                                                aTGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGACGGCCTCCCGGTGGAGGCCGGCTGATCAACTACTTTGATAAGCCCCGGTACGATAAGACTACGCGAAAATCTGGtg                                                       >  2:611474/1‑142 (MQ=255)
                                                                 tGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGGCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTGCGGTAGGCTGACGCGTCATTTTGGTTa                                                      >  2:39472/1‑143 (MQ=255)
                                                                 tGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTa                                                      >  1:22494/1‑143 (MQ=255)
                                                                  gacgacCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTAc                                                     <  1:477632/143‑1 (MQ=255)
                                                                       cccaGCCGGCCGCCGACTGCGTGGCAAGTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTgg                                                <  2:419853/140‑1 (MQ=255)
                                                                         aaGCCGTACGCCGATTGCGGGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGAGCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCgg                                             <  1:214642/143‑1 (MQ=255)
                                                                          aGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCgg                                             <  1:667587/143‑1 (MQ=255)
                                                                                aCGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCaa                                       >  1:331442/1‑143 (MQ=255)
                                                                                 cGCCGATGGCGTGGAAGTTGCGGGTTTCGGGGACGGTGGAGATAGCGCGGCGGCGTCCCGGTGGAGGCCCGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAg                                      <  1:399891/143‑1 (MQ=255)
                                                                                     gATTGCGTGGGAATTGCGGGTTGCGGTGGCGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATGAGCGCGCGTTGGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTgg                                  <  2:538698/143‑1 (MQ=255)
                                                                                           gTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTCCGCTGGCGGGGGCAAGATGGCTTTGc                            >  2:110514/1‑143 (MQ=255)
                                                                                              gacaTTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCg                         <  1:236452/140‑1 (MQ=255)
                                                                                                acgtggCGGGTTCCGGTGACGTTGGAGATAGCGCGGCCGCCTCCCGGTGGAGGCCTGCTGATCAACTTCGTTGATCAGCCCGCGTTCGACAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCACCCTGGCTGGGCTCgg                        <  2:525112/138‑1 (MQ=255)
                                                                                                     gCGGGTTTCGGGGACGATGGAGATAGCGCGGCGGCCGCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATc                   <  2:8962/143‑1 (MQ=255)
                                                                                                                tGACGGTGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCGGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCga        <  2:258898/143‑1 (MQ=255)
                                                                                                                    gATGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTGCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTGCGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGGtt    >  2:633926/1‑143 (MQ=255)
                                                                                                                     aTGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTg   <  2:107721/143‑1 (MQ=255)
                                                                                                                      tGGAGATAGCGCGGCGGCCTCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATggg  >  2:495192/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                           
GCAGCGGGTAATGTTCTTCTGCCAGCACTTTGCTGCTGAGATCAAACAGAGTGATGGTGGCGTCATGACGACCAAGCCGTACGCCGATTGCGTGGAAATTGCGGGTTTCGGTGACGATGGAGATAGCGCGGCGGCCCCCGGTGGAGGCCTGCTGATCAACTTCTTTGATCAGCCCGCGTTCGATAAGCTGACGCGTAATTTTGGTTACGCTGGCGGGGGCAAGCTGGCTTTGCTCGGCAATCTGAATCCGCGAGATTGG  >  NC_000913/701243‑701501

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: