Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,505,541 T→A 12.3% V218D (GTC→GAC)  petB → cytochrome b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,505,5410TA12.3% 99.4 / 3.6 39V218D (GTC→GAC) petBcytochrome b
Reads supporting (aligned to +/- strand):  ref base T (21/13);  new base A (2/3);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGG  >  NC_002947/1505397‑1505684
                                                                                                                                                |                                                                                                                                               
ttCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGAc                                                                                                                                             <  1:164166/149‑1 (MQ=255)
       tgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCaggg                                                                                                                                      <  6:312258/149‑4 (MQ=255)
              cTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                        >  1:216261/1‑139 (MQ=255)
              cTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                        <  2:216261/140‑2 (MQ=255)
                 aTCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                       >  7:53933/1‑136 (MQ=255)
                 aTCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                       <  8:53933/138‑3 (MQ=255)
                                 tGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCg                                                                                                            >  1:38940/1‑149 (MQ=255)
                                   aaCCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCt                                                                                                          >  4:217191/1‑149 (MQ=255)
                                     ccGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTg                                                                                                         >  4:172924/1‑148 (MQ=255)
                                        cttcttCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATACTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCATGATGGCGTCGACATCAAGAAGAAAAAGGAc                                                                                                                           >  3:146529/1‑127 (MQ=255)
                                           cttcGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTa                                                                                                  >  2:180761/1‑149 (MQ=255)
                                               gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaag                                                                                                                                       >  8:209597/1‑108 (MQ=255)
                                               gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaag                                                                                                                                       <  7:209597/108‑1 (MQ=255)
                                                cccTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCg                                                                                             >  8:52824/1‑149 (MQ=255)
                                                     cATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTcc                                                                                         <  1:180761/148‑1 (MQ=255)
                                                               cccTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                           <  8:324762/136‑1 (MQ=255)
                                                               cccTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                           >  7:324762/1‑136 (MQ=255)
                                                                  tGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACAcc                                                                           >  1:317984/1‑149 (MQ=255)
                                                                     cAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTc                                                                                                              <  3:290761/111‑1 (MQ=255)
                                                                     cAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTc                                                                                                              >  4:290761/1‑111 (MQ=255)
                                                                     cAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGtt                                                                        <  3:217191/149‑1 (MQ=255)
                                                                                 tcGGCCTGCTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCATGATGACGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTc                                                            >  6:349992/1‑149 (MQ=255)
                                                                                    gCCTGGTTGTATTCCACATCCTCGCCTTGCACGAAGTGGGATCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGc                                                         <  8:274148/149‑1 (MQ=255)
                                                                                                             cTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTCCCTGTTCGTGTTCt                                >  1:321426/1‑149 (MQ=255)
                                                                                                                gCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTg                               <  2:38940/147‑1 (MQ=255)
                                                                                                                  aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGc                                                                                                           <  4:330630/69‑1 (MQ=255)
                                                                                                                  aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGc                                                                                                           >  3:330630/1‑69 (MQ=255)
                                                                                                                        tGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtggtg                     >  3:239394/1‑149 (MQ=255)
                                                                                                                          ggTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGATATTGTCGGCGTGGTGGTGTTCCTGt                                         >  1:63477/1‑127 (MQ=255)
                                                                                                                          ggTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGATATTGTCGGCGTGGTGGTGTTCCTGt                                         <  2:63477/127‑1 (MQ=255)
                                                                                                                            ttCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGACGGGGTGttct                 >  8:333848/1‑149 (MQ=255)
                                                                                                                                aacaacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttcttctt             >  4:264940/1‑149 (MQ=255)
                                                                                                                                      ccTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTg        >  3:293146/1‑148 (MQ=255)
                                                                                                                                        tGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtg                        <  4:245046/130‑1 (MQ=255)
                                                                                                                                        tGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtg                        >  3:245046/1‑130 (MQ=255)
                                                                                                                                            ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtgg                                                  <  8:103642/100‑1 (MQ=255)
                                                                                                                                            ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtgg                                                  >  7:103642/1‑100 (MQ=255)
                                                                                                                                            ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATgg  <  4:293146/148‑1 (MQ=255)
                                                                                                                                               gTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGGTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATGGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttg                  >  7:54511/1‑128 (MQ=255)
                                                                                                                                               ctCGACATCAAGAAGAAAAAGCACGAAAACGGTATTCCGCTGTATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttc                  <  8:54511/128‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGG  >  NC_002947/1505397‑1505684

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: