Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,584 T→G 9.8% S97A (TCC→GCC)  petA → ubiquinol‑cytochrome c reductase iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,504,5840TG 9.8% 136.7 / 8.0 51S97A (TCC→GCC) petAubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (28/18);  new base G (3/2);  total (31/20)
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.83e-01

CCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGC  >  NC_002947/1504440‑1504725
                                                                                                                                                 |                                                                                                                                             
ccGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCaa                                                                                                                                           <  1:161279/149‑1 (MQ=255)
 cGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAg                                                                                                                                          <  2:294168/149‑1 (MQ=255)
        ggTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCgg                                                                                                                                   >  7:259064/1‑149 (MQ=255)
         gTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGAGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCgg                                                                                                                                   >  1:75964/1‑148 (MQ=255)
             aCATCGCCAAGGACGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCAAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACTGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAcag                                                                                                                              >  4:93548/1‑149 (MQ=255)
              cATCGCCAAGGTCGAGGCAGGTCAGCAGATGGTGGCCGAGTGGCGGGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGCGACCTGTCCGACCCCGGGTCCAAGGCCTCGGTAcagc                                                                                                                             <  8:239494/149‑1 (MQ=255)
              cATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCAGACCCCGAGTCCAAGGCCTCGGTAcagc                                                                                                                             >  8:294629/1‑149 (MQ=255)
                  gCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAcagccagc                                                                                                                         <  1:47725/149‑1 (MQ=255)
                      aGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTAc                                                                                                                     <  2:75964/149‑1 (MQ=255)
                         tCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGGCCTTGGACTCggg                                                                                                                                  <  1:177933/133‑2 (MQ=255)
                         tCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGGCCTTGGACTCggg                                                                                                                                  >  2:177933/1‑132 (MQ=255)
                                   gcagcaGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAgg                                                                                                        >  7:157612/1‑149 (MQ=255)
                                        agATGGTGTCCGAGTGGCGTGGCCAGCCCGTATTCATCTTGGTGCGAACGGAGGAAATCCTT‑GGCTATCTGAAAAAAATCAGGGGTTAGCTGTGCGAGCCCGAGTCCTAGGCCTCGGTAGTGCCAGCGTACGTCGACGCTGAGGTTCGc                                                                                                   <  7:88607/149‑1 (MQ=255)
                                            ggtggCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGTGGAAATCCTTCGGCAATGTGAAAAAAAGCACGGGTGACCTGTCCGCCCCCGAGTCCAAGGCCGGGGTACCGGCGGCGTACGTCGCCCCTGAGGTTcgcgcg                                                                                                >  1:89677/1‑149 (MQ=255)
                                              tggCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGAt                                                                                              >  4:131039/1‑148 (MQ=255)
                                                  cGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGc                                                                                         <  8:259064/149‑1 (MQ=255)
                                                      tggcgtggcCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCaa                                                                                           >  2:239189/1‑143 (MQ=255)
                                                      tggcgtggcCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCaa                                                                                           <  1:239189/143‑1 (MQ=255)
                                                               cAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCGCGGTCCACCCGGCGTACTGCCACCCTGCGGTTCGCGCGATCACGCCGGAAattctca                                                                             >  4:31136/1‑148 (MQ=255)
                                                                 gCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑CGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTAGATCAAGCCCGAAataatcaat                                                                           >  5:22126/1‑141 (MQ=255)
                                                                   ccGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAattctcattctc                                                                        >  8:11562/1‑149 (MQ=255)
                                                                   ccGTATTAATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAAACAGAAAAAAATCACGGGTCACCTGTCCGCCCCCGAGTCCAAGGCCACGGTACACCCAGCGTACGACGACACTGAGGTTCGCTCGATAAAGCCGGAAattaacattc                                                                          >  4:82080/1‑147 (MQ=255)
                                                                           aTCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACATGTCCGACACCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCt                                                                >  4:275266/1‑149 (MQ=255)
                                                                              gTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCtgt                                                              <  7:11562/148‑1 (MQ=255)
                                                                                       aCGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGATCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGAACGCGCCt                                                    >  1:309339/1‑149 (MQ=255)
                                                                                         ggaggaAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTgg                                                  >  5:282620/1‑149 (MQ=255)
                                                                                                 tCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAca                                                                                                                               <  2:5470/64‑1 (MQ=255)
                                                                                                 tCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAca                                                                                                                               >  1:5470/1‑64 (MQ=255)
                                                                                                  ccTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGc                                         >  6:117118/1‑149 (MQ=255)
                                                                                                   cTT‑GGCAATCTGAAAAAAATCACGTGTTACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGc                                         <  3:93548/148‑1 (MQ=255)
                                                                                                    tt‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGcc                                        >  3:141513/1‑148 (MQ=255)
                                                                                                       ggCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGccc                                       <  7:294629/147‑1 (MQ=255)
                                                                                                       ggCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCAc                                     >  2:33588/1‑149 (MQ=255)
                                                                                                              tGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTggg                                                 <  1:134150/130‑1 (MQ=255)
                                                                                                              tGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTggg                                                 >  2:134150/1‑130 (MQ=255)
                                                                                                              tGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCgg                              >  8:65915/1‑149 (MQ=255)
                                                                                                               gAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGc                             >  2:272783/1‑149 (MQ=255)
                                                                                                                aaaaaaaTCACGGGTGACCTGTCCGACCCCGAGTCCAATGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGc                             <  3:82080/148‑1 (MQ=255)
                                                                                                                aaaaaaaTCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGc                             <  8:157612/148‑1 (MQ=255)
                                                                                                                       tCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTg                     >  5:176480/1‑149 (MQ=255)
                                                                                                                       tCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTg                     >  4:304461/1‑149 (MQ=255)
                                                                                                                         aCGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTgg                    <  5:207865/148‑1 (MQ=255)
                                                                                                                            ggTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGcgcc                >  4:163095/1‑149 (MQ=255)
                                                                                                                             gTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCCTTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGcgccc               >  5:191096/1‑149 (MQ=255)
                                                                                                                              tGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGcgcccg              <  4:141513/149‑1 (MQ=255)
                                                                                                                                         aCCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTggg   >  2:56617/1‑149 (MQ=255)
                                                                                                                                          taccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGc  >  6:32005/3‑149 (MQ=255)
                                                                                                                                           cccTAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGATCTCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTc                                 <  7:82474/117‑1 (MQ=255)
                                                                                                                                           cccGAGTCAAAGGCCTCGGTACAGCCAGCGTAAGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCACGTCGGCCTGTGAAAGCACCTCGGCTGCTCGCCCACCTTc                                 >  8:82474/1‑117 (MQ=255)
                                                                                                                                            ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTg     <  2:204368/144‑1 (MQ=255)
                                                                                                                                            ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTg     >  1:204368/1‑144 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGC  >  NC_002947/1504440‑1504725

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