Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,040,289 C→T A271V (GCG→GTG)  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,040,2890CT96.2% 78.9 / NA 26A271V (GCG→GTG) msrPprotein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/11);  minor base A (1/0);  total (15/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 = 1.00e+00

GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACAC‑TGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTG  >  NC_000913/2040077‑2040527
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
aaaaTATCCTTATGTCAAAGGTTTGCGTCTCGAAGAACAAATGCATCCGCTCCCAC‑TGATCCCCGTAGGTGTTTAGGGCAAGGTGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGt                                                                                                                                                                                                          <  1:44597/250‑1 (MQ=255)
 aaaTATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACAC‑TGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGt                                                                                                                                                                                                          >  1:522135/1‑250 (MQ=255)
        cTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACAC‑TGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGCTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGAt                                                                                                                                                                                                  >  2:209140/1‑251 (MQ=255)
             gccGAAGGATTGCGCCTCGACAAAGCAATGCACCCGCTCACAC‑TGATGACCGTAGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCAccc                                                                                                                                                                                             <  1:308301/249‑1 (MQ=255)
                  aGGATTGCGCCCCGACGAGGCAATGCACCCGC‑CACCCCTGATGCCCGTAGTTTTTTATGGAAGGGTTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCt                                                                                                                                                                                        <  1:20133/251‑1 (MQ=255)
                              cgacgaAGCAATGCATCCGCTCACAC‑TGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCt                                                                                                                                                                             >  1:550375/1‑250 (MQ=255)
                                                         tgatgaCCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTAAGTATTAAGCTGACCCGCGAGCGTCCGCCAAACACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGGACGATTTATTGGTTCAggcgg                                                                                                                                                   >  2:226898/1‑250 (MQ=255)
                                                            tgaCTGTAGGTGGTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCATGGAAATATGGCTTTACAGGGATTAAATCGATATTCAGTCTTAAGCTGACCCGCGAGCGTCCACCAAACACCTGGAATCTGGCAGAGCCTGACGAATACGGTGTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGATCGATATAATGGATCAGGCGGCATc                                                                                                                                               >  2:305319/1‑251 (MQ=255)
                                                                    ggtttctatggagaggcgtatccGCCACAAAATGGCGCGCCGGTGTGACTGATTGTGCCGTGGAAATATGGCTTTAAAGTGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTTACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGAGTTATTGGTTCAGGCGGCATCCTCGATGt                                                                                                                                       <  1:258787/230‑1 (MQ=255)
                                                                               cAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCgccagc                                                                                                                            >  1:567235/1‑251 (MQ=255)
                                                                                           gCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGt                                                                                                                >  1:50818/1‑251 (MQ=255)
                                                                                               caAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTa                                                                                                             >  1:170716/1‑250 (MQ=255)
                                                                                                         cgcCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTAcgcc                                                                                                  <  1:542307/251‑1 (MQ=255)
                                                                                                             ggTGTGACTGATTGTGCCTTGGAAATAAGGCTTTAAAGGGATGAAATCGTTAGTCAGTATTTAGCTGACCCGCGAGCGTCCGCCAACCACCGGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGAGCACCCGCGCTGGTCACAGGCTACCGAACGAGTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTAcgccgacc                                                                                              <  2:330978/251‑4 (MQ=255)
                                                                                                                   aCTGATTGGTCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGGCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTgg                                                                                        <  1:250450/251‑1 (MQ=255)
                                                                                                                                             ttAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGAtt                                                                <  1:555536/249‑1 (MQ=255)
                                                                                                                                               aaaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCg                                                            >  1:99029/1‑251 (MQ=255)
                                                                                                                                                                   cAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTc                                        >  1:323249/1‑251 (MQ=255)
                                                                                                                                                                         tAAGCGGACCCGCTATCGTCCGCCAACCACCTGGAATATGGCAGTGCCTGTCTAATACGGTTTTTACGCCAACGTTAATCCGGATGTTGATCACCCGCGCTGGTCACAGGTTACCGAACGAGTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACTCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACa                                   <  2:99029/250‑1 (MQ=255)
                                                                                                                                                                            gCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCa                                >  1:544756/1‑250 (MQ=255)
                                                                                                                                                                               gACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGACTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATAGTGGCCTGGATTTGCGGGGGAATTTCTAAATGCGGCTGACAGCaaaa                             >  2:432140/1‑250 (MQ=255)
                                                                                                                                                                                aCCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACa                           <  1:209140/251‑1 (MQ=255)
                                                                                                                                                                                        gCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATTCGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACAt                   <  2:50818/251‑1 (MQ=255)
                                                                                                                                                                                          gTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATgg                 >  1:426073/1‑251 (MQ=255)
                                                                                                                                                                                                    caccTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTAAAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGATGTATCGTGGCCTGGATTTGAGGGAGAATTTATAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGttt       >  1:373601/1‑251 (MQ=255)
                                                                                                                                                                                                         ggAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCAACCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTg  >  1:380870/1‑251 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACAC‑TGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTG  >  NC_000913/2040077‑2040527

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