Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,989,488 (A)8→9 coding (185/483 nt) yqeJ → protein YqeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,989,4801.A91.7% 35.6 / ‑2.7 12coding (177/483 nt)yqeJprotein YqeJ
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (4/7);  total (4/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

ACATGCACTGCAAATTGGGTAATTTTTAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAAC‑‑AAAAAAAATAAAACAGGAACTGTTGCGTTGAGTGGAACCTGGCAACAAGGTAACAGAGAAAG  >  NC_000913/2989409‑2989542
                                                                        |                                                               
aCATGCACTGCAAATTTGGTAATTTTTAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAAC‑‑AAAAAAAAt                                                       <  2:543212/81‑1 (MQ=255)
       cTGCAAATTGGGTAATTTTTAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACa                                                 <  1:543213/81‑1 (MQ=255)
         gCAAATTGGGTAATTTTTAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAg                                                >  2:543231/1‑80 (MQ=255)
              ttGGGTAATTTTTAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTg                                          >  2:558598/1‑81 (MQ=255)
                      tttttAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTGTTGCGTTg                                  <  2:543214/81‑1 (MQ=255)
                      tttttAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTGTTGCGTTg                                  <  2:543215/81‑1 (MQ=255)
                            aTGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTGTTGCGTTGAGTgg                             >  1:543244/1‑80 (MQ=255)
                            aTGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTGTTGCGTTGAGTgg                             <  2:543216/80‑1 (MQ=255)
                               aCCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTGTTGCGTTGAGTGGAAc                          <  1:543217/80‑1 (MQ=255)
                                ccAGGGACGAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTGTTGCGTTGAGTGGAACCt                        <  2:543218/81‑1 (MQ=255)
                                        gAGCAAATCTTACCATCGACTTTATGTATAACAAAAAAAAATAAAACAGGAACTGTTGCGTTGAGTGGAACCTGGcaacaa                <  2:543219/81‑1 (MQ=255)
                                                       tCGACTTTATGTATAACAAAAAAAAAATAAAACAGGAACTGTTGCGTTGAGTGGAACCTGGCAACAAGGTAACAGAGAAAg  >  1:543239/1‑81 (MQ=255)
                                                                        |                                                               
ACATGCACTGCAAATTGGGTAATTTTTAATGACCAGGGACGAGCAAATCTTACCATCGACTTTATGTATAAC‑‑AAAAAAAATAAAACAGGAACTGTTGCGTTGAGTGGAACCTGGCAACAAGGTAACAGAGAAAG  >  NC_000913/2989409‑2989542

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