Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,052,693 G→C 8.2% H139Q (CAC→CAG)  PP_3572 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,052,6930GC 8.2% 153.4 / 3.4 49H139Q (CAC→CAG) PP_3572hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base C (2/2);  total (25/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.13e-01

GGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGG  >  NC_002947/4052548‑4052828
                                                                                                                                                  |                                                                                                                                       
ggTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTgg                                                                                                                                      >  7:278123/1‑149 (MQ=255)
 gTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGgc                                                                                                                                     >  6:280129/1‑149 (MQ=255)
  tCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGgcc                                                                                                                                    >  7:190760/1‑149 (MQ=255)
     cTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGgccgcc                                                                                                                                 >  4:269764/1‑149 (MQ=255)
            cccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGt                                                                                                                          >  4:296466/1‑149 (MQ=255)
             ccAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGtt                                                                                                                         <  3:269764/149‑1 (MQ=255)
                        gcAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAg                                                                                                               >  5:77294/1‑148 (MQ=255)
                        gcAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAg                                                                                                               >  3:287270/1‑148 (MQ=255)
                           ggTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCAc                                                                                                            >  6:242561/1‑148 (MQ=255)
                            gTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCa                                                                                                          <  3:85858/149‑1 (MQ=255)
                            gTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCa                                                                                                          >  5:155714/1‑149 (MQ=255)
                             tCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAt                                                                                                         <  2:139046/149‑1 (MQ=255)
                             tCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCAt                                                                                                         <  8:190760/149‑1 (MQ=255)
                                accaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑CCTGgccgac                                                                                                                                 <  7:104860/122‑3 (MQ=255)
                                accaGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑CCTGgccgac                                                                                                                                 >  8:104860/1‑120 (MQ=255)
                                           ggTCAGCGCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGTCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCt                                                                                           <  4:216096/149‑1 (MQ=255)
                                            gTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTg                                                                                          <  3:296466/149‑1 (MQ=255)
                                            gTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTg                                                                                          >  3:112499/1‑149 (MQ=255)
                                               aGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTggg                                                                                        <  8:278123/148‑1 (MQ=255)
                                                    gCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAg                                                                                   <  6:77294/148‑1 (MQ=255)
                                                      cAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGa                                                                                 >  3:89173/1‑148 (MQ=255)
                                                      cAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGa                                                                                 >  4:5151/1‑148 (MQ=255)
                                                       aGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAca                                                                               <  4:89173/149‑1 (MQ=255)
                                                         cAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcagc                                                                             <  6:155714/149‑1 (MQ=255)
                                                         cAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGAcag                                                                              >  4:96836/1‑148 (MQ=255)
                                                                     ccgccgAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAgcc                                                                 >  7:123005/1‑149 (MQ=255)
                                                                            gTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑CGGGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCgcgc                                                          <  5:12276/149‑1 (MQ=255)
                                                                             tGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCa                                                         >  6:105295/1‑149 (MQ=255)
                                                                                aCGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTg                                                      >  7:280179/1‑149 (MQ=255)
                                                                                 cGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGt                                                     >  2:279938/1‑149 (MQ=255)
                                                                                 cGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGt                                                     <  5:280129/149‑1 (MQ=255)
                                                                                  gggCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCt                                                                                           <  7:177219/110‑1 (MQ=255)
                                                                                  gggCCCGCCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCt                                                                                           >  8:177219/1‑110 (MQ=255)
                                                                                    gCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCg                                                  <  3:291208/149‑1 (MQ=255)
                                                                                               ccACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTg                                        >  7:190950/1‑148 (MQ=255)
                                                                                                cACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGc                                                                      >  7:228319/1‑117 (MQ=255)
                                                                                                cACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGc                                                                      <  8:228319/117‑1 (MQ=255)
                                                                                                 aCCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTCTTCGCCGGCCATTGCAt                                     >  2:6426/1‑149 (MQ=255)
                                                                                                  cccGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATg                                    <  4:287270/149‑1 (MQ=255)
                                                                                                    cGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGgt                                  <  8:263631/149‑1 (MQ=255)
                                                                                                           aTCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGcc                           <  8:123005/149‑1 (MQ=255)
                                                                                                              gccgcGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAg                        <  5:105295/149‑1 (MQ=255)
                                                                                                                cgcgTGCATGGCTGCCAGCAGGCGGTCGGCCATGGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGt                       >  1:34469/1‑149 (MQ=255)
                                                                                                                    tGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTc                  <  1:6426/149‑1 (MQ=255)
                                                                                                                      cATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTccc                >  7:298644/1‑149 (MQ=255)
                                                                                                                           cTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTcccgcccg           <  3:5151/149‑1 (MQ=255)
                                                                                                                                    ggCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAgg  <  3:96836/149‑1 (MQ=255)
                                                                                                                                     gccggcGGCCA‑CCTGGCCGACCGCCTGCTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTg        >  2:238344/6‑142 (MQ=255)
                                                                                                                                     gccggcGGCCA‑CCTGGCCGACCGCCTGCTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTg        <  1:238344/137‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
GGTCGCTGTCGGCCCAATCGACGCGCAGGTCCACCAGCGGCCAGGTCAGCTCGCCAGCAATCTTCAGCGCCGCCGAGTGCACGGGCCCGGCTTCGCCACCCGCTGCCATCGCCGCGTGCATGGCTGCCAGCAGGCGGTCGGCCA‑GGTGGCCGCCGGCCTGTTCAAAAGCCAGCACCATGGCTTCGATCACCTGGGTCGAGGACAGCAGGTTGCCAGCCGCCGCGCATTGTTCGCCGGCCATTGCATGGTGTGTGCCCAGTGCTTCCCGCCCGGTGAACAGG  >  NC_002947/4052548‑4052828

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