Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,664,878 T→A 9.5% Y4* (TAT→TAA nuoI → NADH‑quinone oxidoreductase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,664,8780TA 9.5% 115.4 / 3.8 42Y4* (TAT→TAAnuoINADH‑quinone oxidoreductase subunit I
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base A (2/2);  total (23/19)
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.45e-01

CCCACGCTATGACCAGGTGATGGACTTCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCC  >  NC_002947/4664735‑4665026
                                                                                                                                                 |                                                                                                                                                      
cacACGCTATGACCATGTGATGGACTTCAGCTGGAAGTTCTGCCTGCCGCTGACGATGATCTATTTTCTGGTGACCTCTGCGATCGTGCTCTACATCACTCCAGCCGTCGCTGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCgg                                                                                                                                                   <  3:179589/147‑1 (MQ=255)
 ccACGCTATGACCAGGTGATGGACTTCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGc                                                                                                                                                  >  6:255715/1‑149 (MQ=255)
       tATGACCAGGTGATGGACTTCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATc                                                                                                                                           >  4:112893/1‑149 (MQ=255)
             cAGGTGATGGACTTCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAg                                                                                                                                     <  3:112893/149‑1 (MQ=255)
                  gATGGACTTCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGggcac                                                                                                                                <  7:132806/149‑1 (MQ=255)
                         ttCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATc                                                                                                                                           <  8:71992/131‑1 (MQ=255)
                         ttCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATc                                                                                                                                           >  7:71992/1‑131 (MQ=255)
                                   aGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCa                                                                                                              <  5:123490/149‑1 (MQ=255)
                                                 cTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATg                                                                                                   >  6:269779/1‑146 (MQ=255)
                                                 cTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATg                                                                                                   <  5:269779/146‑1 (MQ=255)
                                                   gACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCCCCCAGCTGCGCAGCCTGGCAATGgtggt                                                                                              >  8:8030/1‑147 (MQ=255)
                                                      ccTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCCCACAGTGAGGATTTGACCCATGTACAAGTATATCGGCGACATCGTAAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGtt                                                                                              >  1:249548/1‑148 (MQ=255)
                                                         gATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCa                                                                                        >  1:64166/1‑149 (MQ=255)
                                                           tCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCAc                                                                                       <  5:255715/148‑1 (MQ=255)
                                                                    tGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATg                                                                                                   <  6:39167/127‑1 (MQ=255)
                                                                    tGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATg                                                                                                   >  5:39167/1‑127 (MQ=255)
                                                                      gTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCaa                                                                            <  1:12423/148‑1 (MQ=255)
                                                                        gACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAgcg                                                                         >  1:247813/1‑149 (MQ=255)
                                                                           cGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGa                                                                      >  3:247033/1‑149 (MQ=255)
                                                                           cGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGa                                                                      >  4:277510/1‑149 (MQ=255)
                                                                             cTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTAAATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTATCCAACGGGTTCAGCCAGCGCGaca                                                                    >  6:158430/1‑149 (MQ=255)
                                                                             cTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTAAATCGGCGAAATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGCGTTCCGCAAgcgctcca                                                                    >  6:162934/1‑145 (MQ=255)
                                                                                 gATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGCCATCGTTTAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGCGTTCCGCACGCGCGACCCCCt                                                                >  4:171401/1‑149 (MQ=255)
                                                                                 gATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATgg                                                                                                  <  3:221104/115‑1 (MQ=255)
                                                                                 gATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATgg                                                                                                  >  4:221104/1‑115 (MQ=255)
                                                                                  aTCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTAAATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTg                                                               <  5:158430/149‑1 (MQ=255)
                                                                                  aTCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTAAATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTg                                                               <  5:162934/149‑1 (MQ=255)
                                                                                              aacaCGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCgaa                                                   <  6:163748/149‑1 (MQ=255)
                                                                                                  cGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAAc                                               <  3:277510/149‑1 (MQ=255)
                                                                                                     cAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACAGACGTTAAGGGCACCGGCCCCCAGCTCCGCAGCATGGCAATGGTGTTGTCCCACGGGTTCCGCAAGCCCGACACCATGAAATACCCCGAACAAccc                                             >  1:241882/1‑149 (MQ=255)
                                                                                                      aGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCgt                                           <  2:249548/149‑1 (MQ=255)
                                                                                                      aGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCgt                                           <  2:64166/149‑1 (MQ=255)
                                                                                                         cGTCGAGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACACAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACACGTGTa                                        >  7:287220/1‑149 (MQ=255)
                                                                                                                        tGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCgcgc                           <  4:247033/149‑1 (MQ=255)
                                                                                                                            gATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTAcc                       >  7:260778/1‑149 (MQ=255)
                                                                                                                              tttGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACcgc                     >  6:148229/1‑149 (MQ=255)
                                                                                                                                 gACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATAccc                                                       >  4:1709/1‑112 (MQ=255)
                                                                                                                                 gACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATAccc                                                       <  3:1709/112‑1 (MQ=255)
                                                                                                                                   cccATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCg                >  6:200729/1‑149 (MQ=255)
                                                                                                                                            aaGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCGAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTcc        >  1:64980/1‑149 (MQ=255)
                                                                                                                                              gTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTgc                                   <  6:122503/119‑1 (MQ=255)
                                                                                                                                              gTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTgc                                   >  5:122503/1‑119 (MQ=255)
                                                                                                                                                 tatCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCAccc  >  7:2498/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
CCCACGCTATGACCAGGTGATGGACTTCAGCTGGAAGTTCTGCCTGCCGCTGACCCTGATCAATTTGCTGGTGACCGCTGCGATCGTGCTCTACAACACGCCAGCCGTCGCGGCC‑‑CAGTGAGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGTGTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCC  >  NC_002947/4664735‑4665026

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: