Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,588 A→T 9.5% K98M (AAG→ATG) ‡ 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,5880AT 9.5% 131.8 / 5.4 53K98M (AAG→ATG) ‡petAubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base A (29/19);  new base T (3/2);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

CCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGG  >  NC_002947/1504440‑1504734
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                                           >  2:177933/1‑132 (MQ=255)
                         tCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGGCCTTGGACTCggg                                                                                                                                           <  1:177933/133‑2 (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                                                                                                                                        >  1:5470/1‑64 (MQ=255)
                                                                                                 tCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAca                                                                                                                                        <  2:5470/64‑1 (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              >  1:204368/1‑144 (MQ=255)
                                                                                                                                            ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTg              <  2:204368/144‑1 (MQ=255)
                                                                                                                                                   cAATGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATggg  <  3:31136/149‑1 (MQ=255)
                                                                                                                                                   cAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGcc          >  8:190192/1‑141 (MQ=255)
                                                                                                                                                   cAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGcc          <  7:190192/141‑1 (MQ=255)
                                                                                                                                                    aaGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATggg  >  7:96700/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGG  >  NC_002947/1504440‑1504734

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