Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,210,520 C→T 9.5% F255F (TTC→TTT PP_1953 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,210,5200CT 9.5% 114.2 / 2.3 42F255F (TTC→TTTPP_1953short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base T (2/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

AGGCTCTGTCGACACACCATTGATGCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATCGCGCGCTCGCACGAGTCGAGGCTGCCA  >  NC_002947/2210376‑2210664
                                                                                                                                                |                                                                                                                                                
aGGCTCTGTCGACACACCATTGATGCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTtcatca                                                                                                                                              <  5:73210/149‑1 (MQ=255)
    tctGTCGACACACCATTGATGCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTTATCACCg                                                                                                                                           >  6:252084/1‑148 (MQ=255)
        tctacacacCATTGATGCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGa                                                                                                                                         <  6:36982/143‑1 (MQ=255)
                 cATTGATGCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTtct                                                                                                                             >  6:358807/1‑149 (MQ=255)
                   ttGATGCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTcc                                                                                                                           <  8:318315/149‑1 (MQ=255)
                        gCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAAcg                                                                                                                      >  6:310985/1‑149 (MQ=255)
                              gCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCggg                                                                                                                 <  5:310985/148‑1 (MQ=255)
                              gCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCAAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGa                                                                                                                >  5:341910/1‑149 (MQ=255)
                               cAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGAt                                                                                                               <  2:115989/149‑1 (MQ=255)
                                aaCCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCCGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGGGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATg                                                                                                              >  1:132407/1‑149 (MQ=255)
                                   caccGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTTATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGcac                                                                                                           <  5:252084/149‑1 (MQ=255)
                                   caccGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTTATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGca                                                                                                            >  5:237384/1‑148 (MQ=255)
                                     ccGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTTATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACAt                                                                                                         <  6:237384/149‑1 (MQ=255)
                                     ccGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGa                                                                                                                                         >  6:208300/1‑117 (MQ=255)
                                     ccGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGa                                                                                                                                         <  5:208300/117‑1 (MQ=255)
                                           gTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAAtt                                                                                                   <  6:341910/149‑1 (MQ=255)
                                              gTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCAc                                                                                                <  5:358807/149‑1 (MQ=255)
                                                 cAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGTATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCt                                                                                              <  1:158376/148‑1 (MQ=255)
                                                      ccTGGTCAGCATGCAACACATCAAGCGGCTGGCCCAACGCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCGTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCtt                                                                                         <  2:132407/148‑1 (MQ=255)
                                                        tGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGcc                                                                                      >  5:361898/1‑149 (MQ=255)
                                                              gCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTtcatc                                                                                                                                               >  2:136671/1‑86 (MQ=255)
                                                                  gCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATtat                                                                               >  6:107868/1‑146 (MQ=255)
                                                                  gCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATaa                                                                             <  5:107868/148‑1 (MQ=255)
                                                                         aTCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTTTAAGAATGt                                                                       >  3:262412/1‑147 (MQ=255)
                                                                                ggCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTc                                                               >  4:353154/1‑148 (MQ=255)
                                                                                  cTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTcaac                                                            <  7:134140/149‑1 (MQ=255)
                                                                                        cAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATg                                                      >  4:58204/1‑149 (MQ=255)
                                                                                          aCCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTc                                                    <  6:361898/149‑1 (MQ=255)
                                                                                                 gagCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCaaa                                             <  2:229693/149‑1 (MQ=255)
                                                                                                         ccAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGAtt                                     >  3:328842/1‑149 (MQ=255)
                                                                                                         ccAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGAtt                                     >  1:206220/1‑149 (MQ=255)
                                                                                                          cAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAAAGGCGGGATGCACATGTAATTCACCTGTCTTGACAAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGAAGGAttt                                    >  1:56526/1‑149 (MQ=255)
                                                                                                            aaGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTtga                                  <  2:206220/149‑1 (MQ=255)
                                                                                                             aGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTtgat                                 <  3:58204/149‑1 (MQ=255)
                                                                                                                      ttCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATGAAAGCAGGATTTGATGATcgcgcg                        >  1:365885/1‑149 (MQ=255)
                                                                                                                          ttGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAg                                                                            <  8:103252/93‑1 (MQ=255)
                                                                                                                          ttGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAg                                                                            >  7:103252/1‑93 (MQ=255)
                                                                                                                               ggTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGt                                                     <  6:38320/111‑1 (MQ=255)
                                                                                                                                        tCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATCGCGCGCTCGCAc                 >  8:298518/1‑138 (MQ=255)
                                                                                                                                        tCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATCGCGCGCTCGCAc                 <  7:298518/138‑1 (MQ=255)
                                                                                                                                           tCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATc                             >  2:9461/1‑123 (MQ=255)
                                                                                                                                           tCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATc                             <  1:9461/123‑1 (MQ=255)
                                                                                                                                            ccTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATCGCGCGCTCGCACGAGTCGAGGCTGCCa  >  2:106063/1‑149 (MQ=255)
                                                                                                                                            ccTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATCGCGCGCTCGCACGAGTCGAGGCTGCCa  >  3:199464/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGGCTCTGTCGACACACCATTGATGCATGAGCAACCACCGGAAGTCGTTCAGTTCCTGGTCAGCATGCAACCCATCAAGCGGCTGGCCCAACCCGAGGAGCTTGCCCAAAGCATCCTGTTCCTTGCCGGTGAGCATTCGTCCTTCATCACCGGACAGACGCTTTCTCCCAACGGCGGGATGCACATGTAATTCACCTGTCTTGCCCAATTATAAGAATGTAAGGGTTCAACAACATGTCTATCAAAGCAGGATTTGATGATCGCGCGCTCGCACGAGTCGAGGCTGCCA  >  NC_002947/2210376‑2210664

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