Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,533,452 C→G 7.3% R57P (CGT→CCT)  mrpF ← K(+)/H(+) antiporter subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,533,4520CG 7.3% 167.3 / 3.4 55R57P (CGT→CCT) mrpFK(+)/H(+) antiporter subunit F
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base G (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TGACCAGCTCCAGCCAGAACGGCAGTACCAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCGAAGATGAACAGGCTGGCGAGGACG  >  NC_002947/2533304‑2533598
                                                                                                                                                    |                                                                                                                                                  
ttaCCAGCTCAATACAGAACGGCATTACCAGGGCTTCAGTTATTCGATAACCTCACCGATCATCAGAAACATGGCCAAGGCAAACGAGACGACGAAACAAAACAGTGGTATCTGCAGGGGACCTTCGAAGTAGGTGTCACTGGCATAAc                                                                                                                                                    <  8:324359/147‑1 (MQ=255)
  aCCAGCTCCAGCCAGAACGGCAGTACCAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACgg                                                                                                                                                  <  1:72274/149‑1 (MQ=255)
   ccAGCTCCAGCCAGAACGGCAGTACCAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACCAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGa                                                                                                                                                 <  5:255491/149‑1 (MQ=255)
       cTCCAGCCAGAACGGCAGTACCAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATAcc                                                                                                                                             >  2:342343/1‑149 (MQ=255)
            gccagAACGGCAGTACCAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACccagca                                                                                                                                        <  5:351641/149‑1 (MQ=255)
                        gTACCAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGa                                                                                                                            <  2:91217/149‑1 (MQ=255)
                           ccAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGc                                                                                                                          <  6:107790/148‑1 (MQ=255)
                            cAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACccagcacca                                                                                                                                     <  3:174568/136‑1 (MQ=255)
                            cAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTAGAAGTAGGTGTCAATGGCATAACGGATACccagaacaa                                                                                                                                     >  4:174568/1‑134 (MQ=255)
                             aGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATg                                                                                                                       >  2:36940/1‑149 (MQ=255)
                              gggCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATgg                                                                                                                      >  2:369921/1‑149 (MQ=255)
                                    cAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGa                                                                                                                >  3:215631/1‑149 (MQ=255)
                                     aGTCATTCGATCACCTCCCCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGa                                                                                                                <  7:154243/148‑1 (MQ=255)
                                     aGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGAt                                                                                                               >  2:301314/1‑149 (MQ=255)
                                            cGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAg                                                                                                         <  4:83733/148‑1 (MQ=255)
                                             gATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAgg                                                                                                        <  4:215631/148‑1 (MQ=255)
                                                 acctcaccGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTc                                                                                                   <  3:9448/149‑1 (MQ=255)
                                                 acctcaccGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTc                                                                                                   >  8:157286/1‑149 (MQ=255)
                                                   ctcaccGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGTCGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCATCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTc                                                                                                   >  2:129726/1‑147 (MQ=255)
                                                   ctcaccGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTc                                                                                                   >  6:137110/1‑147 (MQ=255)
                                                      accGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAggg                                                                                              >  1:98699/1‑149 (MQ=255)
                                                                agAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCAcgcg                                                                                    >  7:206393/1‑149 (MQ=255)
                                                                    aCTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCt                                                                                >  6:251424/1‑149 (MQ=255)
                                                                      ttGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTa                                                                                                            >  7:79087/1‑119 (MQ=255)
                                                                      ttGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTa                                                                                                            <  8:79087/119‑1 (MQ=255)
                                                                         gCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCCCCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTggg                                                                            <  1:342343/148‑1 (MQ=255)
                                                                          ccAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAgg                                                                                               >  5:103521/1‑128 (MQ=255)
                                                                          ccAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAgg                                                                                               <  6:103521/128‑1 (MQ=255)
                                                                            aaGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGa                                                                        <  1:36940/149‑1 (MQ=255)
                                                                               gCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACgg                                                                     <  5:137110/149‑1 (MQ=255)
                                                                                      aGCCTACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACgg                                                              <  7:193899/149‑1 (MQ=255)
                                                                                      aGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACg                                                               >  5:305200/1‑148 (MQ=255)
                                                                                       gCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCAAGCACGCGGTCTTGGGCGGACGGGCCACgg                                                              >  4:116337/1‑148 (MQ=255)
                                                                                        ccGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGaa                                                            <  1:301314/149‑1 (MQ=255)
                                                                                            cgaAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGc                                                        <  7:157286/149‑1 (MQ=255)
                                                                                                aCCAAACAGTGCATTCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGa                                                    >  6:81889/1‑149 (MQ=255)
                                                                                                     aCAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAgc                                               >  1:21775/1‑149 (MQ=255)
                                                                                                            aaTCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCgg                                                     >  5:121444/1‑136 (MQ=255)
                                                                                                            aaTCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCgg                                                     <  6:121444/136‑1 (MQ=255)
                                                                                                              tCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCt                                       >  5:251016/1‑148 (MQ=255)
                                                                                                                agcagGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGccatg                                    <  8:206393/149‑5 (MQ=255)
                                                                                                                     ggCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGc                                 >  1:125409/1‑147 (MQ=255)
                                                                                                                         ccTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGc                           <  5:81889/149‑1 (MQ=255)
                                                                                                                          cTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCa                                     <  7:276524/138‑1 (MQ=255)
                                                                                                                          cTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCa                                     >  8:276524/1‑138 (MQ=255)
                                                                                                                           ttCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCg                          >  7:220777/1‑148 (MQ=255)
                                                                                                                              gAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCg                          >  1:360945/1‑145 (MQ=255)
                                                                                                                              gAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCGAAGa                      <  6:305200/149‑1 (MQ=255)
                                                                                                                                aGTAGGTGTCACTGGCATAAGGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATg                                   >  2:13762/1‑134 (MQ=255)
                                                                                                                                aGTAGGTGTCACTGGCATAAGGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATg                                   >  2:19512/1‑134 (MQ=255)
                                                                                                                                aGTAGGTGTCACTGGCATAAGGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATg                                   <  1:19512/134‑1 (MQ=255)
                                                                                                                                aGTAGGTGTCACTGGCATAAGGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATg                                   <  1:13762/134‑1 (MQ=255)
                                                                                                                                aGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCGAAGAt                     <  5:251424/148‑1 (MQ=255)
                                                                                                                                    ggTGTCACTGGCATAACGGATACCCAGCATCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCGAAGATGAACa                <  1:129726/149‑1 (MQ=255)
                                                                                                                                          aCTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCGAAGATGAACAGGCTgg          >  8:274729/1‑149 (MQ=255)
                                                                                                                                                  aaCGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCGAAGATGAACAGGCTGGCGAGGACg  >  8:326533/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGACCAGCTCCAGCCAGAACGGCAGTACCAGGGCTTCAGTCATTCGATCACCTCACCGCGCAGCAGAAACTTGGCCAAGGCAAACGAGCCGACGAAACCAAACAGTGCAATCAGCAGGGCACCTTCGAAGTAGGTGTCACTGGCATAACGGATACCCAGCACCAGCATGGTGAGCATGGCCAGGATGTACAGGTAGTCCAGGGCCAGCACGCGGTCTTGGGCGGACGGGCCACGGAACAGCCGGATCAGCGCCAGGCTCATGGCCAGGGCGAAGATGAACAGGCTGGCGAGGACG  >  NC_002947/2533304‑2533598

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