Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,505,553 A→C 12.3% D219A (GAC→GCC)  petB → cytochrome b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,505,5530AC12.3% 116.0 / 11.5 49D219A (GAC→GCC) petBcytochrome b
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base C (2/4);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

TTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGG  >  NC_002947/1505406‑1505693
                                                                                                                                                   |                                                                                                                                             
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                                                                                                                                 >  2:500212/1‑150 (MQ=255)
           tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaa                                                                                                                                 >  1:77014/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                                                                                                                                         >  1:972122/1‑113 (MQ=255)
                                        cttcttCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                         <  2:972122/114‑2 (MQ=255)
                                               gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCt                                                                                                          >  1:9778/1‑137 (MQ=255)
                                               gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCt                                                                                                          <  2:9778/137‑1 (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)
                                                                                                                                                   |                                                                                                                                             
TTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGG  >  NC_002947/1505406‑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: