Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,039,555 A→G L26L (TTA→TTG msrP → protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,039,5550AG83.3% 58.1 / 1.1 31L26L (TTA→TTGmsrPprotein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (16/9);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

GCCTTTATTACTTAGCCAATATGGGTATTCAACCTATCGTGGCAGTTAATGCATAAATATGTAAAAGCAGATCTCTGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTAAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTGGATAAAGCCGATCCCGCCGCTAATGCTGGT  >  NC_000913/2039321‑2039792
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
gCCTTTATTACTTAGCCAATATGGGTATTCAACCTATCGTGGCAGTTAATGCATAAATATGTAAAAGCAGATCTCTGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATc                                                                                                                                                                                                                                <  1:393529/250‑1 (MQ=255)
    ttattaCTTAGCCAATATGGGTATTCACCCTATCGTGGCAGTTAATGCATAAATATGTAAAAGCAGATCTCTGCTACTGTTAGGGAGGACCATCCCTGCCAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGACAAAGAATCAATTTTTAAAAGACTCAAATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGCAATCACTGGGCCTCAGTGc                                                                                                                                                                                                                           >  1:418095/1‑251 (MQ=255)
              gcaaaTATGGGTATTCAACCTATCGTGGCAGTTAATGCATAAATATGTAAAAGCAGATCTCTGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCActt                                                                                                                                                                                                                 <  2:12822/248‑1 (MQ=255)
                    atGGGTATTCAACCAATCGTGGCAGTTAATGCATAAATATGTAAAAGCAGATATCTGCTACTGTTAGGGAGGAAACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGACGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCActttctttt                                                                                                                                                                                                           >  1:186039/1‑250 (MQ=255)
                                     catGGCAGTTAATGCATAAATATGTAAAAGGAGATCTCTGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCATGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAGTCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATg                                                                                                                                                                                           <  1:100676/248‑1 (MQ=255)
                                                        ttatGTACAAGCAGATCTCTGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTACTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGGTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCCCATGCCGACATGCTTAGCTGtttt                                                                                                                                                                       >  2:96037/2‑251 (MQ=255)
                                                         tatGTAAAAGCAGATCTCTGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGttt                                                                                                                                                                       >  1:20064/1‑250 (MQ=255)
                                                                         tctGCTACTGTTAGGGAGGACACTCTCTGACAGATTAACAGTAAACGGATCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTAAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCCGCGCATGCCGATCTGCTTTGCTGGTTTAAATGGacagaccgccc                                                                                                                                                      >  2:159641/1‑241 (MQ=255)
                                                                         tctGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGCTCTGCTTAGCTGGTTTAAAGGGAACGATcgccc                                                                                                                                                      >  1:287668/1‑251 (MQ=255)
                                                                                    tAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGcccgccccc                                                                                                                                            <  1:115633/250‑1 (MQ=255)
                                                                                                       cAGATTAACAGTAAACGGCCCTTGCTGGCTAACGACAAAACAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATAAGATGTTACGGCCGAGTCGGTATTCTTTATGGAGCGTCGGCAGGTGTTAAAAGCACTGGGCATCCGGGAAAATGCACTTTCTTTGCCTAACGCTGCGCATGAGGAACTAGTTCGCAGGTTTAAAGGGGAAGATCGCCCGCCCGACCCCGCCGGCAAAGCGGTGGAgg                                                                                                                         >  2:289863/1‑249 (MQ=255)
                                                                                                             aaCAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTAAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCaa                                                                                                                  <  1:159641/251‑1 (MQ=255)
                                                                                                              aCAGTAAACGTCTCTGGCTGGCTAACGACAAAAAAGTGTGAGGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATGTTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACGGCACTTTCTTTGCCTCACGCTGCGCAGGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATcgcccgcccg                                                                                                                                                 <  2:382986/219‑1 (MQ=255)
                                                                                                              aCAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCCGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATcgcccgcccg                                                                                                                                                 >  1:382986/1‑219 (MQ=255)
                                                                                                                aGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGcc                                                                                                               >  1:449693/1‑251 (MQ=255)
                                                                                                                         acacgactctcttccgtcataaaaTTGTGATGGCTTATGAAAAAGAACTAATTTTAAAAAGAATCAGATTTTACGGCTGAGCCGGTATTCTTTATGAAGCTTCGTCAGGTGTTGAAAGCACTGGGCATCAGTGCAACGGCTCTTTCTTTGCCTCACGCTGCGCATGCCGTTCTGCTTAGCTGGTTTAAAGGGAACGATCCCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTgccgcatg                                                                                                      <  2:420683/231‑4 (MQ=255)
                                                                                                                                 ggCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCaaaa                                                                                                <  1:65535/249‑1 (MQ=255)
                                                                                                                                         acaaaaacatgggATGGCTTATGAAAAAGAATCAATTTTTAAACGAATCAGATGTTACGGCCGAGTCGGTATTATTTATAAAGCGACGGAAGGTTTTGAAAGCACTGGTCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCCAGACGATCTGCTTAGCTGGATTACAAGGAACGAACGCACGCCCCCCCCCGACGCAGAAACGTTTGAGTTCAGCACGCCTGCCGCCTGGCCCAATAACCTGcc                                                                                       >  2:442259/13‑250 (MQ=255)
                                                                                                                                          cAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGCGTCGGTATTCATTATGAAACATCGTCAGGAGTTAAAAGAACTGGGACTAAGTGCAACGGCACTTTCTTAGCCTAACGCTGAGCATGCCGATCTGATTAACTGGTTTAATGGGTAATATCGCCCGCACGCCCCCGCCGGATAAGCGATGGAGTTCAGAAAGCCTGTCTCCTGGCaaacaaccctgccgc                                                                                     >  2:316301/1‑239 (MQ=255)
                                                                                                                                          aaaagaggtgtGAGGGGTTAGGAAAAAAAATCAATTTTAAAAAAACGTAGATGTTGCGGCCTGGTGGTTATTCTTTATGGGGCGTCGGCGGGTTTGGAATGCATTGGCCATCGGGGCATCGGCGCTTTCTTGGCCTCCCGCGGCGCATGCCGAGCTGTTTAGGTGTTTTAAAGGGAAAGATCCCCCGCCCGCCCCCGCCGGAAAATCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCAc                                                                                     <  2:335108/244‑1 (MQ=255)
                                                                                                                                             aaaaGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTg                                                                                                      >  1:420683/1‑231 (MQ=255)
                                                                                                                                                       tGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGATGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCACGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCcagc                                                                           >  1:93818/1‑248 (MQ=255)
                                                                                                                                                                     aaaTCAATTTTAAAAAGAAGCAGATGTTCCGGCCGAGTGGGTATTCTTTATGAAGCGTCGGCGGGTGTGAAAAGCCCTGGGCAGCAGGGAACCTGCACTTTCTTTGCCTCCCGCTGCGCATGCCGATCTGCTTAGCGGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGtt                                                          <  1:291544/250‑1 (MQ=255)
                                                                                                                                                                          actttttAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCCTCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTAAGCCAGCCTGCCCCCTGGCAAAATAACCTGCCCCTGACCCCAGCAGATAAAGTTTCCGGTTttaa                                                      >  2:137057/3‑250 (MQ=255)
                                                                                                                                                                            tttttAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCGGGCGTTGAAAGCACTGGGCATCAGTGAAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGCTCTGCTTAGCTGGTTTAAAGGGAACGATCGCACGCCCGCCCCCGCCGGAAAAGCGCTGGCGTTCAGCAAGCATGCCGCCTGGCAAAATAACCTGCCACTGCCGCCAGCAGTTCACGCTTCCAGTTataata                                                    >  1:286449/1‑250 (MQ=255)
                                                                                                                                                                                tAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGGGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTCGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCCGCAGCAAAAACTTCCGGTTAtaataaatac                                               >  1:75468/1‑247 (MQ=255)
                                                                                                                                                                                 aaaaGTATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAAAGTCGGCTGGTGTTTCAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTATATGGAACGATCCCCCGCCCGACCCCGCCGGCACAGCGCTGGAGTTCAGCAAACCTGCCGCCTGGCAAAATAACCTGACACTGACGCCCGCAGATCAAGTTTCCGTTTATAATAACTTCt                                              >  2:74571/1‑251 (MQ=255)
                                                                                                                                                                                  aaaGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGACAGCAGATAAAGTTTCCGGTTATAATAACTTCt                                              >  2:312135/1‑250 (MQ=255)
                                                                                                                                                                                                           ggttATTCTTTATGATGCGTCGGCAGGTGTTGAAACCACGGGGCTTCGTTCCACCTGCACTTCCTTCGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGACCAACCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTGGATAAAGCCGAt                    <  1:135490/248‑1 (MQ=255)
                                                                                                                                                                                                                ttCTTTATGAAGCGTCGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCCTGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCAACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTGGATAAAGCCGCTCccgc               >  1:153755/1‑251 (MQ=255)
                                                                                                                                                                                                                               cGGCAGGTGTTGAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGATATAATAACTTCTATGAATTAGGGCTGAATAAAGCCGATCCCGCCGCTAATGCTgga  >  1:348913/1‑248 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
GCCTTTATTACTTAGCCAATATGGGTATTCAACCTATCGTGGCAGTTAATGCATAAATATGTAAAAGCAGATCTCTGCTACTGTTAGGGAGGACACTCCCTGACAGATTAACAGTAAACGGCTCTTGCTGGCTAACGACAAAAAAGTGTGATGGCTTATGAAAAAGAATCAATTTTTAAAAGAATCAGATGTTACGGCCGAGTCGGTATTCTTTATGAAGCGTCGGCAGGTGTTAAAAGCACTGGGCATCAGTGCAACTGCACTTTCTTTGCCTCACGCTGCGCATGCCGATCTGCTTAGCTGGTTTAAAGGGAACGATCGCCCGCCCGCCCCCGCCGGAAAAGCGCTGGAGTTCAGCAAGCCTGCCGCCTGGCAAAATAACCTGCCACTGACGCCAGCAGATAAAGTTTCCGGTTATAATAACTTCTATGAATTCGGGCTGGATAAAGCCGATCCCGCCGCTAATGCTGGT  >  NC_000913/2039321‑2039792

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