Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 700,948 A→C I216S (ATT→AGT)  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,9480AC98.0% 164.7 / NA 52I216S (ATT→AGT) nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (23/27);  total (24/28)
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.99e-01

CAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTT  >  NC_000913/700808‑701072
                                                                                                                                            |                                                                                                                            
cAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTg                                                                                                                            <  1:114069/143‑1 (MQ=255)
  acgctGCCGCAGTTGGAGGGTGCCCCCGGCGGGTCGACCTGAATATGGCCAATTTCCCCGACGTTGCCGTTGTGGCCAATAAAAATGCGTCCGTTAGAGATAGTCCGGGCACCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  2:273087/139‑1 (MQ=255)
  aaGTTGTCGCAGTGGCAGCTTTCACCCAGCGGTTCGGCCTGAAGATGGCCAATTTCACCGACGTTGCCGTGGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  2:834041/143‑1 (MQ=255)
  aaGTTGCCGCAGTTGCAGCGTTCGCCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCAAATATAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  2:516119/143‑1 (MQ=255)
  aaGCTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTGCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  2:299305/143‑1 (MQ=255)
    gTTGCCGCGGTGGCAGCGTTCACACAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGAGGCCAATAAAAATAGGCCCGTTAGAGATAATCCCGGAACCGGTTCGGCGATGTACAGGCACCAGAGTgtcga                                                                                                                        >  2:32906/1‑139 (MQ=255)
         cGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGCGATAATCCCGGCCCCGGTTCCGCGATGGGCACGCGCCAGACTGGACGCTTCg                                                                                                                   >  1:320419/1‑143 (MQ=255)
              tGGCGGCGTTCACCCGGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATc                                                                                                              <  2:655253/143‑1 (MQ=255)
                gCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACCCACCAGACTGGCGTCTTCGCAATCCt                                                                                                            >  1:196744/1‑143 (MQ=255)
                  aGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGGTAATCCCGGCCCCGGTTCCGCGATGGACCCGCACCAGACTGGAGTCTTCGCAATCCTGa                                                                                                          >  1:416100/1‑143 (MQ=255)
                  aGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAACGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGa                                                                                                          <  1:744059/143‑1 (MQ=255)
                   gCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAAGGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGAc                                                                                                         <  1:555982/143‑1 (MQ=255)
                   cctctCACCCAGCGGTCCGACCTGAATATGGCCAATTTCCCCGACGTTGCCGTTGCGGCCAATAAAAATGCGCACGATAAAGAGAATCCCTGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGAc                                                                                                         <  1:400682/139‑1 (MQ=255)
                     gTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACtt                                                                                                       >  1:418102/1‑143 (MQ=255)
                           ccAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCCATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGCGTCCTCGCAATCCTGACTTGccccc                                                                                                 >  2:783303/1‑142 (MQ=255)
                           ccAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGTCCCAGGTTCCGCGATGGACACGCACCAGACGGGAGGCTTCGCAATCCTGACTTGCCCCg                                                                                                 >  1:949514/1‑143 (MQ=255)
                             aGCGGTTCGACCTGTATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGaa                                                                                               <  1:17036/143‑1 (MQ=255)
                             aGCGGTGCGACCTGAAGATGGCCACTGTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGaa                                                                                               <  1:982414/143‑1 (MQ=255)
                                gtttccgccTTAATATGTCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAagta                                                                                            <  2:74012/136‑1 (MQ=255)
                                gtttCGCCCTGAATATGGCCAATGTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAagta                                                                                            <  2:751135/141‑1 (MQ=255)
                                ggTTCGACCTGAATATGGCCAATTTCACCGACGTGGCCGTTGCGGCCAATAAAAAGGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAagta                                                                                            <  1:431747/143‑1 (MQ=255)
                                       gctgggtatGGCTAGTTTGGCCGTGGTGGCCGGTGCGGCCAATCAAATGGCGCGCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTcc                                                                                     <  2:9986/137‑1 (MQ=255)
                                        cTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCAAATAAAAATGCTCTCGGGAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTTcg                                                                                    >  1:833901/1‑143 (MQ=255)
                                         tGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgc                                                                                   >  1:765980/1‑143 (MQ=255)
                                             tatGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGCTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCACTCCCGACTTGCACCGAGGGAGTGCCCcgccagc                                                                               >  2:82021/1‑143 (MQ=255)
                                             tatGGCCAACTTCACCGACGGTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagc                                                                               <  1:315295/143‑1 (MQ=255)
                                                 gCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgcca                                                                           <  1:317320/143‑1 (MQ=255)
                                                 gCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGCGCCCcgccagcccca                                                                           >  2:999616/1‑143 (MQ=255)
                                                  ccAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgccag                                                                          <  1:304745/143‑1 (MQ=255)
                                                        tCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACg                                                                    <  1:567932/143‑1 (MQ=255)
                                                          aCCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGa                                                                  >  1:910096/1‑143 (MQ=255)
                                                           ccGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGac                                                                 >  1:801891/1‑142 (MQ=255)
                                                            cGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGGCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGGCTACGGcta                                                                >  2:88226/1‑143 (MQ=255)
                                                               cgttgccgttgcGGCCAATAAAAATGCGCCCGTTAGAGACAATCCCGGCCCCGGCTCCGCGATGGACACGCACCAGACTGGAGTCTGCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCg                                                             <  1:547484/143‑1 (MQ=255)
                                                                     cgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCACTCCTGACTTGCACCCAAGTACTGCTCCGCCAGCGCTGTACTAGGGATATCCTGACCa                                                       >  2:258507/1‑143 (MQ=255)
                                                                     cgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGAAACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTCCGGATATCGCGACCa                                                       >  1:375888/1‑143 (MQ=255)
                                                                        tgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTGCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACg                                                    <  1:45780/143‑1 (MQ=255)
                                                                           ggCGAATAAAAACGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACATGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGaaa                                                 <  2:323346/143‑1 (MQ=255)
                                                                                aTAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGt                                            >  1:92984/1‑143 (MQ=255)
                                                                                                  gagaTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttcttc                          <  1:891105/143‑1 (MQ=255)
                                                                                                  gagaTAACCCCGGCCCCGGATCCGCGATGGACACGCACCAGACTGGAGTCTGCGCAATCCTGACCTGCACCCAAGTAGTGCTCCGCCAGCGCCATACTACGGATATCGTGACCAACGAAACACGTCACTTTAAAACGttcttc                          <  1:108279/143‑1 (MQ=255)
                                                                                                   agaTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGGCACTTTAAAACGTTCTTcc                         >  1:778060/1‑143 (MQ=255)
                                                                                                    tatcaTCCCGGCTCCGGATCTGCGATGGACACGAGGCAGACTGGAGTCTTCCCAATCCTGACTTGCACTTCAGTAGTGCTCCGCCAGCGCCATACTACGGATATCGTGACCATCGAAACAGGTCACTTTAAAACGTTCTTCCa                        <  2:71126/139‑1 (MQ=255)
                                                                                                    cagacccccGGGCCCGGCTCCGCGCTGGACACCCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCa                        <  1:464609/137‑1 (MQ=255)
                                                                                                     aTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCTGCCAGCGCCAGACTACGGGTATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCag                       >  1:806994/1‑143 (MQ=255)
                                                                                                         tCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGCGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGCGCt                   >  1:497210/1‑143 (MQ=255)
                                                                                                         tCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGCGCt                   >  1:928280/1‑143 (MQ=255)
                                                                                                           ccGCCCCCGGTGCCGCGATGGACACGCACCAGAATGGATTCTTCGCAACACACACTTATAGGTGTGGAGTGCTCCGCCAGCGTCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTc                 <  1:774600/143‑1 (MQ=255)
                                                                                                              gCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTAc              <  2:348068/143‑1 (MQ=255)
                                                                                                                cccGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAACAACGAAACAGGTCACTTTAAAACGTTCTTACAGAGATTCGACCa            >  1:477196/1‑143 (MQ=255)
                                                                                                                        cgcgGTGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCCGACTACGGATATCGTGACCAACGAAACAGGTCACTGTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCac    >  2:291296/1‑142 (MQ=255)
                                                                                                                          cgATGGACACGCACCAGACTAGAGTCTTCCCAAACCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATTGTGACCAACGACAAAAGTCACCTTAAAACGCTGTGACAGAGCTACTACCAGCCcacacgt  >  2:640371/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                            
CAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTT  >  NC_000913/700808‑701072

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: ‑