Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,491 T→C D192G (GAT→GGT)  yjjI ← DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,4910TC89.1% 154.1 / 9.0 55D192G (GAT→GGT) yjjIDUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (5/1);  new base C (28/21);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

ATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGG  >  NC_000913/4616397‑4616587
                                                                                              |                                                                                                
aTCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAg                                                                                             <  2:649102/100‑1 (MQ=255)
 tCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGa                                                                                            <  1:198460/100‑1 (MQ=255)
   aTGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATa                                                                                          >  1:670664/1‑100 (MQ=255)
     gCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATaaa                                                                                        >  2:457783/1‑100 (MQ=255)
      cACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAg                                                                                       <  1:1129251/100‑1 (MQ=255)
      cACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAg                                                                                       <  1:1129315/100‑1 (MQ=255)
       aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAgg                                                                                      <  2:725141/100‑1 (MQ=255)
       aCCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAgg                                                                                      <  2:725134/100‑1 (MQ=255)
         cGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGt                                                                                     >  2:402777/1‑99 (MQ=255)
          ggACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCa                                                                                   >  1:924121/1‑100 (MQ=255)
          ggACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCa                                                                                   <  1:689235/100‑1 (MQ=255)
             ccGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGt                                                                                <  1:651776/100‑1 (MQ=255)
                  ggCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGt                                                                           <  1:1276264/100‑1 (MQ=255)
                   gCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTg                                                                          >  2:104781/1‑100 (MQ=255)
                         tgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCt                                                                    >  1:1298532/1‑100 (MQ=255)
                         tgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCt                                                                    >  2:705324/1‑100 (MQ=255)
                         tgCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCt                                                                    >  1:1117576/1‑100 (MQ=255)
                             gTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCtt                                                                >  1:1290754/1‑100 (MQ=255)
                             gTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCtt                                                                >  2:656977/1‑100 (MQ=255)
                              tttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTc                                                               <  2:939124/100‑1 (MQ=255)
                              tttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTc                                                               >  2:749659/1‑100 (MQ=255)
                               ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTcn                                                              <  1:540407/100‑2 (MQ=255)
                               ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCa                                                              >  2:1275523/1‑100 (MQ=255)
                               ttGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCa                                                              >  1:62131/1‑100 (MQ=255)
                                  cTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAAct                                                           >  2:1152174/1‑100 (MQ=255)
                                   tACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAActc                                                          >  1:300265/1‑100 (MQ=255)
                                    aCATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAActct                                                         >  2:1047222/1‑100 (MQ=255)
                                      aTTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgc                                                       >  2:831689/1‑100 (MQ=255)
                                      aTTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgc                                                       >  1:106011/1‑100 (MQ=255)
                                       ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgca                                                      >  2:413993/1‑100 (MQ=255)
                                       ttCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACtctgca                                                      <  1:1033771/100‑1 (MQ=255)
                                                gTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGt                                             >  1:585203/1‑100 (MQ=255)
                                                   cttcGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAg                                          >  2:263805/1‑100 (MQ=255)
                                                       gCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATcg                                      >  2:3659/1‑100 (MQ=255)
                                                           cTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATcgcgcg                                  <  2:1003791/100‑1 (MQ=255)
                                                             tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg                                >  2:1150936/1‑100 (MQ=255)
                                                             tCCAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCggg                                <  2:260905/100‑1 (MQ=255)
                                                              ccAGCAGCAGGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGt                               >  1:754401/1‑100 (MQ=255)
                                                                 gcagcagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAAt                            <  1:522472/100‑1 (MQ=255)
                                                                    gcagGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATggg                         >  2:644244/1‑100 (MQ=255)
                                                                     cagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGc                        <  2:168014/100‑1 (MQ=255)
                                                                     cagGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGc                        <  2:167652/100‑1 (MQ=255)
                                                                      agGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCg                       >  1:577082/1‑100 (MQ=255)
                                                                      agGTCATCAGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCg                       <  1:635465/100‑1 (MQ=255)
                                                                         tcatcaGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAAt                    <  1:543821/100‑1 (MQ=255)
                                                                          catcaGGGGTGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATc                   >  2:91597/1‑100 (MQ=255)
                                                                                gggTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACgg             >  1:238528/1‑100 (MQ=255)
                                                                                   tGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGcc          <  2:55885/100‑1 (MQ=255)
                                                                                   tGATTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGcc          <  1:1150947/100‑1 (MQ=255)
                                                                                     aTTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGa        >  1:522895/1‑100 (MQ=255)
                                                                                     aTTTCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGa        >  1:1219804/1‑100 (MQ=255)
                                                                                      tttCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGat       >  1:126365/1‑100 (MQ=255)
                                                                                        tCAGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGatat     >  2:245722/1‑100 (MQ=255)
                                                                                          aGGACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTg   <  1:749664/100‑1 (MQ=255)
                                                                                           ggACCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTgg  >  1:1232328/1‑100 (MQ=255)
                                                                                              |                                                                                                
ATCATGCACCGGACCGTTGGCGATGTGCGGTTTGCTACATTCACAGATGTTCTTCGCCACTTCCAGCAGCAGGTCATCAGGGGTGATTTCAGGATCGTAGATAAAGGTCAGGTTCGGTGAAACCTGCTTCAACTCTGCATCTGCACGTAAGATCGCGCGGGTAATGGGCGAATCAGACGGGCCGATATTGG  >  NC_000913/4616397‑4616587

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