Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,505,547 G→T 11.6% G217V (GGC→GTC)  petB → cytochrome b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,505,5470GT11.6% 131.7 / 8.7 52G217V (GGC→GTC) petBcytochrome b
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base T (2/4);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

CCAGTGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGG  >  NC_002947/1505398‑1505693
                                                                                                                                                     |                                                                                                                                                   
ccAGTGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACCTCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATgg                                                                                                                                                     <  1:1152784/150‑1 (MQ=255)
 cAGTGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGc                                                                                                                                                    >  1:207866/1‑150 (MQ=255)
   gTGGATTCGTGGTGACTACCTGATCTCGGGCATCACACTGAACCGCTTCTTCGCCATGCATGTGGTTGCGCTGCCAATCGTGATTATCGGCATGGTGGTATTGAACATCCTCGCCTTGCACGAAGTGGGTTCGAACAAACATGATGGGGt                                                                                                                                                  >  2:436842/1‑150 (MQ=255)
    tGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTc                                                                                                                                                 <  1:348164/150‑1 (MQ=255)
        ttCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACa                                                                                                                                             >  2:727309/1‑150 (MQ=255)
          cGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATc                                                                                                                                           <  1:727309/150‑1 (MQ=255)
            tggtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                         <  2:800028/150‑2 (MQ=255)
             ggtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                        <  1:365259/150‑3 (MQ=255)
                   tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaa                                                                                                                                 >  1:77014/1‑150 (MQ=255)
                   tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaa                                                                                                                                 >  2:500212/1‑150 (MQ=255)
                   tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTATCGGTCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaa                                                                                                                                 >  1:572050/1‑150 (MQ=255)
                       tGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAgg                                                                                                                             >  1:676224/1‑150 (MQ=255)
                        gATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGa                                                                                                                            <  2:207866/150‑1 (MQ=255)
                                 cATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGCACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACgg                                                                                                                   <  1:280466/150‑1 (MQ=255)
                                       ccTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGAAGCACACGTGAGAAAGGACGAAAACGGTATTcc                                                                                                             >  1:604919/1‑150 (MQ=255)
                                             ccGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACAAAAACGGTATTCCGCTGGa                                                                                                       <  2:676224/150‑1 (MQ=255)
                                               gCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATg                                                                                                     >  1:347067/1‑150 (MQ=255)
                                                cttcttCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                         <  2:972122/114‑2 (MQ=255)
                                                cttcttCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                         >  1:972122/1‑113 (MQ=255)
                                                       gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCt                                                                                                          <  2:9778/137‑1 (MQ=255)
                                                       gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCt                                                                                                          >  1:9778/1‑137 (MQ=255)
                                                         ccTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                           <  2:347067/150‑1 (MQ=255)
                                                         ccTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGAATGGTATTCCGtt                                                                                           >  2:510571/1‑150 (MQ=255)
                                                         ccTCCATGTGGTTGCCCTGCCAATCGTCATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTTGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                           <  1:436842/150‑1 (MQ=255)
                                                                         cTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGCTGGAGTCGACATCAAGAAGAGAGGTGCGGAGGTCATGAGTCCGCTGGATGGTATTCCGTTCCACCCGTACTACAcc                                                                           <  1:510571/150‑1 (MQ=255)
                                                                         cTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACAcc                                                                           >  2:857837/1‑150 (MQ=255)
                                                                           gCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGt                                                                         >  1:1104423/1‑150 (MQ=255)
                                                                               aTCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCCTACTACACCGTTaaa                                                                     <  1:339188/150‑1 (MQ=255)
                                                                                   tGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                        >  2:605138/1‑78 (MQ=255)
                                                                                   tGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                        <  1:605138/80‑3 (MQ=255)
                                                                                                 gtggtATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtg                                                   >  1:509573/1‑150 (MQ=255)
                                                                                                     tATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGgtgt                                               >  1:40054/1‑150 (MQ=255)
                                                                                                      aTTGCACATCCTCGCCGTGCACGAAATGGGTTCGAACAACGCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCACTGAATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGCCATTGTCGGCGTGGTGGTGtt                                              <  2:77014/150‑1 (MQ=255)
                                                                                                                     cTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTg                               >  1:1186590/1‑150 (MQ=255)
                                                                                                                         cACGAAGTGTGTTCGAACAACCCGGATGGCGTCGACATAAAGAATAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGcc                           <  2:509573/150‑1 (MQ=255)
                                                                                                                         cACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTa                                                                                                                 <  1:511713/64‑1 (MQ=255)
                                                                                                                         cACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTa                                                                                                                 >  2:511713/1‑64 (MQ=255)
                                                                                                                          aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCTGCGGGGGGGTGTGCCTGTTCGTGTTCTGCGACg                          >  2:197592/1‑150 (MQ=255)
                                                                                                                            gAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtg                                                      >  2:197689/1‑120 (MQ=255)
                                                                                                                            gAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtg                                                      <  1:197689/120‑1 (MQ=255)
                                                                                                                              aGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtggt                      <  1:857837/150‑1 (MQ=255)
                                                                                                                                  ggTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttc                  <  2:40054/150‑1 (MQ=255)
                                                                                                                                    ttCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAATAAAGTTATCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttctt                >  1:676895/1‑150 (MQ=255)
                                                                                                                                    ttCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttctt                >  2:1029416/1‑150 (MQ=255)
                                                                                                                                        aacaacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCgg                                                          >  2:178650/1‑104 (MQ=255)
                                                                                                                                        aacaacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCgg                                                          <  1:178650/104‑1 (MQ=255)
                                                                                                                                           aacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGAt                                                                                                      <  2:236694/57‑1 (MQ=255)
                                                                                                                                           aacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGAt                                                                                                      >  1:236694/1‑57 (MQ=255)
                                                                                                                                             cccTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGa       >  1:73053/1‑150 (MQ=255)
                                                                                                                                              ccTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGaa      >  2:655440/1‑150 (MQ=255)
                                                                                                                                                 gATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATg   <  1:197592/150‑1 (MQ=255)
                                                                                                                                                  aTGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATgg  <  1:500212/150‑1 (MQ=255)
                                                                                                                                                  aTGGCGGCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATgg  <  2:1104423/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCAGTGGATTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGG  >  NC_002947/1505398‑1505693

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: