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.A87.5% 66.8 / 2.6 24coding (744/1293 nt)yjcFpentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base A (9/12);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATATTAGAGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAG  >  NC_000913/4282088‑4282556
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                
ccATATTAGCGTCGGAACAAATAACCATTTCGAGGAACGCATTTGTTAACCGCGTAAAAAAGAGTTTGTTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAGTTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGATATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑ATAA‑TTTg                                                                                                                                                                                                                              >  2:151911/1‑250 (MQ=255)
 cATATTAGAGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTAAGCCATGTATGCCATAT‑AA‑A‑AAATTTGa                                                                                                                                                                                                                             >  2:91786/1‑250 (MQ=255)
                       aCCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTGACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTaa                                                                                                                                                                                                      <  1:217594/251‑1 (MQ=255)
                            ttCGAGGAAAGTATTTGTTAACCGCGTAAAAAAGATTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAAATCAGCATGTTCAAAATTAATTAAATCAAGTTTAGTGAAGCTAAGGTCTGCTCTAAATAAATTAAGATATGGCGCTATAACTTTATAAAGATTTACCTCAGCCATTTATGCCATAT‑AA‑A‑AAATTTGAATAATCATGGAAGGCATTTTTAatagt                                                                                                                                                                                                  <  2:259484/250‑3 (MQ=255)
                                     cGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTGTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGGCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTc                                                                                                                                                                                        <  1:42897/251‑1 (MQ=255)
                                          ttGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCAttg                                                                                                                                                                                    <  1:91786/250‑1 (MQ=255)
                                                     gTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGt                                                                                                                                                                           >  1:310373/1‑248 (MQ=255)
                                                       aaaaaaGAGTTTGCTATCAACTAAGTTTATATTCTGGATGGTGGCTTGATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAACTAATGTCGGCTCTAAATAAATTAATATATGGCGATATTACTTTATAAAGATTTACCCCAGCCATGTACGCCATAT‑AA‑A‑AAATTGGAAGAATCAATGAAGGTAGTTGTAATATTGGCTATATTCATTGGTGCATTATAGAAt                                                                                                                                                                      <  2:264350/251‑1 (MQ=255)
                                                         aaaaGAGTTTGCTATCAATTAAGTTTATATACTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTg                                                                                                                                                                    <  2:287920/251‑1 (MQ=255)
                                                           aaGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGaa                                                                                                                                                                  <  1:40263/251‑1 (MQ=255)
                                                                                 ttATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTTAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAAt                                                                                                                                             <  1:5082/250‑1 (MQ=255)
                                                                                          ggAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGACTTTAGAGTTcc                                                                                                                                   >  2:74101/1‑251 (MQ=255)
                                                                                                                                                                gaAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCAt                                                                                                                                                                                      >  1:190487/1‑130 (MQ=255)
                                                                                                                                                                gaAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCAt                                                                                                                                                                                      <  2:190487/130‑1 (MQ=255)
                                                                                                                                                                gaAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCAtt                                                             >  1:214032/1‑251 (MQ=255)
                                                                                                                                                                                ctAAATAAATTAATATATGGCGCTATTACGTTATAAAGATTTACCTAAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATAAAGGAAGGTATTGTTAATATTGGCTCTATTCATTGTTGCATTGTAGCATTGAACATGCATAGCGGAACATGAATTTAGAGTTCAGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGACTACAATTATCATAAAGTATTTTTTCCATTATGGAGCATTTAAAAt                                             >  2:150129/1‑251 (MQ=255)
                                                                                                                                                                                 taaataaaTTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAACCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAAAATGAATTTAGAGTTCCGACAGTAATTTTTATCTGAGCGAAATTTGAGTCATCAAGAATACAATTTTCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAAtt                                            >  2:369045/1‑251 (MQ=255)
                                                                                                                                                                                  aaataaatTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTaaaa                                              <  1:313590/248‑1 (MQ=255)
                                                                                                                                                                                        atTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAACTTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGtttt                                      >  1:300034/1‑250 (MQ=255)
                                                                                                                                                                                        atTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTta                                     <  1:176438/251‑1 (MQ=255)
                                                                                                                                                                                             tatatGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGACATAT‑AACA‑AA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGCAACATGAATTTAGAGTTCCGTCAGTCCTTTTTATCTGAGCGAAATTTGAGTCATAGAGAATACCATTATCAAAAAGTATTTTTTCCATTATGGCAAAATTAAACTTCGTTTtataaat                                >  2:23531/1‑248 (MQ=255)
                                                                                                                                                                                                                ataAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAAt              >  1:6148/1‑250 (MQ=255)
                                                                                                                                                                                                                  aaaGATTTAACTCAACCATGTATGCCATATCAA‑A‑AA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCGTTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTACTGAAAAATAGc           >  2:172756/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccTCAGCCATGTATGCCATAT‑AA‑A‑AAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGAAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAg  <  1:412602/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                
CCATATTAGAGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATAT‑AA‑A‑AA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAG  >  NC_000913/4282088‑4282556

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