Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,505,544 A→C 8.1% 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,5440AC 8.1% 174.9 / 3.6 59D219A (GAC→GCC) petBcytochrome b
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base C (2/3);  total (27/32)
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.00e+00

GTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAA  >  NC_002947/1505400‑1505681
                                                                                                                                                 |                                                                                                                                         
gtggtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCAt                                                                                                                                        >  3:46019/1‑148 (MQ=255)
 tggtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCa                                                                                                                                      >  7:86647/1‑149 (MQ=255)
  ggtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCATGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATc                                                                                                                                       <  6:9391/147‑1 (MQ=255)
  ggtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCATGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATc                                                                                                                                       >  5:9391/1‑147 (MQ=255)
   gtgACTACCTGATCTCGTGCATCACCCTGAACCGCTGCTTCGCCCTGGATGTGGTTGCCCTGCGAATCGTGATTCTGGGCCTTGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCGGATGGCGTCGAGTTCaag                                                                                                                                    <  8:168247/149‑1 (MQ=255)
   gtgACTACCTGATCTCGGGCATCTCCCTGAACCGCTTCTTCGCGCTGCATGTGGTTGCCCTGCCAATGGTGATTCTCGGCGTGGTGCTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGATCAACCCTGATGTCGACGCCATCagg                                                                                                                                    <  3:166110/149‑3 (MQ=255)
   gtgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaa                                                                                                                                     <  2:249965/148‑1 (MQ=255)
    tgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCaggg                                                                                                                                   <  4:46019/149‑4 (MQ=255)
    tgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCaggg                                                                                                                                   <  3:62875/149‑4 (MQ=255)
     gacacccTGATCTCGGCCATCACCCTTAACCGCATCTTCGCCCTGCATGTGGATGCCCTGCCAATGGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTAGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaagac                                                                                                                                  <  7:207401/144‑2 (MQ=255)
     gACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTAGGCCTGGTGGTATTGCACATCCTCGCATTGCAAGAAGTGTGTTCGAACAACCCTGATGGCGTCGACATCaagaa                                                                                                                                  >  8:156017/1‑149 (MQ=255)
       cTACCTGATCTCGGGCATCAACCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCATGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaagaag                                                                                                                                 >  5:251848/1‑148 (MQ=255)
        tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCGATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACAt                                                                                                                                        >  1:62340/1‑141 (MQ=255)
        tACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACAt                                                                                                                                        <  2:62340/140‑1 (MQ=255)
          ccTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGACCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGaaa                                                                                                                              <  6:114497/148‑1 (MQ=255)
             gATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCCCATCCTCGCCTTGCACGAAGTGGGTTCGCACAACCCTGATGGCGTCGCCCTCAAGAAGAAAAAg                                                                                                                           >  5:330865/1‑148 (MQ=255)
                   gggCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGAc                                                                                                                        <  8:137378/145‑1 (MQ=255)
                    ggCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAGGGACGaaaa                                                                                                                   >  2:151896/1‑149 (MQ=255)
                    ggCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACCAAGTGGGTTAGCTCAACACTGATGGCGTCGTCATCAAGAAGAAAAAGGACGaaaa                                                                                                                   >  8:57599/1‑149 (MQ=255)
                     gCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAAc                                                                                                                  >  4:90511/1‑149 (MQ=255)
                                aaCCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGc                                                                                                        <  4:303599/148‑1 (MQ=255)
                                aaCCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGc                                                                                                        <  7:57599/148‑1 (MQ=255)
                                           cGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTAtt                                                                                             <  3:90511/148‑1 (MQ=255)
                                            gCCCTGCATGTGTTTGCCCTGCCAATCGTGATTCTCGGCCTGGTTGTATTGCACATCCTCTCCTTGCACGAAGTGGTTTCGAACAACCCTGATGTCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTc                                                                                            <  7:156017/148‑1 (MQ=255)
                                             cccTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCg                                                                                          <  3:70467/149‑1 (MQ=255)
                                               cTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                        >  1:267033/1‑149 (MQ=255)
                                                    aTGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCAc                                                                                    >  2:194051/1‑149 (MQ=255)
                                                       tgGTTGCCCTGCCAATCGTGATTCTCGGCCTGCTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCg                                                                                 >  1:40167/1‑149 (MQ=255)
                                                        ggTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCg                                                                                 >  3:11216/1‑148 (MQ=255)
                                                              ccTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCCTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGAACTaca                                                                          >  2:37399/1‑149 (MQ=255)
                                                               cTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCg                                                                                                         >  2:217357/1‑117 (MQ=255)
                                                               cTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCg                                                                                                         <  1:217357/117‑1 (MQ=255)
                                                                tGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACAcc                                                                        >  4:180302/1‑149 (MQ=255)
                                                                  ccAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCg                                                                       <  4:11216/148‑1 (MQ=255)
                                                                  ccAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCg                                                                       <  2:267033/148‑1 (MQ=255)
                                                                     aTCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTaa                                                                   <  7:57362/149‑1 (MQ=255)
                                                                        gTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGa                                                                >  5:115927/1‑149 (MQ=255)
                                                                         tGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGa                                                                >  3:135004/1‑148 (MQ=255)
                                                                                   ccTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTc                                                                                                           >  4:71560/1‑95 (MQ=255)
                                                                                   ccTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTc                                                                                                           <  3:71560/95‑1 (MQ=255)
                                                                                       gtggtATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCAc                                                                                    >  4:241933/1‑114 (MQ=255)
                                                                                       gtggtATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCAc                                                                                    <  3:241933/114‑1 (MQ=255)
                                                                                             ttGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGgtgt                                            >  8:259533/1‑148 (MQ=255)
                                                                                               tcacaTCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAATGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTTGTGTTc                                          <  6:251848/147‑1 (MQ=255)
                                                                                               gCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACTAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTc                                          <  6:87932/148‑1 (MQ=255)
                                                                                               gCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTc                                          <  4:135004/148‑1 (MQ=255)
                                                                                                 acaTCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGAATTCCGCTGGATGGTAATCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGCTGGTGTTccag                                       >  1:150314/1‑147 (MQ=255)
                                                                                                  caTCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGt                                      <  2:40167/149‑1 (MQ=255)
                                                                                                    tCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTc                                    <  3:180302/149‑1 (MQ=255)
                                                                                                      cTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCgt                                  <  6:115927/149‑1 (MQ=255)
                                                                                                      cTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCg                                   <  2:28251/148‑1 (MQ=255)
                                                                                                            ttGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGtac                                                                              <  4:48451/99‑1 (MQ=255)
                                                                                                                   aaGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtg                     >  5:82781/1‑149 (MQ=255)
                                                                                                                       gggTTCTAACAAACCTGATGGCGTCGACATCAAGAAGAATAAGGTCGAAAACGGTATTCCGCTGGTTGGTTTTGCGTTCGAGCAGTACTACAGCGTTAAAGAGTTTGTCGGCGTGGTGGTGTTTCTGTTCGTGTTGTGCGCCgtggtg                  <  2:150314/148‑1 (MQ=255)
                                                                                                                          ttCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttct              >  1:175342/1‑149 (MQ=255)
                                                                                                                             gAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttcttc            <  7:259533/148‑1 (MQ=255)
                                                                                                                                      tGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGaaa  >  2:185853/1‑149 (MQ=255)
                                                                                                                                          ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGtt                <  1:259370/131‑1 (MQ=255)
                                                                                                                                          ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACAACGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGtt                >  2:259370/1‑131 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
GTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAA  >  NC_002947/1505400‑1505681

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