Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,021,365 C→G 12.2% *314S (TGA→TCA)  PP_3547 ← NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,021,3650CG12.2% 113.2 / 8.1 41*314S (TGA→TCA) PP_3547NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (26/10);  new base G (3/2);  total (29/12)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGT  >  NC_002947/4021226‑4021512
                                                                                                                                           |                                                                                                                                                   
tGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGcgcccgc                                                                                                                                            >  7:240795/1‑149 (MQ=255)
 gTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGGCGGGCGCTGAGcgcccgt                                                                                                                                            >  4:94927/1‑147 (MQ=255)
  ttGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCccgccg                                                                                                                                          >  5:185896/1‑149 (MQ=255)
     aTCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGGCGGGCGCTGAGcgcccgtcaagg                                                                                                                                       <  3:94927/149‑7 (MQ=255)
       cACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAg                                                                                                                                     >  6:174361/1‑149 (MQ=255)
        aCCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGACAACCGCTGCAAGGTACGTTGGTGTTGAACAAGCCGTTCGAAATGGAAGAACTGGCAGAGCTGGTTCGCGAGCTGCTGGAACCTTGACTGGCGCTCAGCGCCCGCCGCCCAgg                                                                                                                                    >  3:57061/1‑149 (MQ=255)
          cGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTc                                                                                                                                  >  7:149333/1‑149 (MQ=255)
                 gATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAgg                                                                                                                           <  5:174361/149‑1 (MQ=255)
                   tGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCt                                                                                                                         >  5:340604/1‑149 (MQ=255)
                      gagTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGTTACGTTGGTGTTGACCAAGCCGTTCTAACTGGAAGCACTGGCAGAGCGGGTTCGCTAGCTGCTGTAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGc                                                                                                                       <  7:284869/148‑1 (MQ=255)
                          gcggcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAGcgcccgtca                                                                                                                                          >  6:258392/1‑122 (MQ=255)
                          gcggcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAGcgcccgtca                                                                                                                                          <  5:258392/125‑4 (MQ=255)
                           cggcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTgg                                                                                                                 >  6:268836/1‑149 (MQ=255)
                           cggcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTgg                                                                                                                 >  3:171440/1‑149 (MQ=255)
                           cggcgTACAGCGCCGGCCAACCGCTGCAAGGTACTTTGGTGTTGACCAAGCCGTTCCACCGGTAAGCACTGGCAGAGCGGGTTCGCTAGCGGCGGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTgg                                                                                                                 <  7:264663/149‑1 (MQ=255)
                                            cAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCGGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATaggaggatgacta                                                                                                 >  8:318180/1‑141 (MQ=255)
                                                     cAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGTAACCTTCCCGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAAGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAAc                                                                                       <  5:170098/149‑1 (MQ=255)
                                                         gTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGTCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCa                                                                                     >  8:21224/1‑147 (MQ=255)
                                                         gTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCAGAAAGCAACCa                                                                                     >  7:317216/1‑147 (MQ=255)
                                                             gTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGAt                                                                               <  6:185896/149‑1 (MQ=255)
                                                                ggTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGc                                                                            >  5:97660/1‑149 (MQ=255)
                                                                  tgtTGACCAAGCAGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCATCGCCCGCAGCCCAGGTCGATGAGGCTGCCGGTGGCTTAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGTCCt                                                                          >  7:82247/1‑149 (MQ=255)
                                                                   gtTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTc                                                                         >  5:202187/1‑149 (MQ=255)
                                                                                   gAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCATCGACCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGACCGAAAGCAACCAGAGGATGGCATCGGCCACTTCCTCGGGa                                                         >  2:288701/1‑149 (MQ=255)
                                                                                   gAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGa                                                         >  1:181914/1‑149 (MQ=255)
                                                                                   gAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGa                                                         >  4:323361/1‑149 (MQ=255)
                                                                                         gAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGcccc                                                   >  2:29986/1‑149 (MQ=255)
                                                                                         gAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGcccc                                                   >  8:5939/1‑149 (MQ=255)
                                                                                             cACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTc                                               >  5:57397/1‑149 (MQ=255)
                                                                                              aCTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCg                                              >  1:124851/1‑149 (MQ=255)
                                                                                                  gCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGcc                                           >  1:228392/1‑148 (MQ=255)
                                                                                                          ggTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCaa                                  <  6:202187/149‑1 (MQ=255)
                                                                                                                  aGCTGCTGCAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGTCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGtt                          <  7:21224/149‑1 (MQ=255)
                                                                                                                   gctgctGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTc                         >  6:152094/1‑149 (MQ=255)
                                                                                                                     tgctgGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCATAGGATGGCGTCGGCAACTTCATCGGGCCGCCCCCCTCGGCCCATGGGCAACCCAGGtgcct                       >  5:224476/1‑145 (MQ=255)
                                                                                                                      gctgGAACCTGCCCGGGCGCCCACCGCCCGCCGCCCAGGGCGATGACGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAg                      <  8:71206/149‑1 (MQ=255)
                                                                                                                          gAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTg                  >  4:135249/1‑149 (MQ=255)
                                                                                                                                 gACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAcg            <  3:135249/148‑1 (MQ=255)
                                                                                                                                 caggggCGCTCTCGGCCCGGCGCCCATGTCGATGTGGCTTCCGGTGGCATTCGTGGGTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAcg            <  6:224476/145‑1 (MQ=255)
                                                                                                                                   cGGGCGCTGAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCgg         >  7:280244/1‑149 (MQ=255)
                                                                                                                                          tCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGt  >  8:197491/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
TGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGT  >  NC_002947/4021226‑4021512

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: ‑