Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,201,027 A→G R154R (AGA→AGG)  yjaH → DUF1481 domain‑containing protein YjaH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,201,0270AG94.7% 125.2 / NA 39R154R (AGA→AGG) yjaHDUF1481 domain‑containing protein YjaH
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (20/16);  minor base T (0/2);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGATCTGAGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGAGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGC  >  NC_000913/4200897‑4201150
                                                                                                                                  |                                                                                                                           
gcgGTGATCTGAGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCATATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGtgctt                                                                                                                    >  2:548561/1‑139 (MQ=255)
gcgGTGATCTGAGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGtgctg                                                                                                                    >  1:663613/1‑140 (MQ=255)
      aTCTGAGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAg                                                                                                              <  1:716901/140‑1 (MQ=255)
      aTCTGAGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAg                                                                                                              <  1:320916/140‑1 (MQ=255)
          gAGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGgc                                                                                                          >  2:122270/1‑140 (MQ=255)
           aGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGgcg                                                                                                         <  2:7571/140‑1 (MQ=255)
                 aTGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCTGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCa                                                                                                   <  1:470105/140‑1 (MQ=255)
                 aTGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCa                                                                                                   >  2:760158/1‑140 (MQ=255)
                  tGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCAt                                                                                                  >  2:109835/1‑140 (MQ=255)
                  tGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCAt                                                                                                  >  1:550759/1‑140 (MQ=255)
                   gCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATg                                                                                                 <  2:2771/140‑1 (MQ=255)
                     aaCGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATgtg                                                                                               <  1:122270/140‑3 (MQ=255)
                     aaCGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCg                                                                                               <  2:550759/140‑1 (MQ=255)
                           aaTCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGaa                                                                                         >  1:430233/1‑140 (MQ=255)
                           aaTCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGaa                                                                                         <  2:663613/140‑1 (MQ=255)
                               cGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGa                                                                                     <  2:430233/140‑1 (MQ=255)
                                 aTGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCTAGGGGCGCTGGCATGCGATGGAACAGAcc                                                                                   >  1:272729/1‑140 (MQ=255)
                                   ggggAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTCACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCa                                                                                                   <  2:644953/122‑1 (MQ=255)
                                   ggggAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTCACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCa                                                                                                   >  1:644953/1‑122 (MQ=255)
                                              gcagcTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGCGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCg                                                                      <  1:657611/140‑1 (MQ=255)
                                               cagcTTTCTAACGACC‑AATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATACGAAGGAACAGACGGTGACCACCTGCGaa                                                                    >  1:237576/1‑140 (MQ=255)
                                                      cTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCaa                                                              <  2:702946/140‑1 (MQ=255)
                                                            aCCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCAtt                                                        <  1:124280/140‑1 (MQ=255)
                                                            aCCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGGGATCACCTGCGAAGGGCAAACCAtt                                                        >  1:276887/1‑140 (MQ=255)
                                                            aCCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGAGAACACCTGCGAAGGGCAAACCAtt                                                        >  1:293930/1‑140 (MQ=255)
                                                                               ttATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGc                                     >  2:675703/1‑140 (MQ=255)
                                                                               ttATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACATGCGAAGCGCAGACAATTAAACCTGATTTAGATTCGc                                     >  1:268027/1‑140 (MQ=255)
                                                                                tATCGTGCTGACCAGATACGCCCGATTAGCGCTGCCTTACGTCAGGGGCGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCa                                    <  2:237576/140‑1 (MQ=255)
                                                                                       cTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGGAAAGCATTAAACATGATTTAGATTCGCAGGCGATa                             >  2:268539/1‑140 (MQ=255)
                                                                                          aTCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCTAACCATTAAACATGGTTTAGATTCGCAGGCGATAgcg                          >  1:391805/1‑140 (MQ=255)
                                                                                          aTCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGAt                              >  2:725205/1‑136 (MQ=255)
                                                                                          aTCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGAt                              <  1:725205/136‑1 (MQ=255)
                                                                                          aTCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAgcg                          >  2:566976/1‑140 (MQ=255)
                                                                                           tCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAgcgc                         <  1:760158/140‑1 (MQ=255)
                                                                                                 ccGCCAGATTAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCg                   >  2:686725/1‑140 (MQ=255)
                                                                                                         ttAGCGATGCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGcgccgc           >  1:724468/1‑140 (MQ=255)
                                                                                                                 gCCTTACGTCAGGGGAGTGTCGTGCTGCGCCAGGGTCGCCGGCATGCGGTGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCg   <  2:276887/140‑1 (MQ=255)
                                                                                                                 gCCTTACGTCAGGGGAGTGTCGTGCTGCGCCAGGGCCGCGGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCg   <  2:293930/140‑1 (MQ=255)
                                                                                                                 gCCTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCg   <  1:109835/140‑1 (MQ=255)
                                                                                                                 gCCTTACGTCAGGGGAGAGTCGGGCTGCGCCAGGGGCGCAGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCg   <  2:272729/140‑1 (MQ=255)
                                                                                                                  ccTTACGTCAGGGGAGGGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGTGCCGCCAGAGCCGc  >  1:192509/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                           
GCGGTGATCTGAGCTTTATGCAACGTGAATCCGATGGGGAAAAGCAGCAGCTTTCTAACGACCAAATCGATTTATACCGTTATCGTGCTGATCAGATCCGCCAGATTAGCGATGCCTTACGTCAGGGGAGAGTCGTGCTGCGCCAGGGGCGCTGGCATGCGATGGAACAGACCGTGACCACCTGCGAAGGGCAAACCATTAAACCTGATTTAGATTCGCAGGCGATAGCGCATATCGAGCGCCGCCAGAGCCGC  >  NC_000913/4200897‑4201150

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