Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,269,345 A→C I311L (ATC→CTC)  yhaC → uncharacterized protein YhaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,3450AC84.2% 52.4 / 3.0 19I311L (ATC→CTC) yhaCuncharacterized protein YhaC
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (9/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCTTGAACAGAATTAAATTAACAGAGATGAAAGAATTCTCATCAGAAAAAATATATGATTACATCGATATAATCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGATTTAAACATA  >  NC_000913/3269271‑3269425
                                                                          |                                                                                
aTTGCTTGAACAGAATTAAATTAACAGAGATGAAAGAATTCTCATCAGAAAAAATATATGATTACATCGATATACTCAtt                                                                             >  2:549569/1‑80 (MQ=255)
                           agaTGAAAGAATTCTCATCAGAAAAAATATATGATTACATCGATATAATCATTGAAGATTATGAAAATACCAAAGAAAATg                                                 <  2:549539/81‑1 (MQ=255)
                             aTGAAAGAATTCTCATCAGAAAAAATATATGATTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCt                                               <  1:549541/81‑1 (MQ=255)
                                 aaGAATTCTCATCAGAAAAAATATATGATTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATc                                           >  1:549560/1‑81 (MQ=255)
                                    aaTTCTCATCAGAAAAAATATATGATTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTgg                                        >  1:549569/1‑81 (MQ=255)
                                          catcaGAAAAAATATATGATTACATCGATATAATCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTcc                                  >  2:549552/1‑81 (MQ=255)
                                           atcaGAAAAAATATATGATTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGgtcgtc                                   <  2:549541/79‑1 (MQ=255)
                                                aaaaaaTATATGATTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCaaa                            <  1:549543/81‑1 (MQ=255)
                                                         atGATTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATtata                   >  1:549574/1‑81 (MQ=255)
                                                         atGATTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATtata                   >  2:549571/1‑81 (MQ=255)
                                                            aTTACATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTa                >  2:549564/1‑81 (MQ=255)
                                                               aCATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATg              >  1:549568/1‑80 (MQ=255)
                                                               aCATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATg              <  1:549544/80‑1 (MQ=255)
                                                                cATCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATgg             <  1:549545/80‑1 (MQ=255)
                                                                 aTCGATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGAt           >  2:549565/1‑81 (MQ=255)
                                                                    gATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGAttt         <  2:549545/80‑1 (MQ=255)
                                                                    gATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGAttt         >  2:549570/1‑80 (MQ=255)
                                                                    gATATACTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGATTTa        <  2:549546/81‑1 (MQ=255)
                                                                      tataATCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGATTTaaa      >  2:549572/1‑81 (MQ=255)
                                                                          aTCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGATTTAAACATa  <  1:549548/81‑1 (MQ=255)
                                                                          |                                                                                
ATTGCTTGAACAGAATTAAATTAACAGAGATGAAAGAATTCTCATCAGAAAAAATATATGATTACATCGATATAATCATTGAAGATTATGAAAATACCAAAGAAAATGCTTATCTGGTCGTCCCCCAAATTAATTATACTATGGATTTAAACATA  >  NC_000913/3269271‑3269425

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