Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,314,327 C→T 7.0% A229A (GCG→GCA)  msrP ← methionine sulfoxide reductase catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,314,3270CT 7.0% 161.3 / 2.5 57A229A (GCG→GCA) msrPmethionine sulfoxide reductase catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCTTCAGCGTTTCGAAACGCACATAACGCGCACTGGAA  >  NC_002947/5314183‑5314471
                                                                                                                                                |                                                                                                                                                
cccTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCga                                                                                                                                              >  4:187243/1‑149 (MQ=255)
 ccTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACcgc                                                                                                                                                >  5:210127/1‑146 (MQ=255)
       cAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGc                                                                                                                                       <  8:35700/149‑1 (MQ=255)
       cAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGc                                                                                                                                       <  5:252078/149‑1 (MQ=255)
       cAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGc                                                                                                                                       <  6:202244/149‑1 (MQ=255)
              taCCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACggg                                                                                                                                >  3:182602/1‑149 (MQ=255)
                 cGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTg                                                                                                                              >  1:92701/1‑148 (MQ=255)
                 cGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTg                                                                                                                              >  2:308315/1‑148 (MQ=255)
                      tgctCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGc                                                                                                                        >  7:210056/1‑149 (MQ=255)
                      tgctCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGc                                                                                                                        <  6:210127/149‑1 (MQ=255)
                                        aTGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAAccc                                                                                                      >  3:2702/1‑149 (MQ=255)
                                        aTGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAAccc                                                                                                      <  2:92701/149‑1 (MQ=255)
                                           cgcACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTtc                                                                                                   >  7:171125/1‑149 (MQ=255)
                                            gcACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTtct                                                                                                  <  8:210056/149‑1 (MQ=255)
                                            gcACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTtc                                                                                                   >  6:74554/1‑148 (MQ=255)
                                                aaTCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATc                                                                                                                          >  6:125748/1‑121 (MQ=255)
                                                aaTCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATc                                                                                                                          <  5:125748/121‑1 (MQ=255)
                                                  tCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCtat                                                                                              >  7:57581/1‑147 (MQ=255)
                                                    gACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAg                                                                                           <  2:1336/148‑1 (MQ=255)
                                                              tctTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCtat                                                                                              >  7:170603/1‑135 (MQ=255)
                                                              tctTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCtat                                                                                              <  8:170603/135‑1 (MQ=255)
                                                                  gAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGTCCGTACATACCCACCGCGAGAATTGCCAACGGTTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCaa                                                                              <  5:74554/147‑1 (MQ=255)
                                                                     gCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCg                                                                          >  3:277434/1‑148 (MQ=255)
                                                                        aTACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGc                                                                      >  2:155803/1‑149 (MQ=255)
                                                                              ccAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTTCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGaa                                                                 <  3:187243/148‑1 (MQ=255)
                                                                              ccAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAAc                                                                >  1:31586/1‑149 (MQ=255)
                                                                                aaGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGt                                                              <  4:277434/149‑1 (MQ=255)
                                                                                   gCACCACCAAACGCAACGGCGCACCGTTCTGGTTCGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGc                                                           >  5:147764/1‑149 (MQ=255)
                                                                                       caccaAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCg                                                       >  4:194374/1‑149 (MQ=255)
                                                                                        accaAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCgg                                                      >  5:150909/1‑149 (MQ=255)
                                                                                         ccaAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGc                                                     <  1:308315/149‑1 (MQ=255)
                                                                                         ccaAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGc                                                     <  4:182602/149‑1 (MQ=255)
                                                                                          caAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATaa                                                                                            <  3:291483/109‑1 (MQ=255)
                                                                                          caAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATaa                                                                                            >  4:291483/1‑109 (MQ=255)
                                                                                                 aaCGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATAtgctg                                             >  6:22782/1‑149 (MQ=255)
                                                                                                  aCGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATAtgctgc                                            <  3:194374/149‑1 (MQ=255)
                                                                                                  aCGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATAtgctgc                                            <  1:155803/149‑1 (MQ=255)
                                                                                                       gcACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGTCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATc                                       <  2:236974/149‑1 (MQ=255)
                                                                                                          ccGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCtt                                    <  2:31586/149‑1 (MQ=255)
                                                                                                            gTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCTTCa                                  <  6:147764/149‑1 (MQ=255)
                                                                                                                  gTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCTTCAGCGttt                            >  5:139540/1‑149 (MQ=255)
                                                                                                                  gTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCTTCAGCGt                              <  8:171125/147‑1 (MQ=255)
                                                                                                                     gggCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCAAAGCCAGAACCTTGCCCAGGCATATGATCCTGATCCTTCCGCGTTTCg                          >  4:224925/1‑148 (MQ=255)
                                                                                                                        cAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCTTCAGCGTTTCGAAAc                      <  8:57581/149‑1 (MQ=255)
                                                                                                                             cccGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGc                                                     <  2:309536/113‑1 (MQ=255)
                                                                                                                                 gCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCAGCATATGCTGCGGATCCTTCAGCGTTTCGAAACGCACATa               >  2:199822/1‑147 (MQ=255)
                                                                                                                                   cGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCTTCAGCGTTTCGAAACGCACATAAcgc           <  6:356581/149‑1 (MQ=255)
                                                                                                                                      aCATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGAt                                        >  8:58474/1‑117 (MQ=255)
                                                                                                                                      aCATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGAt                                        <  7:58474/117‑1 (MQ=255)
                                                                                                                                          aCCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTg                                                            >  2:194151/1‑93 (MQ=255)
                                                                                                                                          aCCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTg                                                            <  1:194151/93‑1 (MQ=255)
                                                                                                                                            ccACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCAGCATATGCTGCGGATCCTTCAGCGTTTCGAAACGCACATAACGCGCACTGGaa  <  1:199822/149‑1 (MQ=255)
                                                                                                                                             cACTGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAg                                                                             >  2:228726/1‑73 (MQ=255)
                                                                                                                                             cACTGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAg                                                                             >  2:223632/1‑73 (MQ=255)
                                                                                                                                             cACTGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAg                                                                             <  1:228726/73‑1 (MQ=255)
                                                                                                                                             cACTGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAg                                                                             <  1:223632/73‑1 (MQ=255)
                                                                                                                                             cACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGAt                                                                                 <  4:89000/69‑1 (MQ=255)
                                                                                                                                             cACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGAt                                                                                 >  3:89000/1‑69 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAGGCCAGTCGATCAAGGCGAAGCCAGAACGTTGCCCCGGCATATGCTGCGGATCCTTCAGCGTTTCGAAACGCACATAACGCGCACTGGAA  >  NC_002947/5314183‑5314471

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