Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,314,235 G→C 7.4% S260W (TCG→TGG)  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,2350GC 7.4% 164.3 / 4.4 54S260W (TCG→TGG) msrPmethionine sulfoxide reductase catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (28/22);  new base C (2/2);  total (30/24)
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.83e-01

CGCTGGGCAAACGACGTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAG  >  NC_002947/5314089‑5314382
                                                                                                                                                  |                                                                                                                                                   
cGCTGGGCAAACGACGTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTAATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACTATCGAc                                                                                                                                                   <  5:260265/149‑1 (MQ=255)
 tcttggCAAACGACGTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGAc                                                                                                                                                   <  5:195459/144‑1 (MQ=255)
     ggCAAACGACGTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGat                                                                                                                                                   >  6:195459/1‑143 (MQ=255)
        aaaCGACGTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGc                                                                                                                                            >  7:218404/1‑148 (MQ=255)
          acgacgTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGctct                                                                                                                                         >  7:155923/1‑149 (MQ=255)
            gacgTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCtt                                                                                                                                        >  3:126625/1‑148 (MQ=255)
                   gcgcgcCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCtt                                                                                                                                        <  4:63954/141‑1 (MQ=255)
                   gcgcgcCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCtt                                                                                                                                        >  3:63954/1‑141 (MQ=255)
                   gcgcgcCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCa                                                                                                                                 <  8:335408/148‑1 (MQ=255)
                   gcgcgcCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCa                                                                                                                                 >  7:335408/1‑148 (MQ=255)
                         cTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTc                                                                                                                           >  1:10778/1‑148 (MQ=255)
                            gCTCCAGCGGGGGTGGTTGACTTTCGGGTTCACATTTGCATATAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCaa                                                                                                                        <  8:218404/148‑1 (MQ=255)
                            gCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCAt                                                                                                                                <  6:339094/140‑1 (MQ=255)
                            gCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCAt                                                                                                                                >  5:339094/1‑140 (MQ=255)
                                cAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGcac                                                                                                                   >  1:204301/1‑149 (MQ=255)
                                    gagggTGGTCGACGGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTACCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCAGAATCGACTTGATGCTGTTGAGGCCATACGTCCACGGcgccag                                                                                                                >  4:118217/3‑147 (MQ=255)
                                       ggtggtCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGctc                                                                                                                                          >  4:64908/1‑119 (MQ=255)
                                            tCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCaa                                                                                                       <  6:236631/149‑1 (MQ=255)
                                                 tGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGgcg                                                                                                  >  4:346358/1‑149 (MQ=255)
                                                       gTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGt                                                                                            >  1:90147/1‑149 (MQ=255)
                                                           acaTTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGaa                                                                                                                                     <  5:81058/104‑1 (MQ=255)
                                                           acaTTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGaa                                                                                                                                     >  6:81058/1‑104 (MQ=255)
                                                            caTTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTgg                                                                                       >  1:263725/1‑149 (MQ=255)
                                                                ggCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTg                                                                                    <  2:204301/148‑1 (MQ=255)
                                                                 gCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTgg                                                                                   <  3:346358/148‑1 (MQ=255)
                                                                    tAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAg                                                                               <  4:126625/149‑1 (MQ=255)
                                                                       aaaCCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCCACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTc                                                                            >  3:299536/1‑149 (MQ=255)
                                                                       aaaCCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCCACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTc                                                                            >  3:304462/1‑149 (MQ=255)
                                                                        aaCCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCACCCGCAACGGCGCACCGTTCTGGTTGGGCAGTTcc                                                                           <  8:155923/149‑1 (MQ=255)
                                                                         aCCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCaaa                                                                                                            >  4:299548/1‑115 (MQ=255)
                                                                         aCCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCaaa                                                                                                            <  3:299548/115‑1 (MQ=255)
                                                                          ccATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCg                                                                         >  4:328038/1‑149 (MQ=255)
                                                                          ccATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCg                                                                         >  3:294237/1‑149 (MQ=255)
                                                                           cATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCgg                                                                        <  2:263725/149‑1 (MQ=255)
                                                                            atatTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCCACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGc                                                                       <  4:299536/149‑1 (MQ=255)
                                                                            atatTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCCACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGc                                                                       <  4:304462/149‑1 (MQ=255)
                                                                               ttCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGt                                                                    <  2:10778/149‑1 (MQ=255)
                                                                                 gaTCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTa                                                                   >  8:285017/2‑148 (MQ=255)
                                                                                 cATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTa                                                                   <  2:90147/148‑1 (MQ=255)
                                                                                  aTCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACa                                                                 >  2:313142/1‑149 (MQ=255)
                                                                                     cGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATAc                                                              >  5:21330/1‑149 (MQ=255)
                                                                                             gCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACcgc                                                       >  1:214629/1‑148 (MQ=255)
                                                                                                    ccAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAAtt                                                <  7:285017/148‑1 (MQ=255)
                                                                                                    ccAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTg                                               <  1:313142/149‑1 (MQ=255)
                                                                                                     cAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGc                                              >  6:137799/1‑149 (MQ=255)
                                                                                                                  gctgctCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATc                                 >  7:28002/1‑149 (MQ=255)
                                                                                                                         cGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTc                            <  3:328038/147‑1 (MQ=255)
                                                                                                                               cAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCAt                                  <  6:81955/135‑1 (MQ=255)
                                                                                                                               cAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCAt                                  >  5:81955/1‑135 (MQ=255)
                                                                                                                               cAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGgc                    <  6:21330/149‑1 (MQ=255)
                                                                                                                                 ggCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAg                      >  3:257726/1‑145 (MQ=255)
                                                                                                                                          gcACAATCGACTTGATGCTCTTGAAGCCACACTTGCAAGGCACCACCATTCGGATCGGCGCAGCGTTCTGGTTGGGGAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTTCATCGCCTCATCCAGGCGCAACCCTtct         <  5:317842/149‑1 (MQ=255)
                                                                                                                                              aaTCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCta      >  4:341596/1‑148 (MQ=255)
                                                                                                                                                tCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCtat     >  7:247773/1‑147 (MQ=255)
                                                                                                                                                  gACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAg  <  2:214629/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCTGGGCAAACGACGTTCGCGCGCCTGGCTCCAGCGCGGGTGGTCGACTGTCGGGTTCACATTGGCATAGAAACCATATTCATCCGGCGCCAGCCCTTCCCAGGTAGTACCGGGCTGCTCCGCTACCAGGCTGATGCGCACAATCGACTTGATGCTCTTGAAGCCATACTTCCAAGGCACCACCAAACGCAACGGCGCACCGTTCTGGTTGGGCAGTTCCCGGCCGTACATACCCACCGCGAGAATTGCCAACGGGTGCATCGCCTCATCCAGGCGCAACCCTTCTCTATAAG  >  NC_002947/5314089‑5314382

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: ‑