Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,474,953 G→A 11.0% intergenic (‑733/+68) algD ← / ← PP_1289 GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,474,9530GA11.0% 93.1 / 2.1 36intergenic (‑733/+68)algD/PP_1289GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (14/18);  new base A (2/2);  total (16/20)
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.41e-01

GATATTTTTATGCCAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATC  >  NC_002947/1474805‑1475101
                                                                                                                                                     |                                                                                                                                                    
gaTATTTTTATGCCAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTg                                                                                                                                                      >  8:158946/1‑149 (MQ=255)
 atatTTTTATGCCAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCGTGTTg                                                                                                                                                      >  2:151782/1‑148 (MQ=255)
          tGCCAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACTGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGCAACAGGGCac                                                                                                                                              <  1:96706/148‑2 (MQ=255)
            ccAGCATCTTAAAAATTTTTTCATAATTCGTCAAAAAAGGTATTAACTGTATTATTTCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGTATTAGTGGCGTTATGGCGGTTTGCAGAGTGCCCTGCAACAGGGCacgct                                                                                                                                           <  6:321097/148‑5 (MQ=255)
             cAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCg                                                                                                                                         <  2:140415/149‑1 (MQ=255)
                cATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGAt                                                                                                                                      >  5:5836/1‑149 (MQ=255)
                        aaaTTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCAcggcgg                                                                                                                              <  3:3500/149‑1 (MQ=255)
                           ttttttCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTTCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCt                                                                                                                            >  5:183660/1‑148 (MQ=255)
                              tttCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGAc                                                                                                                         <  1:160253/148‑1 (MQ=255)
                               ttCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATATAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGGTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACct                                                                                                                       >  1:80245/1‑149 (MQ=255)
                                    aaTTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCAcggcgg                                                                                                                              <  2:191022/137‑1 (MQ=255)
                                    aaTTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCAcggcgg                                                                                                                              >  1:191022/1‑137 (MQ=255)
                                       tCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCgg                                                                                                               >  2:13073/1‑149 (MQ=255)
                                        cGTCAAAAAAGGTATGCACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCgg                                                                                                               <  7:158946/148‑1 (MQ=255)
                                          tCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCtt                                                                                                            <  7:281863/149‑1 (MQ=255)
                                             aaaaaGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAg                                                                                                         >  7:237295/1‑149 (MQ=255)
                                             aaaaaGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAg                                                                                                         >  7:233026/1‑149 (MQ=255)
                                              aaaaGGTATGAACGGTCTTATGGCGTTTTGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGc                                                                                                        >  6:25230/1‑149 (MQ=255)
                                                    taggAACGGTCTGATGGCGTTTGGAGAAAAAGGGAGGTTCGAC‑GGTGTTTTGGCGTCGGATTCGTCGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGCATCACGGCGTCTGACCTCTTCGCGGCTGACGCCGCGaa                                                                                                  <  7:168446/146‑1 (MQ=255)
                                                        aaCGGTCTTATGGCGTTTTGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAagag                                                                                               <  5:25230/148‑1 (MQ=255)
                                                        aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAgg                                                                                              <  6:198829/149‑1 (MQ=255)
                                                                       tttGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCAc                                                                               >  5:132772/1‑149 (MQ=255)
                                                                             gATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCcgcg                                                                         <  6:132772/149‑1 (MQ=255)
                                                                              aTAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGTTTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCTGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCgg                                                                        <  8:11017/149‑1 (MQ=255)
                                                                                        gTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATcgcg                                                              >  2:78895/1‑149 (MQ=255)
                                                                                             aC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCCGCTTCAGCCCCGAAGAGGGCGTGACAGGCACACTCCGCGGGGTGATCGCCCAgg                                                          >  8:230260/1‑148 (MQ=255)
                                                                                              c‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGaa                                                        >  6:16741/1‑149 (MQ=255)
                                                                                                         gCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATg                                               <  7:230260/148‑1 (MQ=255)
                                                                                                           gTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTc                                            <  1:13073/149‑1 (MQ=255)
                                                                                                              ggATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTgg                                         <  5:16741/149‑1 (MQ=255)
                                                                                                                 ttCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCtt                                      >  2:234243/1‑149 (MQ=255)
                                                                                                                   cGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTg                                     <  2:80245/148‑1 (MQ=255)
                                                                                                                    gTGGCGTGTTGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTGAGCAGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATGCCGCAGGAACACCAAATGGTATGGCTTgg                                    >  7:66847/1‑148 (MQ=255)
                                                                                                                     tGGCGTGATGGCTGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGAc                                  <  8:237295/149‑1 (MQ=255)
                                                                                                                     gggctggATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGAc                                  <  8:233026/143‑1 (MQ=255)
                                                                                                                         gTGATGGCGGGTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCATATGGTCTGGCTTGGACTGCt                              <  2:294486/149‑1 (MQ=255)
                                                                                                                                     ccgcctctGCCCTGCAACAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCCCTgtag                   >  5:132368/8‑148 (MQ=255)
                                                                                                                                        cctctGCCCTGCAACAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGCCAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgtagtag                >  1:304798/5‑148 (MQ=255)
                                                                                                                                                     gCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATAGCGCAGGACCACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATc  >  5:226346/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATATTTTTATGCCAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATC  >  NC_002947/1474805‑1475101

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