Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,474,952 T→A 10.5% intergenic (‑732/+69) 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,9520TA10.5% 100.5 / 2.8 38intergenic (‑732/+69)algD/PP_1289GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (16/18);  new base A (2/2);  total (18/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 = 9.76e-01

GATATTTTTATGCCAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAG  >  NC_002947/1474805‑1475087
                                                                                                                                                    |                                                                                                                                       
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)
                                                                                                                                                    |                                                                                                                                       
GATATTTTTATGCCAGCATCTGAAAAATTTTTTCAGAATTCGTCAAAAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGAC‑GGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAG  >  NC_002947/1474805‑1475087

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