Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,282,332 (A)5→6 coding (744/1293 nt) yjcF ← pentapeptide repeat‑containing protein YjcF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,282,3321.A94.1% 46.3 / ‑3.2 17coding (744/1293 nt)yjcFpentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (5/11);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑A‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAA  >  NC_000913/4282124‑4282541
                                                                                                                                                                                                                  |                                                                                                                                                                                                                 
aCGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTggn                                                                                                                                                                                <  1:146372/245‑2 (MQ=255)
        gTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTAATATTCTGGAGGGTTGCTTTATTCAGAATGACACGAGACAGACCAGAATGTACACCGTTAATTTAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTACATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATACTCAAGGAAGGTATTTTTAATATTGGCTCTATTCAttgttg                                                                                                                                                                  >  1:447148/1‑251 (MQ=255)
                           gTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATGTACCTCAGCCATGTATGCCATATAAAA‑AATTGGAATAATCaaggaagn                                                                                                                                                                                                <  1:150834/202‑2 (MQ=255)
                           gTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGATAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATGACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTAGTTTTAATAtt                                                                                                                                                                                   <  1:448364/215‑1 (MQ=255)
                                 tATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTGGc                                                                                                                                                                                <  1:79483/212‑1 (MQ=255)
                                                 atTCGGGAGGGTTGCTGTATTCAGGGTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGATTAGAGAAGATAATGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTCTTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTag                                                                                                                          <  1:34989/250‑1 (MQ=255)
                                                      ggAGGGTTGCTTTAGTCAGGTTGACACGAGACAGATCAGCATGTGCAAAGTTAATTAAATCAAGTGTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACATTATAAAGATTTACCTCAGGCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTGTTAATATTGGCTCTATTCAGTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTc                                                                                                                     <  1:382557/250‑1 (MQ=255)
                                                       gAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGt                                                                                                                                                            <  1:86912/210‑1 (MQ=255)
                                                                         gTTGACACGAGACAAATCAGCATGTTAAAAGGTAATTAAATCAAGGTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGGTATTACTTTATAAAGATTTACCTCAGCCATGTGTGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTAttttt                                                                                                                                                                                         <  1:508215/163‑1 (MQ=255)
                                                                            gACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTGGctc                                                                                                                                                                              >  1:138638/1‑171 (MQ=255)
                                                                                   gaCAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCaa                                                                                                                                                                                                      >  1:227224/1‑140 (MQ=255)
                                                                                                                               gCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTAtn                                                                                                    <  1:148244/194‑2 (MQ=255)
                                                                                                                                  aaGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCTATTATgg                                          <  1:134638/249‑1 (MQ=255)
                                                                                                                                              aaataaatTAATACATGGCGCTATTCCTTTATCCAGATTTCCCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAAtt                                                                                                                                                      >  1:155728/1‑129 (MQ=255)
                                                                                                                                               aataaatTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCCTATAACCCA‑TTTGACTACTCAAGGAAGGTATTTTTAATATTGGCTCTCTTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTAt                                                                                                     >  1:282636/1‑177 (MQ=255)
                                                                                                                                                                          ttATAAAGATTTACCTCAGCCATGTATTCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATAGAACTTGCAGATGGGATCATGAAGTTAGCGTTCCGTCAGTCATTTTTATCGGAGCGAAATTTTGGTTAACGAGAATACAATTATCAAAAAGTATTTTTGCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGaaa  <  1:497970/249‑1 (MQ=255)
                                                                                                                                                                                 naTTTACCTCAGCCATGTATGCCATATAAAA‑AATTTGAATAATCAAGGAAGGTATTTTTAATGTTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCaa                                                               >  1:304835/2‑181 (MQ=255)
                                                                                                                                                                                                          taaaaa‑a‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTg     <  1:331812/211‑1 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                 
ACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAA‑A‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAA  >  NC_000913/4282124‑4282541

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