Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,997,767 T→C G195G (GGT→GGC yigB → 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,997,7670TC96.1% 165.1 / ‑2.8 51G195G (GGT→GGCyigB5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (30/19);  total (31/20)
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.86e-01

TTTGTGCTGCGCGCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAG  >  NC_000913/3997636‑3997903
                                                                                                                                   |                                                                                                                                        
tttGTGCTGCGCGCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAAt                                                                                                                                  >  2:562886/1‑140 (MQ=255)
   gtgCTGCGCGCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCg                                                                                                                               >  1:136730/1‑140 (MQ=255)
   gtgCTGCGCGCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCg                                                                                                                               >  1:626952/1‑140 (MQ=255)
   gtgCTGCGCGCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCg                                                                                                                               >  1:41870/1‑140 (MQ=255)
              tGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGc                                                                                                                    >  2:159985/1‑140 (MQ=255)
                gCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACATGGGCGGGGCAATTCGCAGCGGAATGCAg                                                                                                                  >  1:364709/1‑140 (MQ=255)
                     cACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCttg                                                                                                             >  2:601983/1‑140 (MQ=255)
                           cgcTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGCGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGAt                                                                                                       <  2:592637/140‑1 (MQ=255)
                                aaaaCCGTTCAGCGAAATGTACTTTTTGGTTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGAc                                                                                                  <  1:405727/140‑1 (MQ=255)
                                         cAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATg                                                                                         <  2:41870/140‑1 (MQ=255)
                                                        tttGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCaaa                                                                          >  1:163316/1‑140 (MQ=255)
                                                          tGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAAcc                                                                        >  2:23124/1‑140 (MQ=255)
                                                           ggCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCt                                                                       >  2:662999/1‑140 (MQ=255)
                                                             cTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTgg                                                                     >  2:810846/1‑140 (MQ=255)
                                                               gCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGa                                                                   <  2:656928/140‑1 (MQ=255)
                                                               gCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCAGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGa                                                                   <  1:885645/140‑1 (MQ=255)
                                                               gCGGAAAAACTCAACATGCCGAGCGCCGATATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGa                                                                   <  1:413538/140‑1 (MQ=255)
                                                                  gAAAAACTCAACGTGCCGATCGTCGGGATATTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAg                                                                <  1:406515/140‑1 (MQ=255)
                                                                            aCGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTg                                                      >  2:26609/1‑140 (MQ=255)
                                                                             cGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCGGACCGGAACATGGCGATCTGATGCATACCTTGGACAGCCGTTTACTGc                                                     >  2:406515/1‑140 (MQ=255)
                                                                                gCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGc                                                  >  1:697530/1‑140 (MQ=255)
                                                                                gCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGc                                                  >  2:507650/1‑140 (MQ=255)
                                                                                gCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGc                                                  >  2:66998/1‑140 (MQ=255)
                                                                                   gATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATc                                               <  1:159985/140‑1 (MQ=255)
                                                                                         cGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGCTGCAAACCTGGGACAGGCGTTCACTGCCGCATCTGGaaa                                         >  2:432075/1‑140 (MQ=255)
                                                                                                     tGTTGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTgg                             >  1:518512/1‑140 (MQ=255)
                                                                                                       ttGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTcc                                    >  2:49775/1‑131 (MQ=255)
                                                                                                       ttGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTcc                                    <  1:49775/131‑1 (MQ=255)
                                                                                                       ttGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTgg                             <  2:619332/138‑1 (MQ=255)
                                                                                                       ttGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTgg                             >  1:619332/1‑138 (MQ=255)
                                                                                                        tGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAt                          >  2:600934/1‑140 (MQ=255)
                                                                                                        tGGGGACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCCTTTGCTGCCGCAGCTGGGAATTTCCCGGTTGGCAt                          >  2:265967/1‑140 (MQ=255)
                                                                                                         ggggACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCCGGGACAGTCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAtc                         >  2:576858/1‑140 (MQ=255)
                                                                                                          gggACGATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAtct                        >  2:669710/1‑140 (MQ=255)
                                                                                                               gATCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGAc                   <  2:697530/140‑1 (MQ=255)
                                                                                                                aTCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGAcc                  >  1:946849/1‑140 (MQ=255)
                                                                                                                aTCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGAcc                  >  1:68809/1‑140 (MQ=255)
                                                                                                                aTCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGAc                   <  2:538524/139‑1 (MQ=255)
                                                                                                                aTCTCACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGAc                   >  1:538524/1‑139 (MQ=255)
                                                                                                                 tctcACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGa                    >  1:696476/1‑137 (MQ=255)
                                                                                                                 tctcACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGa                    <  2:696476/137‑1 (MQ=255)
                                                                                                                  ctcACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTc                <  1:600934/140‑1 (MQ=255)
                                                                                                                   tcACCACTGACGTGGGCTGGGCCATTCGCAGAGGATTGCGGGCTTGTTCGATCAGACCGGAAAATCGCGATCTGATGCAAACCTGGAACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGATTGGCATCTCTGACCTCg               <  2:497261/140‑1 (MQ=255)
                                                                                                                   tcACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCg               <  1:810846/140‑1 (MQ=255)
                                                                                                                   tcACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCg               <  2:626952/140‑1 (MQ=255)
                                                                                                                   tcACCACTGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCg               <  2:136730/140‑1 (MQ=255)
                                                                                                                      ccacTGACGTGGGCGTGGCAATTCGTAGCGGAATGCAAGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTg            <  1:432075/140‑1 (MQ=255)
                                                                                                                         cTGACGTGGGCGGGGCATTTCGCAGCGGAATGCAGACTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGata         <  2:408190/140‑1 (MQ=255)
                                                                                                                          tGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGat          >  2:930259/1‑138 (MQ=255)
                                                                                                                          tGACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGat          <  1:930259/138‑1 (MQ=255)
                                                                                                                           gACGTGGGCGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGAtata       >  1:789829/1‑140 (MQ=255)
                                                                                                                           gACGTGGGCGGGGCAATTCGCAGCGGAAAGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGAtata       >  2:611133/1‑140 (MQ=255)
                                                                                                                            aCGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATaa      >  2:142806/1‑140 (MQ=255)
                                                                                                                               tgggcgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAAATTCCCGGTTGGTATCTCTGACCTCGCTGATATAATCa   >  1:319727/1‑140 (MQ=255)
                                                                                                                                gggcgggGCAATTCGCAGCGGAAAGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAg  <  1:611133/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                        
TTTGTGCTGCGCGCTGGCCCGCACGGGCGCTCAAAACCGTTCAGCGATATGTACTTTTTGGCTGCGGAAAAACTCAACGTGCCGATCGGCGAGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAG  >  NC_000913/3997636‑3997903

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: