Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,591 C→A 9.5% A99D (GCC→GAC)  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,5910CA 9.5% 140.9 / 6.3 53A99D (GCC→GAC) petAubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base C (30/18);  new base A (3/2);  total (33/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 = 5.15e-01

GGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGGTAGGT  >  NC_002947/1504448‑1504739
                                                                                                                                                |                                                                                                                                                    
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)
                                                                                                                                           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)
                                                                                                                                                ccTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGGTAggt  >  1:157165/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTT‑GGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGGTAGGT  >  NC_002947/1504448‑1504739

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