Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,854,559 G→A 8.1% A141T (GCA→ACA)  PP_1658 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,854,5590GA 8.1% 138.4 / 2.1 50A141T (GCA→ACA) PP_1658hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/20);  new base A (2/2);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

AGCACCGAGGATTACTACAACTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGC  >  NC_002947/1854415‑1854707
                                                                                                                                                |                                                                                                                                                    
aGCACCGAGGATTACTACAACTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCa                                                                                                                                                    <  3:62318/147‑1 (MQ=255)
     cGAGGATTACTACAACTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTg                                                                                                                                                 >  4:69723/1‑145 (MQ=255)
          aTTACTACAACTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAg                                                                                                                                        <  5:80282/149‑1 (MQ=255)
              ctacAACTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAgc                                                                                                                                    >  4:188916/1‑149 (MQ=255)
              ctacAACTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAgc                                                                                                                                    >  8:277404/1‑149 (MQ=255)
                    cTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGAcc                                                                                                                              >  2:21164/1‑149 (MQ=255)
                    cTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGAcc                                                                                                                              <  2:196084/149‑1 (MQ=255)
                    cTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGAcc                                                                                                                              >  7:236851/1‑149 (MQ=255)
                        gtcgACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTgg                                                                                                                          <  3:188916/149‑1 (MQ=255)
                                   gaagGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTctgc                                                                                                               >  8:186053/1‑149 (MQ=255)
                                   gaagGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTctgc                                                                                                               <  3:69723/149‑1 (MQ=255)
                                   gaagGTCAAGCGCAGTGCCGTCAACGCGCTGATGCGCACCAAGCTCAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTctgc                                                                                                               <  8:156918/149‑1 (MQ=255)
                                       gTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTgg                                                                                                                          >  8:226115/1‑134 (MQ=255)
                                       gTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTgg                                                                                                                          <  7:226115/134‑1 (MQ=255)
                                           aGCGCATTGCCGTCAACGCGCTGATGCTGACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCGCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGATGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAAc                                                                                                       <  4:268666/149‑1 (MQ=255)
                                           aGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAAc                                                                                                       <  8:236851/149‑1 (MQ=255)
                                            gcgcGTTGACGGCAATGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCaa                                                                                                        >  7:102061/1‑147 (MQ=255)
                                            gcgcATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACa                                                                                                      <  7:277404/149‑1 (MQ=255)
                                              gcATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGc                                                                                                    <  2:108357/149‑1 (MQ=255)
                                                aTTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCACATTGAAGATCCAGGAGCGCGACCCTGGcc                                                                                                                        >  2:302769/1‑127 (MQ=255)
                                                aTTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCACATTGAAGATCCAGGAGCGCGACCCTGGcc                                                                                                                        <  1:302769/127‑1 (MQ=255)
                                                    ccGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGt                                                                                               <  7:186053/148‑1 (MQ=255)
                                                         aaCGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAgg                                                                                          >  7:258257/1‑148 (MQ=255)
                                                                  aTGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACg                                                                                 >  2:286104/1‑148 (MQ=255)
                                                                  aTGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACg                                                                                 <  8:102061/148‑1 (MQ=255)
                                                                   tGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACg                                                                                 >  1:256806/1‑147 (MQ=255)
                                                                        aCCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAgacga                                                                          <  8:258257/149‑1 (MQ=255)
                                                                        aCCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAgacga                                                                          >  7:303397/1‑149 (MQ=255)
                                                                          cAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGAcc                                                                        >  2:77754/1‑149 (MQ=255)
                                                                                    aaCGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTgg                                                                                                                          >  4:185888/1‑89 (MQ=255)
                                                                                    aaCGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTgg                                                                                                                          <  3:185888/89‑1 (MQ=255)
                                                                                     aCGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCt                                                             <  1:21164/149‑1 (MQ=255)
                                                                                            gcTGGCAGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCAGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTAGAGGAAACGGACGAAGACGAACCGTACGCGGCTTACATGa                                                      >  6:361694/1‑149 (MQ=255)
                                                                                              tGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATc                                                    >  8:278770/1‑149 (MQ=255)
                                                                                                gCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAgacgac                                                                         <  4:9892/126‑1 (MQ=255)
                                                                                                gCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAgacgac                                                                         >  3:9892/1‑126 (MQ=255)
                                                                                                      gTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGCCGACCCGTCCGCGGCTTCAAAGAGCCGtgatga                                            >  3:208138/1‑149 (MQ=255)
                                                                                                         gCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCt                                         >  5:249462/1‑149 (MQ=255)
                                                                                                          cccATGCCCGTGGTTACGACGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCTCCCCTGGAACTATTATGCTGATCAACAGCCATGTCGAGGAAACGGACGAAGCCCACCCGTACGCGGCTTCCAAGATCACTGATGATCt                                         >  6:300067/1‑148 (MQ=255)
                                                                                                           ccATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGa                                       >  1:334184/1‑149 (MQ=255)
                                                                                                            cATGCCGGTTGTTACTAGGTTATCCCCCGTGAGGCCTCATTGAAGATCCAGGAGCGCTACCCTGGCAGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGAt                                      <  5:361694/149‑1 (MQ=255)
                                                                                                            cATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGa                                                      <  1:182490/133‑1 (MQ=255)
                                                                                                            cATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGa                                                      >  2:182490/1‑133 (MQ=255)
                                                                                                                ccGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGATGATCCAGGAGCGCGACCCTTGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGAtgtg                                   <  5:300067/148‑1 (MQ=255)
                                                                                                                           gAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGaa                        >  6:209708/1‑148 (MQ=255)
                                                                                                                                tATCCCCCCTGAGGCCGCATTGAACCTCCAGGCCCGCGACCCTGGCCGTATTCCGCCGCTAAACACCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATccc                  <  4:208138/149‑1 (MQ=255)
                                                                                                                                 aTCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCg                 >  7:41801/1‑149 (MQ=255)
                                                                                                                                     cccGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCtt             >  3:94693/1‑149 (MQ=255)
                                                                                                                                          gAGGCCACATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATc                    >  8:25983/1‑137 (MQ=255)
                                                                                                                                          gAGGCCACATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATc                    <  7:25983/137‑1 (MQ=255)
                                                                                                                                                gCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGc  >  3:162953/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGCACCGAGGATTACTACAACTTCGTCGACGACAAGAAGGTCAAGCGCATTGCCGTCAACGCGCTGATGCGCACCAAGCTGAGCAACGGCGCGCTGGCCGTGGTAGCCCATGCCGGTGGTTACGAGGTTATCCCCCGTGAGGCCGCATTGAAGATCCAGGAGCGCGACCCTGGCCGTATTCTGCTGCTCAACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGC  >  NC_002947/1854415‑1854707

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