Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,615,011 (C)5→6 coding (533/864 nt) yjjW ← putative glycyl‑radical enzyme activating enzyme YjjW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,615,0061.C90.2% 131.0 / 7.4 41coding (538/864 nt)yjjWputative glycyl‑radical enzyme activating enzyme YjjW
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (18/19);  total (20/21)
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.09e-01

GTGTTGCAAATAATCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAG  >  NC_000913/4614872‑4615140
                                                                                                                                        |                                                                                                                                       
gtgtTGCAAATAATCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTccccc                                                                                                                                   >  1:25634/1‑140 (MQ=255)
 tgtTGCAAATAATCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTcccccc                                                                                                                                  >  1:726052/1‑139 (MQ=255)
  gtTGCAAATAATCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTcccccca                                                                                                                                 <  1:166101/140‑3 (MQ=255)
  gtTGCAAATAATCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTcccccca                                                                                                                                 <  1:77312/140‑3 (MQ=255)
    tGCAAATAATCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTccccccacg                                                                                                                               <  1:196799/140‑5 (MQ=255)
         ataatCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTg                                                                                                                           <  2:117461/140‑1 (MQ=255)
            atCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTgaga                                                                                                                        <  1:900684/140‑1 (MQ=255)
             tCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGAt                                                                                                                       >  2:61981/1‑140 (MQ=255)
                aCCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGCTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGa                                                                                                                    >  1:671212/1‑140 (MQ=255)
                      ccAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTg                                                                                                              <  2:371013/140‑1 (MQ=255)
                         gggATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTgcgc                                                                                                           >  2:329170/1‑140 (MQ=255)
                          ggATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGcc                                                                                                          >  1:829362/1‑140 (MQ=255)
                               aCCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGc                                                                                                     >  2:993894/1‑140 (MQ=255)
                               aCCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGc                                                                                                     >  2:455055/1‑140 (MQ=255)
                               aCCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGc                                                                                                     >  2:244997/1‑140 (MQ=255)
                                      aaCGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCgg                                                                                              >  1:193741/1‑140 (MQ=255)
                                       aCGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCggg                                                                                             >  1:423909/1‑140 (MQ=255)
                                                           cgcgCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGt                                                                         >  2:272847/1‑140 (MQ=255)
                                                            gcgcTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGtt                                                                        <  1:206134/140‑1 (MQ=255)
                                                               ctctGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCg                                                                     >  1:646145/1‑140 (MQ=255)
                                                                      aGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACa                                                              <  1:314292/140‑1 (MQ=255)
                                                                         aaaTAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGc                                                           <  2:423909/140‑1 (MQ=255)
                                                                           aTAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCCCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCg                                                         <  1:61981/140‑1 (MQ=255)
                                                                            tAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCCCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGt                                                        >  1:555915/1‑140 (MQ=255)
                                                                                tgctgcGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACt                                                    <  1:125134/140‑1 (MQ=255)
                                                                                      gcTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTccac                                              >  1:252476/1‑140 (MQ=255)
                                                                                       cTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTccacc                                             >  2:793994/1‑140 (MQ=255)
                                                                                         tAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAg                                           <  1:455055/140‑1 (MQ=255)
                                                                                                 tgATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCCCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTga                                   <  2:555915/140‑1 (MQ=255)
                                                                                                        cgcgTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCCCCCCGGGAGAAATTTTCCCCAGCCGGTTTCGCTCACCATGCCGTTACTGTCCACCAGACAGGTGAGATGgcg                            <  2:193741/140‑1 (MQ=255)
                                                                                                        cgcgTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGgcg                            >  1:295190/1‑140 (MQ=255)
                                                                                                        cgcgTCCGGTGAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGgcg                            <  2:829362/140‑1 (MQ=255)
                                                                                                                  gAGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGa                  <  1:890345/140‑1 (MQ=255)
                                                                                                                   aGTTGTTGATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGAt                 <  1:272847/140‑1 (MQ=255)
                                                                                                                       gttgATGACATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGtt             >  1:741345/1‑140 (MQ=255)
                                                                                                                           atgaCATTCGCTCCTCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTg         >  1:392304/1‑140 (MQ=255)
                                                                                                                            tgaCATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGa        <  1:993894/140‑1 (MQ=255)
                                                                                                                            tgaCATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGa        <  2:741345/140‑1 (MQ=255)
                                                                                                                            tgaCATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGa        <  1:17634/140‑1 (MQ=255)
                                                                                                                            tgaCATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGa        <  1:341041/140‑1 (MQ=255)
                                                                                                                               cATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGa        >  2:587906/1‑137 (MQ=255)
                                                                                                                               cATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGa        <  1:587906/137‑1 (MQ=255)
                                                                                                                               cATTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAg     <  1:244997/140‑1 (MQ=255)
                                                                                                                                aTTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCg               >  1:803330/1‑129 (MQ=255)
                                                                                                                                aTTCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCg               <  2:803330/129‑1 (MQ=255)
                                                                                                                                  tCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAg  >  1:339311/1‑140 (MQ=255)
                                                                                                                                  tCGCTCCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAg  <  2:252476/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GTGTTGCAAATAATCCACCTGGCCAGGGATCACCAGCAAACGCAGTTCCGCCAGCTTGCCGCGCTCTGCCAGCAAATAGATGCTGCGCTTAATCTGCTGATTATCGCGTCCGGTGAGTTGTTGATGACATTCGCTCCCCCACGCTTTGAGATCGAGCATTGCGCCGTCGCACACCGGGAGCAATTTTTCCCAGCCGGTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTGATAGCAG  >  NC_000913/4614872‑4615140

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