Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,237,401 G→A V31I (GTA→ATA)  ygjR → putative oxidoreductase YgjR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,237,4010GA96.9% 107.8 / ‑3.8 32V31I (GTA→ATA) ygjRputative oxidoreductase YgjR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/13);  minor base T (0/1);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCTACATTGACGACAGAATCCCTTTATGGAGTATCCACGCGTTATGATACGTTTCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCC  >  NC_000913/3237265‑3237533
                                                                                                                                        |                                                                                                                                    
cGTGCTACATTGACGACAGAATCCCTTTATGGAGTATCCACGCGTTATGATACGTTTCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCatat                                                                                                                                   >  2:310149/1‑140 (MQ=255)
            acgacAGAATCCCTTTATGGAGTATCCACGCGTTATGATACGTTTCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGcc                                                                                                                       >  2:783644/1‑140 (MQ=255)
                 aGAATCCCTTTATGGAGTATCCACGCGTTATGATACGTTTCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGaa                                                                                                                  >  1:624864/1‑140 (MQ=255)
                                      cACGCGTTATGATGCGTTTCGCTGTGATTGGTACGAACTGGATCCCCCGCCAGTTCGTCGAGGCCGCCCATGAGGGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAAt                                                                                             <  2:879351/140‑1 (MQ=255)
                                       aCGCGTTATGATACGTTTCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATg                                                                                            >  1:539788/1‑140 (MQ=255)
                                                  tACGTTTCGCTGTGATTGGTACGAACTGGATCCCTCGCCAGTTCGTCGGGCCCCCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTc                                                                                 <  2:462243/140‑1 (MQ=255)
                                                     gTTTCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCTTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAg                                                                              <  2:109149/140‑1 (MQ=255)
                                                       ttCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCa                                                                            >  1:462243/1‑140 (MQ=255)
                                                       ttCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCa                                                                            >  1:894753/1‑140 (MQ=255)
                                                         cGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATc                                                                          >  1:573567/1‑140 (MQ=255)
                                                                  tGGTACGAACTGGATCACCCGCCAGTTCGGCGAGGCCCCCCATGAGGGCGGAAAATAAAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCt                                                                 <  2:624864/140‑1 (MQ=255)
                                                                        gAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTgg                                                           >  1:238448/1‑140 (MQ=255)
                                                                        gAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTgg                                                           >  2:58961/1‑140 (MQ=255)
                                                                          aCTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGaa                                                         >  2:396354/1‑140 (MQ=255)
                                                                               aTCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGAt                                                    <  2:894753/140‑1 (MQ=255)
                                                                                   cTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCg                                                <  1:97112/140‑1 (MQ=255)
                                                                                      gCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGaa                                             >  1:717824/1‑140 (MQ=255)
                                                                                      gCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGaa                                             >  1:281434/1‑140 (MQ=255)
                                                                                      gCCAGTGCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCAAGCATCTGTTTCCCTCGCTGGAAGCGATGGCGGaa                                             >  2:129786/1‑140 (MQ=255)
                                                                                        cAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAg                                           >  1:1072014/1‑140 (MQ=255)
                                                                                              gtcgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGa                                                     >  1:91037/1‑124 (MQ=255)
                                                                                              gtcgAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGa                                                     <  2:91037/124‑1 (MQ=255)
                                                                                                    cccccccATGAGGGCGGAAAATAAAAGTAAACCGCCATATATTCCCGCAGCCTTGAACAGG‑CCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGAc                              <  2:942963/136‑1 (MQ=255)
                                                                                                     ccccccATGAGGGCGGTAAATACATGTTAGCCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGAc                              <  2:913822/137‑1 (MQ=255)
                                                                                                        cccATGAGAGCGGAAAATACAATTAAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGAcgcg                           <  2:195915/140‑1 (MQ=255)
                                                                                                            tGAGAGCGGTAAATACAAGTTAACCGCCATATATTCCCGTAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGgtgt                       >  1:1099549/1‑140 (MQ=255)
                                                                                                                  cGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTatatct                 >  1:966271/1‑138 (MQ=255)
                                                                                                                  cGGTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTg                 >  1:357763/1‑140 (MQ=255)
                                                                                                                    gTAAATACAAGTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGcc               <  2:717824/140‑1 (MQ=255)
                                                                                                                    gTAAATACAAGTAAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGcc               <  1:310149/140‑1 (MQ=255)
                                                                                                                              gTTAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAAtt     >  2:319256/1‑140 (MQ=255)
                                                                                                                               ttAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTc    <  1:396354/140‑1 (MQ=255)
                                                                                                                                tAACCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTcc   <  2:238448/140‑1 (MQ=255)
                                                                                                                                 aaCCGCCATATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTccc  >  2:881487/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                    
CGTGCTACATTGACGACAGAATCCCTTTATGGAGTATCCACGCGTTATGATACGTTTCGCTGTGATTGGTACGAACTGGATCACTCGCCAGTTCGTCGAGGCCGCCCATGAGAGCGGTAAATACAAGTTAACCGCCGTATATTCCCGCAGCCTTGAACAGGCCCAGCACTTCGCCAATGATTTTTCTGTCGAGCATCTGTTTACCTCGCTGGAAGCGATGGCGGAAAGCGATGCCATTGACGCGGTGTATATTGCCAGCCCGAATTCCC  >  NC_000913/3237265‑3237533

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