Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,178,156 C→T 89.2% D126N (GAT→AAT)  PP_2789 ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,178,1560CT89.2% 93.9 / 3.1 37D126N (GAT→AAT) PP_2789oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (18/15);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GTTTTCAGGGGTTTTATGCCGCAGAACTCCAGGGTGGTCTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATCCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTACACCAATGTCAGGTGTTCGGCCCAGGTAATGTCGTTCTGT  >  NC_002947/3178011‑3178301
                                                                                                                                                 |                                                                                                                                                 
gTTTTCAGGGGTTTTATGCCGCAGAACTCCAGGGTGGTCTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCgg                                                                                                                                                >  5:91159/1‑149 (MQ=255)
   ttCAGGGGTTTTATGCCGCAGAACTCCAGGGTGGTCTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGa                                                                                                                                              <  5:173850/148‑1 (MQ=255)
          gTTTTATGCCGCAGAACTCCAGGGTGGTCTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAg                                                                                                                                      <  6:99099/149‑1 (MQ=255)
                 gCCGCAGAACTCCAGGGTGGTCTTGAGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACAACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGAGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCtt                                                                                                                                 >  3:24418/1‑147 (MQ=255)
                   cGCTGACCTCCACGGTGGTCTTGCGCAATACGTGCACGCCCGGCATCCGATACACCCACCGGTAGTACCACCGCCGCCACACCAAGCACACCAGTAGGTGAGCCGTCCGGCCCATCACCACTATAATCCGGAAGGTCGACCCCTAGcc                                                                                                                              >  8:160496/1‑147 (MQ=255)
                           tCCAGGGTGGTCTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACAACGGCGGCGTGTCCAAGGTCACCAGTAGGTGAGCCGTGCCGCCCTTGAGCAGTTTATTCCGGAAGGGAGAGCCATTGCGCTACATg                                                                                                                       >  6:175828/1‑147 (MQ=255)
                             cAGGGTGGTCTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAg                                                                                                                    <  6:91159/148‑1 (MQ=255)
                                  tggtCTTGCGCATTTTGTGCTGGCCCGGCTAGGGATACTCGCACGGGTAGTACCACGGGGGCGTGTCCATGGTCACCAGTACGTGAGCCGTGCGGCCCTTGAGCAGTTTAATCGGGAAGGGCGAGCCCTTGCGGTACTTGAATGCaaaa                                                                                                              <  7:160496/149‑1 (MQ=255)
                                      cTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCgg                                                                                                          >  2:204910/1‑149 (MQ=255)
                                           gcATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAg                                                                                                                    <  6:238868/134‑1 (MQ=255)
                                           gcATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAg                                                                                                                    >  5:238868/1‑134 (MQ=255)
                                            cATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGaa                                                                                                    <  7:165523/149‑1 (MQ=255)
                                                              cATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAAcc                                                                                  <  1:204910/149‑1 (MQ=255)
                                                                 gCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCt                                                                               >  8:205276/1‑149 (MQ=255)
                                                                     tACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAg                                                                           <  2:149101/149‑1 (MQ=255)
                                                                        acccaccGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCgaaga                                                                                     >  2:61566/1‑136 (MQ=255)
                                                                        acccaccGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCgaaga                                                                                     <  1:61566/136‑1 (MQ=255)
                                                                          ccaccGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACAcgcg                                                                                             <  5:47140/126‑1 (MQ=255)
                                                                          ccaccGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACg                                                                      >  1:210976/1‑149 (MQ=255)
                                                                          ccaccGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGAAGAAACAcgcg                                                                                             >  6:47140/1‑126 (MQ=255)
                                                                                       cACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATCCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATgccgcc                                                         >  3:177508/1‑149 (MQ=255)
                                                                                         cggcggCGTGCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCcca                                                      >  6:59579/1‑149 (MQ=255)
                                                                                           gcggcgTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCccac                                                     <  7:205276/149‑1 (MQ=255)
                                                                                            cggcgTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCccacc                                                    <  2:210976/149‑1 (MQ=255)
                                                                                                 tgtCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCccac                                                     >  7:138917/1‑143 (MQ=255)
                                                                                                      aTGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATCCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTa                                          <  4:177508/149‑1 (MQ=255)
                                                                                                        ggTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTaca                                        >  2:131739/1‑149 (MQ=255)
                                                                                                         gTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTacac                                       >  7:99172/1‑149 (MQ=255)
                                                                                                            aCCAGTAGGTGAGCCGTGCGGCCCTTCATCAGTTTATTCGGGAAGGGCGAGCCGTTGCGGTACTTGAATGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTacac                                       <  5:175828/146‑1 (MQ=255)
                                                                                                               aGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTACACCAATGt                                 <  8:99172/149‑1 (MQ=255)
                                                                                                                   ggTGAGCCGTGCGGCCCTTGAGCAGTTTATCCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAAc                                                                                   <  5:38637/95‑1 (MQ=255)
                                                                                                                   ggTGAGCCGTGCGGCCCTTGAGCAGTTTATCCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGCTCGAAGAAc                                                                                   >  6:38637/1‑95 (MQ=255)
                                                                                                                                 cccTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCggg                                            >  2:93765/1‑120 (MQ=255)
                                                                                                                                 cccTTGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCggg                                            <  1:93765/120‑1 (MQ=255)
                                                                                                                                     tGAGCAGTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTACACCAATGTCAGGTGTTCGGCCCAGGTAATg           >  2:108250/1‑149 (MQ=255)
                                                                                                                                           gTTTATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTACACCAATCTCAGGTGTTCGGCCCAGGTAATGTCGTTc     >  2:54985/1‑149 (MQ=255)
                                                                                                                                              tATTCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTACACCAATGTCAGGTGTTCGGCCCAGGTAATGTCGTTCtgt  >  6:221711/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTTTTCAGGGGTTTTATGCCGCAGAACTCCAGGGTGGTCTTGCGCATTTGGTGCAGGCCCGGCATGCGATACACCCACCGGTAGTACCACGGCGGCGTGTCCATGGTCACCAGTAGGTGAGCCGTGCGGCCCTTGAGCAGTTTATCCGGGAAGGGCGAGCCCTTGCGGTACTTGAAGGCAAAACCGGGCAGAAACACGCGGTCGAAGAACCCCTTGAGCAACGCCGGAATGCCGCCCCACCAGATCGGGTACACCAATGTCAGGTGTTCGGCCCAGGTAATGTCGTTCTGT  >  NC_002947/3178011‑3178301

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: