Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,410,422 T→C V97A (GTT→GCT)  yjfI → conserved protein YjfI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,4220TC95.0% 61.6 / ‑3.6 20V97A (GTT→GCT) yjfIconserved protein YjfI
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (9/10);  total (9/11)
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.54e-01

TCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTG  >  NC_000913/4410344‑4410498
                                                                              |                                                                            
tcttATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGctt                                                                            >  1:883875/1‑81 (MQ=255)
    attaAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTc                                                                          <  1:883851/79‑1 (MQ=255)
              cAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCg                                                              >  2:883886/1‑81 (MQ=255)
                gaagaTGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCt                                                            >  2:883889/1‑81 (MQ=255)
                   gatgatGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTaa                                                         <  2:883852/81‑1 (MQ=255)
                    atgatgCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTaaa                                                        >  2:883883/1‑81 (MQ=255)
                          ccGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAAtcttc                                                   <  2:883853/80‑1 (MQ=255)
                           cGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTtctc                                                 >  1:883881/1‑81 (MQ=255)
                             ttGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCtt                                               >  1:883879/1‑81 (MQ=255)
                                                 cAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCtt                                <  2:883854/76‑1 (MQ=255)
                                                 cAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTgaga                            >  1:883880/1‑80 (MQ=255)
                                                 cAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTg                               <  1:883856/77‑1 (MQ=255)
                                                 cAGCGTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTg                               <  2:883857/77‑1 (MQ=255)
                                                 cAGCGTAAAACAGGAAGGGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTg                               <  2:883855/77‑1 (MQ=255)
                                                     gTACAACAGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAAcc                       <  2:883858/81‑1 (MQ=255)
                                                        caacaGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTc                     <  2:883859/80‑1 (MQ=255)
                                                         aacaGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGc                   >  1:883882/1‑81 (MQ=255)
                                                           caGGAAGAGTATTACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTg                 >  2:883898/1‑81 (MQ=255)
                                                                      ttACATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGc      <  1:883860/81‑1 (MQ=255)
                                                                         cATTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCAtt   <  2:883861/81‑1 (MQ=255)
                                                                          aTTGCTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTg  >  2:883903/1‑81 (MQ=255)
                                                                              |                                                                            
TCTTATTAAGAAATCAGAAGATGATGCCGTTGTCATCGGTAGGGATCTCCAGCGTACAACAGGAAGAGTATTACATTGTTTTCGGGGCGTTATCGCTTAAATCTTCTCTTGAAGATATCCTGCTTGAGATAACCTCGCTGGTAGATAACGCATTG  >  NC_000913/4410344‑4410498

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