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.1% 171.1 / NA 54I216S (ATT→AGT) nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (24/28);  total (25/29)
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                                                                                                                            <  2:114069/143‑1 (MQ=255)
  acgctGCCGCAGTTGGAGGGTGCCCCCGGCGGGTCGACCTGAATATGGCCAATTTCCCCGACGTTGCCGTTGTGGCCAATAAAAATGCGTCCGTTAGAGATAGTCCGGGCACCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  1:273087/139‑1 (MQ=255)
  aaGTTGTCGCAGTGGCAGCTTTCACCCAGCGGTTCGGCCTGAAGATGGCCAATTTCACCGACGTTGCCGTGGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  1:834041/143‑1 (MQ=255)
  aaGTTGCCGCAGTTGCAGCGTTCGCCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCAAATATAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  1:516119/143‑1 (MQ=255)
  aaGCTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTGCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGa                                                                                                                          <  1:299305/143‑1 (MQ=255)
    gTTGCCGCGGTGGCAGCGTTCACACAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGAGGCCAATAAAAATAGGCCCGTTAGAGATAATCCCGGAACCGGTTCGGCGATGTACAGGCACCAGAGTgtcga                                                                                                                        >  1:32906/1‑139 (MQ=255)
         cGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGCGATAATCCCGGCCCCGGTTCCGCGATGGGCACGCGCCAGACTGGACGCTTCg                                                                                                                   >  2:320419/1‑143 (MQ=255)
              tGGCGGCGTTCACCCGGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATc                                                                                                              <  1:655253/143‑1 (MQ=255)
                gCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACCCACCAGACTGGCGTCTTCGCAATCCt                                                                                                            >  2:196744/1‑143 (MQ=255)
                  aGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGGTAATCCCGGCCCCGGTTCCGCGATGGACCCGCACCAGACTGGAGTCTTCGCAATCCTGa                                                                                                          >  2:416100/1‑143 (MQ=255)
                  aGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAACGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGa                                                                                                          <  2:744059/143‑1 (MQ=255)
                   gCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAAGGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGAc                                                                                                         <  2:555982/143‑1 (MQ=255)
                   cctctCACCCAGCGGTCCGACCTGAATATGGCCAATTTCCCCGACGTTGCCGTTGCGGCCAATAAAAATGCGCACGATAAAGAGAATCCCTGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGAc                                                                                                         <  2:400682/139‑1 (MQ=255)
                     gTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACtt                                                                                                       >  2:418102/1‑143 (MQ=255)
                           ccAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCCATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGCGTCCTCGCAATCCTGACTTGccccc                                                                                                 >  1:783303/1‑142 (MQ=255)
                           ccAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGTCCCAGGTTCCGCGATGGACACGCACCAGACGGGAGGCTTCGCAATCCTGACTTGCCCCg                                                                                                 >  2:949514/1‑143 (MQ=255)
                             aGCGGTTCGACCTGTATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGaa                                                                                               <  2:17036/143‑1 (MQ=255)
                             aGCGGTGCGACCTGAAGATGGCCACTGTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGaa                                                                                               <  2:982414/143‑1 (MQ=255)
                                gtttccgccTTAATATGTCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAagta                                                                                            <  1:74012/136‑1 (MQ=255)
                                gtttCGCCCTGAATATGGCCAATGTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAagta                                                                                            <  1:751135/141‑1 (MQ=255)
                                ggTTCGACCTGAATATGGCCAATTTCACCGACGTGGCCGTTGCGGCCAATAAAAAGGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAagta                                                                                            <  2:431747/143‑1 (MQ=255)
                                       gctgggtatGGCTAGTTTGGCCGTGGTGGCCGGTGCGGCCAATCAAATGGCGCGCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTcc                                                                                     <  1:9986/137‑1 (MQ=255)
                                        cTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCAAATAAAAATGCTCTCGGGAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTTcg                                                                                    >  2:833901/1‑143 (MQ=255)
                                         tGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgc                                                                                   >  2:765980/1‑143 (MQ=255)
                                             tatGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGCTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCACTCCCGACTTGCACCGAGGGAGTGCCCcgccagc                                                                               >  1:82021/1‑143 (MQ=255)
                                             tatGGCCAACTTCACCGACGGTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagc                                                                               <  2:315295/143‑1 (MQ=255)
                                                 gCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgcca                                                                           <  2:317320/143‑1 (MQ=255)
                                                 gCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGCGCCCcgccagcccca                                                                           >  1:999616/1‑143 (MQ=255)
                                                  ccAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgccag                                                                          <  2:304745/143‑1 (MQ=255)
                                                        tCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACg                                                                    <  2:567932/143‑1 (MQ=255)
                                                          aCCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGa                                                                  >  2:910096/1‑143 (MQ=255)
                                                           ccGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGac                                                                 >  2:801891/1‑142 (MQ=255)
                                                            cGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGGCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGGCTACGGcta                                                                >  1:88226/1‑143 (MQ=255)
                                                               cgttgccgttgcGGCCAATAAAAATGCGCCCGTTAGAGACAATCCCGGCCCCGGCTCCGCGATGGACACGCACCAGACTGGAGTCTGCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCg                                                             <  2:547484/143‑1 (MQ=255)
                                                                     cgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCACTCCTGACTTGCACCCAAGTACTGCTCCGCCAGCGCTGTACTAGGGATATCCTGACCa                                                       >  1:258507/1‑143 (MQ=255)
                                                                     cgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGAAACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTCCGGATATCGCGACCa                                                       >  2:375888/1‑143 (MQ=255)
                                                                        tgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTGCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACg                                                    <  2:45780/143‑1 (MQ=255)
                                                                           ggCGAATAAAAACGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACATGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGaaa                                                 <  1:323346/143‑1 (MQ=255)
                                                                                ataattattcgcCCGTGAGAGATAATCCCGGCCCCGTTTCCGCCATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGt                                            <  1:255050/134‑1 (MQ=255)
                                                                                aTAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGt                                            >  2:92984/1‑143 (MQ=255)
                                                                                                  gagaTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttcttc                          <  2:891105/143‑1 (MQ=255)
                                                                                                  gagaTAACCCCGGCCCCGGATCCGCGATGGACACGCACCAGACTGGAGTCTGCGCAATCCTGACCTGCACCCAAGTAGTGCTCCGCCAGCGCCATACTACGGATATCGTGACCAACGAAACACGTCACTTTAAAACGttcttc                          <  2:108279/143‑1 (MQ=255)
                                                                                                   agaTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGGCACTTTAAAACGTTCTTcc                         >  2:778060/1‑143 (MQ=255)
                                                                                                    tatcaTCCCGGCTCCGGATCTGCGATGGACACGAGGCAGACTGGAGTCTTCCCAATCCTGACTTGCACTTCAGTAGTGCTCCGCCAGCGCCATACTACGGATATCGTGACCATCGAAACAGGTCACTTTAAAACGTTCTTCCa                        <  1:71126/139‑1 (MQ=255)
                                                                                                    cagacccccGGGCCCGGCTCCGCGCTGGACACCCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCa                        <  2:464609/137‑1 (MQ=255)
                                                                                                     aTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCTGCCAGCGCCAGACTACGGGTATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCag                       >  2:806994/1‑143 (MQ=255)
                                                                                                         tCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGCGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGCGCt                   >  2:497210/1‑143 (MQ=255)
                                                                                                         tCCCGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGCGCt                   >  2:928280/1‑143 (MQ=255)
                                                                                                           ccGGCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGCCTTGCACCGACGTCCTGCTCCGCCAGCGCCAGACTACCGATAAACTCACCAACGCCACAGGTCACCTAAAAACGtccttgaaaacctcc                 >  1:367276/1‑133 (MQ=255)
                                                                                                           ccGCCCCCGGTGCCGCGATGGACACGCACCAGAATGGATTCTTCGCAACACACACTTATAGGTGTGGAGTGCTCCGCCAGCGTCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTc                 <  2:774600/143‑1 (MQ=255)
                                                                                                              gCCCCGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTAc              <  1:348068/143‑1 (MQ=255)
                                                                                                                cccGGTTCCGCGATGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAACAACGAAACAGGTCACTTTAAAACGTTCTTACAGAGATTCGACCa            >  2:477196/1‑143 (MQ=255)
                                                                                                                        cgcgGTGGACACGCACCAGACTGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCCGACTACGGATATCGTGACCAACGAAACAGGTCACTGTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCac    >  1:291296/1‑142 (MQ=255)
                                                                                                                          cgATGGACACGCACCAGACTAGAGTCTTCCCAAACCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATTGTGACCAACGACAAAAGTCACCTTAAAACGCTGTGACAGAGCTACTACCAGCCcacacgt  >  1: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: ‑