Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,316,914 G→C 7.4% intergenic (‑117/+56) ribC ← / ← treSA riboflavin synthase/trehalose synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,316,9140GC 7.4% 168.8 / 4.3 54intergenic (‑117/+56)ribC/treSAriboflavin synthase/trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base G (28/22);  new base C (2/2);  total (30/24)
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.14e-01

GGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCAAA  >  NC_002947/3316767‑3317057
                                                                                                                                                   |                                                                                                                                               
ggCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCgggg                                                                                                                                                >  3:42868/1‑149 (MQ=255)
   ggACGGCCT‑GACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGc                                                                                                                                               >  8:177265/1‑146 (MQ=255)
   ggACGGCCT‑GACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGc                                                                                                                                               <  7:177265/146‑1 (MQ=255)
    gACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAgg                                                                                                                                            <  7:212680/149‑1 (MQ=255)
       ggCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCCCCGCGCAggcgc                                                                                                                                         <  4:24591/149‑2 (MQ=255)
         ccTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGgcggcg                                                                                                                                       >  1:342670/1‑149 (MQ=255)
             gACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCt                                                                                                                                   >  1:149082/1‑149 (MQ=255)
                  tGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTgg                                                                                                                              <  1:107833/149‑1 (MQ=255)
                     cGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCAATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGCCGGCGTGCTTGTGGGc                                                                                                                            >  8:274017/1‑148 (MQ=255)
                        tgtACATGGAGTTGCTCCTGGGCACGGAAAAGAACACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGaaa                                                                                                                        <  3:281109/149‑1 (MQ=255)
                             aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGc                                                                                                                            >  2:236793/1‑140 (MQ=255)
                             aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGc                                                                                                                            <  1:236793/140‑1 (MQ=255)
                             aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGt                                                                                                                   >  5:149005/1‑149 (MQ=255)
                                      cTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACgg                                                                                                          <  4:42868/149‑1 (MQ=255)
                                       tCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACgg                                                                                                          >  6:208696/1‑148 (MQ=255)
                                            ggCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTCCGACGGCATGGc                                                                                                    >  3:218322/1‑149 (MQ=255)
                                                        aaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCGCGCCAGAtt                                                                                         >  2:190603/1‑148 (MQ=255)
                                                        aaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGAtt                                                                                         >  6:156903/1‑148 (MQ=255)
                                                          aaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATCCa                                                                                       >  5:6252/1‑148 (MQ=255)
                                                           aCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCaaa                                                                                     <  6:274427/149‑1 (MQ=255)
                                                             ttCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACa                                                                                   >  2:120511/1‑149 (MQ=255)
                                                             ttCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACa                                                                                   >  5:141726/1‑149 (MQ=255)
                                                                cGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGc                                                                                <  6:149005/149‑1 (MQ=255)
                                                                           cccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGtt                                                                     >  2:278066/1‑149 (MQ=255)
                                                                            ccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTg                                                                    <  2:149082/149‑1 (MQ=255)
                                                                             cGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                   <  1:278066/149‑1 (MQ=255)
                                                                                 gcgcAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGAGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGa                                                                <  5:208696/148‑1 (MQ=255)
                                                                                    gggatgttGCAAGTGCTGGCTTCGTGGCTGGCAAGCCGGCTTGGGTGTGTGGCGGCTGCGGGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATc                                                              <  3:204659/140‑1 (MQ=255)
                                                                                      gTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCa                                                          <  6:141726/149‑1 (MQ=255)
                                                                                            gCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGca                                                      >  3:74636/1‑147 (MQ=255)
                                                                                             cAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCaaa                                                                                     <  7:159475/115‑1 (MQ=255)
                                                                                             cAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCaaa                                                                                     >  8:159475/1‑115 (MQ=255)
                                                                                                 aCCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGcc                                               >  1:203386/1‑149 (MQ=255)
                                                                                                   cAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTc                                             >  6:51631/1‑149 (MQ=255)
                                                                                                          ccTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGGCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACgg                                       <  7:274017/148‑1 (MQ=255)
                                                                                                                tcgcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGAGCTTGTGCGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAgg                                >  5:319365/3‑149 (MQ=255)
                                                                                                                  gcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGtt                              <  2:342670/149‑1 (MQ=255)
                                                                                                                    aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCa                            >  4:95918/1‑149 (MQ=255)
                                                                                                                    aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCa                            <  5:156903/149‑1 (MQ=255)
                                                                                                                    aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTATTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCa                            >  7:329791/1‑149 (MQ=255)
                                                                                                                     aGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                           >  7:177224/1‑149 (MQ=255)
                                                                                                                      gCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCg                                            >  6:367669/1‑131 (MQ=255)
                                                                                                                      gCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCg                                            <  5:367669/131‑1 (MQ=255)
                                                                                                                      gCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                           >  1:210674/1‑148 (MQ=255)
                                                                                                                         ggCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCAt                        >  3:133253/1‑148 (MQ=255)
                                                                                                                           cTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAg                     >  1:35208/1‑149 (MQ=255)
                                                                                                                                gAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGcc                >  7:374646/1‑149 (MQ=255)
                                                                                                                                    acgccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTc            >  6:132089/5‑149 (MQ=255)
                                                                                                                                      gccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCt          <  5:132089/147‑1 (MQ=255)
                                                                                                                                       gcgcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCt          <  5:51631/148‑1 (MQ=255)
                                                                                                                                        cgcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCg        <  2:203386/149‑1 (MQ=255)
                                                                                                                                        cgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGCCGGCATCGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCg        >  7:108717/1‑149 (MQ=255)
                                                                                                                                         gcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAg                                                         <  6:271691/99‑1 (MQ=255)
                                                                                                                                         gcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAg                                                         >  5:271691/1‑99 (MQ=255)
                                                                                                                                              cgcgGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATCCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCaaa  <  6:6252/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCAAA  >  NC_002947/3316767‑3317057

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