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.A84.2% 49.8 / 4.6 19coding (744/1293 nt)yjcFpentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (11/5);  total (13/6)
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.24e-01

GTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTG  >  NC_000913/4282098‑4282581
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         
gTCGGTACATATAAACATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCaagg                                                                                                                                                                                                                                            >  1:307457/1‑251 (MQ=255)
            aaCCATTTCGACGACAGTATTTGTTAACCGCGTAAAAAAGAGCTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAAACAGATAAGCATGTTCAGAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTTGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATACCATCTAAAAAATTTGATTATTCAGGGAAGGTAttttt                                                                                                                                                                                                                                 >  2:304996/1‑250 (MQ=255)
                         aaCGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTAtt                                                                                                                                                                                                                  <  2:28282/251‑1 (MQ=255)
                                 tGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCAttgt                                                                                                                                                                                                            <  1:132660/250‑1 (MQ=255)
                                        cgcgTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGt                                                                                                                                                                                                    <  1:138135/251‑1 (MQ=255)
                                           gTAAAAAATAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGa                                                                                                                                                                                                 <  2:266089/251‑1 (MQ=255)
                                             aaaaaaGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAAt                                                                                                                                                                                               <  1:333551/251‑1 (MQ=255)
                                                                   aaGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATg                                                                                                                                                                         >  1:180277/1‑251 (MQ=255)
                                                                                  aaggTTGCTTTATTCAGGTTGACACGGGACAGATCAGCCTGTTCACAGTTAGTTAAATCATGTTTAGAGAAGCTAAGGGCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATCTACAAAATTTGAATACTCAAGGACGGTATTTTTAATCTTGGCTCTATTCATTGTTGCATTGTAGCATTGAACATGCATAGCGGCGCATGAATTTAGAGTTCCgg                                                                                                                                                          >  2:189704/3‑250 (MQ=255)
                                                                                         cTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCAt                                                                                                                                                   >  1:33609/1‑251 (MQ=255)
                                                                                                                   tCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAACTTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCCTATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCAt                                                                                                                         >  1:204256/1‑251 (MQ=255)
                                                                                                                             cAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATAc                                                                                                                >  2:96411/1‑250 (MQ=255)
                                                                                                                              aaaGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGACTACa                                                                                                               >  1:467059/1‑250 (MQ=255)
                                                                                                                                   taattaaATCACGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCTTCGAGAATACAATTATc                                                                                                         >  2:101347/1‑251 (MQ=255)
                                                                                                                                                       aaGCTAAGGTCGGCTCTAACTAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCAAtt                                                                                      >  1:371794/1‑250 (MQ=255)
                                                                                                                                                              gggcGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATc                                                                                                                        <  2:430801/206‑1 (MQ=255)
                                                                                                                                                              ggTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATc                                                                                                                        >  1:430801/1‑209 (MQ=255)
                                                                                                                                                                                                      ataAAGATTTACCGCAGCCATGTATGCCATATAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGACATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTCAAATTCGTTTTATACATATTACAGTCCTGAAAAAt                                       >  2:133858/1‑249 (MQ=255)
                                                                                                                                                                                                                      gCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGc                      >  1:112361/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tataaaaaaTTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAAtt         >  2:220204/5‑251 (MQ=255)
                                                                                                                                                                                                                                   ataaaaaaTTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTc        <  1:141951/248‑1 (MQ=255)
                                                                                                                                                                                                                                          aaTTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTg  <  1:255496/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         
GTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTTATACATATTACAGTCCTGAAAAATAGCGTTTTGTAGTCTGCAACCTGCAAAATTCGAATTG  >  NC_000913/4282098‑4282581

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