Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 700,805 A→T C264S (TGC→AGC)  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_000913700,8050AT97.9% 166.7 / ‑3.5 48C264S (TGC→AGC) nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (22/25);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GTGCAGTCGTCCAGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCC  >  NC_000913/700708‑700900
                                                                                                 |                                                                                               
gTGCAGTCGTCCAGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACCGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGc                                                                                               >  2:861133/1‑100 (MQ=255)
   cAGTCGTCCAGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCaa                                                                                            >  1:50500/1‑100 (MQ=255)
   cAGTCGTCCAGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCaa                                                                                            >  1:1239350/1‑100 (MQ=255)
     gtcgtcCAGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAg                                                                                          <  1:1131742/100‑1 (MQ=255)
        gtcCAGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTg                                                                                       >  1:1255168/1‑100 (MQ=255)
               ggCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGt                                                                                <  1:426285/100‑1 (MQ=255)
                gCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTg                                                                               <  1:784853/100‑1 (MQ=255)
                    gcgGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAg                                                                           <  1:226710/100‑1 (MQ=255)
                       gCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGt                                                                        <  2:1100639/100‑1 (MQ=255)
                        ctctGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGtt                                                                       >  1:585849/1‑100 (MQ=255)
                            ggTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCAcc                                                                   >  1:872939/1‑100 (MQ=255)
                            ggTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCAcc                                                                   >  2:1244031/1‑100 (MQ=255)
                              tAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCa                                                                 <  1:1244042/100‑1 (MQ=255)
                                gCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGc                                                               >  1:906664/1‑100 (MQ=255)
                                  ccTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCgg                                                             >  1:174308/1‑100 (MQ=255)
                                  ccTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCgg                                                             <  1:77428/100‑1 (MQ=255)
                                   cTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGt                                                            >  1:989010/1‑100 (MQ=255)
                                    tGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGtt                                                           >  1:186586/1‑100 (MQ=255)
                                       tttAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGa                                                        <  1:751051/100‑1 (MQ=255)
                                        ttAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGAc                                                       <  2:906657/100‑1 (MQ=255)
                                         tAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGAcc                                                      >  1:1213871/1‑100 (MQ=255)
                                          aaCAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCt                                                     <  2:1042765/100‑1 (MQ=255)
                                             aGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGaa                                                  <  1:415668/100‑1 (MQ=255)
                                             aGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGaa                                                  <  1:389805/100‑1 (MQ=255)
                                             aGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGaa                                                  <  1:6360/100‑1 (MQ=255)
                                              gATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAat                                                 >  1:90761/1‑100 (MQ=255)
                                               aTTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCGAGGCTGCCAAAGTTGCCTCAGTGGCAGCGGTCACCCAGCGGGTCGACCTTAat                                                 >  2:106970/1‑99 (MQ=255)
                                               aTTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAata                                                >  2:886230/1‑100 (MQ=255)
                                               aTTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAata                                                >  2:1237070/1‑100 (MQ=255)
                                                  cAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATgg                                             <  1:578502/100‑1 (MQ=255)
                                                  cAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATgg                                             >  2:105630/1‑100 (MQ=255)
                                                   aaCACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGc                                            <  2:17757/100‑1 (MQ=255)
                                                    acacCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGcc                                           <  1:715500/100‑1 (MQ=255)
                                                    acacCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGcc                                           <  2:491408/100‑1 (MQ=255)
                                                       cccGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAAt                                        <  1:322352/100‑1 (MQ=255)
                                                         cGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAAttt                                      <  1:507425/100‑1 (MQ=255)
                                                               tCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGa                                >  1:157428/1‑100 (MQ=255)
                                                                cAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGAc                               >  1:316545/1‑100 (MQ=255)
                                                                       aGCGTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGAcgttgccg                        <  2:25848/100‑1 (MQ=255)
                                                                          gTTGGCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGAcgttgccgttg                     >  1:290862/1‑100 (MQ=255)
                                                                              gCAGCGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGc                 <  1:1202142/100‑1 (MQ=255)
                                                                                  cGATAGTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAAt             <  2:174306/100‑1 (MQ=255)
                                                                                     tcgctcccAGGCTGCTCAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGGTCGTCCTGAATAGGGCCAAGTTCAACGGCGGTGCCGTTGCGGCCAATaaa          <  2:343507/94‑1 (MQ=255)
                                                                                       gTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATaaaaa        <  1:1289779/100‑1 (MQ=255)
                                                                                       gTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATaaaaa        >  1:252893/1‑100 (MQ=255)
                                                                                       gTTTCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATaaaaa        <  1:1027033/100‑1 (MQ=255)
                                                                                         ttCCAGGCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATg      <  2:585845/100‑1 (MQ=255)
                                                                                              ggCTGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGcc  >  1:231402/1‑99 (MQ=255)
                                                                                                 |                                                                                               
GTGCAGTCGTCCAGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCC  >  NC_000913/700708‑700900

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