Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,474,996 C→A 14.0% intergenic (‑776/+25) 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,9960CA14.0% 110.8 / 7.2 42intergenic (‑776/+25)algD/PP_1289GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (13/22);  new base A (2/4);  total (16/26)
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.55e-01

AAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGA  >  NC_002947/1474851‑1475142
                                                                                                                                                 |                                                                                                                                                  
aaaaGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGc                                                                                                                                                 >  3:169308/1‑149 (MQ=255)
aaaaGGGAAGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGGGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGc                                                                                                                                                 <  6:231386/149‑1 (MQ=255)
          aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAagag                                                                                                                                        <  5:220356/148‑1 (MQ=255)
           aCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGAAGAgg                                                                                                                                       <  6:202932/148‑1 (MQ=255)
            cGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGCCCTCTTCGCGGCTGAAGCCGCGAAGAggtc                                                                                                                                     <  2:110540/149‑3 (MQ=255)
            cGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGc                                                                                                                                     <  2:231542/149‑1 (MQ=255)
                 ttATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGAc                                                                                                                                >  1:291396/1‑149 (MQ=255)
                   aTGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGa                                                                                                                                 <  2:158622/146‑1 (MQ=255)
                      gCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCt                                                                                                                           <  2:291396/149‑1 (MQ=255)
                      gCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCt                                                                                                                           <  4:169308/149‑1 (MQ=255)
                              tcaTAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACCCcgc                                                                                                                   <  7:270014/147‑1 (MQ=255)
                                 tAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCggg                                                                                                                <  1:18585/149‑1 (MQ=255)
                                 tAAAGGGAGGTGCGACGGATCTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCggg                                                                                                                <  4:216718/149‑1 (MQ=255)
                                  aaaGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCTTGCCCTGTTGCATGGCAGAGGCGGATCACTGCGGCTGACCTCTTCGCGGCTTCAGCCTCGAAGAGGGCGTGACAGGCTCACTCCGCgggg                                                                                                               <  6:149373/149‑1 (MQ=255)
                                   aaGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGt                                                                                                              >  7:1112/1‑149 (MQ=255)
                                         gacctgACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATcgc                                                                                                        <  7:94519/144‑1 (MQ=255)
                                            gCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATcgcg                                                                                                       >  3:311854/1‑147 (MQ=255)
                                             cGACGGATTTTTGGCGTCTGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAg                                                                                                    <  5:5008/149‑1 (MQ=255)
                                             cGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAg                                                                                                    <  5:85892/149‑1 (MQ=255)
                                               aCGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTAGCGGCGTGAGCCGCGAAGAGGGCGTGACACGCTAACTCCGCGGGGTGATcgc                                                                                                        >  8:94519/1‑143 (MQ=255)
                                                 ggATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAac                                                                                                <  4:311854/149‑1 (MQ=255)
                                                            gTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCa                                                                                                                                                   >  8:99857/1‑87 (MQ=255)
                                                            gTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCa                                                                                                                                                   <  7:99857/87‑1 (MQ=255)
                                                                        gCGGGATGGCGGTTTGCAGCGGGCCCTGTTGCAGGGCAGAGGCGGATCCCGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGGGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTg                                                                         <  2:262388/149‑1 (MQ=255)
                                                                                              gCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGc                                                                        >  4:345421/1‑128 (MQ=255)
                                                                                              gCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGc                                                                        <  3:345421/128‑1 (MQ=255)
                                                                                                               agGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTg                                                               >  8:91408/1‑120 (MQ=255)
                                                                                                               agGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTg                                                               <  7:91408/120‑1 (MQ=255)
                                                                                                                        aCGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgtagta                                                          >  6:327441/1‑116 (MQ=255)
                                                                                                                        aCGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgtagta                                                          <  5:327441/116‑1 (MQ=255)
                                                                                                                            gtcacgcccTCTTCGCGGCTGAAGCCGCGAAGAGGTCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTg                            <  4:280343/135‑1 (MQ=39)
                                                                                                                            gtcacgcccTCTTCGCGGCTGAAGCCGCGAAGAGGTCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTg                            >  3:280343/8‑142 (MQ=39)
                                                                                                                               cTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTg                  >  2:11831/1‑149 (MQ=255)
                                                                                                                                 gACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAgg                                                                                                   <  7:225410/66‑1 (MQ=255)
                                                                                                                                  aCCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAg                               <  6:276428/133‑1 (MQ=255)
                                                                                                                                  aCCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACACGCTCCCTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTCCTCGGCACTGTAGTAGTACCCCAGAGCAACGAACTGCTTGAGCTGTTCCGGGGTGCt                >  5:103501/1‑148 (MQ=255)
                                                                                                                                       ttCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgcg          >  7:121896/1‑149 (MQ=255)
                                                                                                                                       ttCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgc           >  8:8492/1‑148 (MQ=255)
                                                                                                                                       ttCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgcg          >  6:104553/1‑149 (MQ=255)
                                                                                                                                         cgcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATgc           <  7:8492/146‑1 (MQ=255)
                                                                                                                                          gcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctg       >  4:217641/1‑149 (MQ=255)
                                                                                                                                          gcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgct        <  8:1112/148‑1 (MQ=255)
                                                                                                                                                 cAGCCGCTAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGa  <  6:103501/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAAAGGTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGA  >  NC_002947/1474851‑1475142

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