Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,077,683 G→T 9.5% R71R (CGC→CGA)  PP_2687 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,077,6830GT 9.5% 125.8 / 2.6 42R71R (CGC→CGA) PP_2687hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base T (2/2);  total (22/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 = 6.06e-01

CGCATTTCGATCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGACCTCGAGCATGCGTG  >  NC_002947/3077537‑3077828
                                                                                                                                                  |                                                                                                                                                 
cgcATTTCGATCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATgcgc                                                                                                                                                <  2:34442/150‑1 (MQ=255)
   aTTTCGATCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCg                                                                                                                                             <  1:1256066/150‑1 (MQ=255)
     ttCGATCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGTGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGgt                                                                                                                                           <  1:728095/150‑1 (MQ=255)
       cGATCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGgtgt                                                                                                                                         <  2:246041/150‑1 (MQ=255)
        gATCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTc                                                                                                                                        >  1:990256/1‑150 (MQ=255)
         aTCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCg                                                                                                                                       >  2:88019/1‑150 (MQ=255)
                   cGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTtgt                                                                                                                             <  1:681982/150‑1 (MQ=255)
                          gTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGCGTTCACTGTGGCGACg                                                                                                                      >  2:1036487/1‑150 (MQ=255)
                                ccACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTa                                                                                                                >  2:810697/1‑150 (MQ=255)
                                ccACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTa                                                                                                                <  2:1215482/150‑1 (MQ=255)
                                ccACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGAAGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTa                                                                                                                >  2:241449/1‑150 (MQ=255)
                                    aTGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCATTGTGGCGACGCAGGTAGTCa                                                                                                            <  2:1230712/150‑1 (MQ=255)
                                     tGTCGGCGGTGATAGTGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATgcgc                                                                                                                                                <  1:403331/113‑1 (MQ=255)
                                     tGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATgcgc                                                                                                                                                >  2:403331/1‑113 (MQ=255)
                                                aTACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGccc                                                                                                <  1:1036487/150‑1 (MQ=255)
                                                          gcACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGATGCTGGGGTTGGTGGCGTCGAGGGTGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCa                                                                                      >  2:627875/1‑150 (MQ=255)
                                                           cACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAg                                                                                     <  2:990256/150‑1 (MQ=255)
                                                              tGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGa                                                                                  >  1:423685/1‑150 (MQ=255)
                                                                       aTTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATTCGCTCGGTGTCGAGTTCCTtg                                                                                                                              <  1:397912/97‑1 (MQ=255)
                                                                       aTTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATTCGCTCGGTGTCGAGTTCCTtg                                                                                                                              <  1:397873/97‑1 (MQ=255)
                                                                       aTTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATTCGCTCGGTGTCGAGTTCCTtg                                                                                                                              >  2:397912/1‑97 (MQ=255)
                                                                       aTTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATTCGCTCGGTGTCGAGTTCCTtg                                                                                                                              >  2:397873/1‑97 (MQ=255)
                                                                           gcccccgcgccCTGCAGGCATCGTTAGGTGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGTAGACTGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCg                                                                     <  1:1043742/144‑1 (MQ=255)
                                                                            tCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCgg                                                                    >  2:1057109/1‑150 (MQ=255)
                                                                               tGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGt                                                                                                                 >  2:454171/1‑102 (MQ=255)
                                                                               tGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGt                                                                                                                 <  1:454171/102‑1 (MQ=255)
                                                                                        gcaggcagCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCGGTTCGGCATGGCCGTCGc                                                        >  2:673311/1‑150 (MQ=255)
                                                                                                     gATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTg                                           >  2:221143/1‑150 (MQ=255)
                                                                                                      aTGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTgg                                          >  2:829884/1‑150 (MQ=255)
                                                                                                          tgctgGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCa                                      >  1:1038527/1‑150 (MQ=255)
                                                                                                                  ttGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGa                                                                                  >  2:225771/1‑98 (MQ=255)
                                                                                                                  ttGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGa                                                                                  <  1:225771/98‑1 (MQ=255)
                                                                                                                     gtggCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGc                           <  1:88019/150‑1 (MQ=255)
                                                                                                                      tggCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTggg                                         <  2:188533/135‑1 (MQ=255)
                                                                                                                      tggCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTggg                                         >  1:188533/1‑135 (MQ=255)
                                                                                                                        gCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAg                        <  1:1057109/150‑1 (MQ=255)
                                                                                                                             gAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCt                   <  1:1170153/150‑1 (MQ=255)
                                                                                                                               gggCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGa                 >  2:995627/1‑150 (MQ=255)
                                                                                                                                ggCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGAc                >  2:1248432/1‑150 (MQ=255)
                                                                                                                                  cGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGACCt              >  1:1206260/1‑150 (MQ=255)
                                                                                                                                     agaagTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGACCTCGa           <  2:423685/150‑1 (MQ=255)
                                                                                                                                         gTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTNATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGACCTCGAGCAt       >  1:550793/1‑150 (MQ=255)
                                                                                                                                           gcagcaTGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGACCTCGAGCATGc     <  2:1038527/150‑1 (MQ=255)
                                                                                                                                             agcaTGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGACCTCGAGCATGCGt   <  2:1206260/150‑1 (MQ=255)
                                                                                                                                              gcaTGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGAGAGTAGCTGACCTCGAGCATGCGTg  >  1:266059/1‑150 (MQ=255)
                                                                                                                                                  gcgcTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGa                         >  2:672350/1‑123 (MQ=255)
                                                                                                                                                  gcgcTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGa                         <  1:672350/123‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGCATTTCGATCCAGGTGGCGTTGATGTTCTCCCACATGTCGGCGGTGATACGGCCGCGCACTGCATGGGCATTGGTCCTGGCGGCCTGCAGGCAGCAGTAGATGCTGCTGGGGTTGGTGGCGTCGAGGGCGAAGAAGTGCAGCATGCGCTCGGTGTCGAGTTCCTTGTGGCGACGCAGGTAGTCATCCAGCATGCCCGAGGCCATCAGCGACATTGCCAGTTCGGCATGGCCGTCGCTGCGCCCGGCCTGGGGCATCAGCGATAGCGAGTAGCTGACCTCGAGCATGCGTG  >  NC_002947/3077537‑3077828

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