Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,076,095 C→A 9.6% intergenic (‑43/+31) PP_2685 ← / ← PP_2686 hypothetical protein/transglutaminase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,076,0950CA 9.6% 108.5 / 4.1 42intergenic (‑43/+31)PP_2685/PP_2686hypothetical protein/transglutaminase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (18/20);  new base A (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

GAACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGTGA  >  NC_002947/3075948‑3076240
                                                                                                                                                   |                                                                                                                                                 
gAACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTcg                                                                                                                                                  >  7:223965/1‑149 (MQ=255)
           tGAACGTGTAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGAc                                                                                                                                       <  8:223965/149‑1 (MQ=255)
            gAACGTGGAGATCTGGGCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCGCGAGTGCGGCAAGAcc                                                                                                                                      <  6:329447/149‑1 (MQ=255)
                    agaTCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGaa                                                                                                                              <  7:52926/149‑1 (MQ=255)
                       tCTGGTCGATGCCGTCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAg                                                                                                                             <  8:225174/147‑1 (MQ=255)
                                  tcGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTACAGGTGCATGGCGACACATTAAGTCATGCGGTTTCGGTTGTGTTGGGTGTGGAAGGGGCTTGCCGCACTCGCGCGTCAAGACCGGTCTGAAGCCAGTCTGTCAc                                                                                                                 <  3:44366/147‑1 (MQ=255)
                                   cGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCGCGAGTGCGGCAAGAcc                                                                                                                                      >  7:236084/1‑126 (MQ=255)
                                   cGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCGCGAGTGCGGCAAGAcc                                                                                                                                      <  8:236084/126‑1 (MQ=255)
                                    ggCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGt                                                                                                                        <  8:136615/139‑1 (MQ=255)
                                    ggCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGt                                                                                                                        >  7:136615/1‑139 (MQ=255)
                                                    tCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGa                                                                                                                                        <  6:60363/107‑1 (MQ=255)
                                                        tgatgaAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCa                                                                                                                                           >  3:133367/1‑100 (MQ=255)
                                                        tgatgaAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCa                                                                                                                                           <  4:133367/100‑1 (MQ=255)
                                                        tgatgaAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGc                                                                                          >  1:102218/1‑149 (MQ=255)
                                                           tgaAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCt                                                                                        <  2:102218/148‑1 (MQ=255)
                                                                  aaGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCCCGCGCGCGGCCAGCCCGGTCTGCAGCCCGTCTGTCACTGGCGATGCCCGGGAACGCGGGGCGGCagggct                                                                                >  6:33697/1‑144 (MQ=255)
                                                                    ggccGTCCTCCAGGTGCATGGCGACACAGTAAGTCTTGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTGTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTc                                                                              <  4:243782/147‑1 (MQ=255)
                                                                      ccGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGAt                                                                                                          >  7:269150/1‑119 (MQ=255)
                                                                      ccGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGAt                                                                                                          <  8:269150/119‑1 (MQ=255)
                                                                       cGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTg                                                                                                               >  4:123298/1‑113 (MQ=255)
                                                                        gTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGgcg                                                                          >  6:81873/1‑149 (MQ=255)
                                                                            tccAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAAc                                                                                                <  6:166696/123‑1 (MQ=255)
                                                                            tccAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAAc                                                                                                >  5:166696/1‑123 (MQ=255)
                                                                                     aTGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCGCTGGCGATGCACGTGAAAGCTGGCCTGGATCGATTCGGCGCCGCCGCCgcggc                                                             >  5:352128/1‑149 (MQ=255)
                                                                                       gtcGACACTGTAATTCAGGGGGTTTCCTTTGTGTTGTTTGTGGATGGGTCTTGCCGAACGCGGGCGTCATGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGt                                                           <  6:58931/147‑1 (MQ=255)
                                                                                             acaGTACGTCATGCGGTTGCCCTTGGGGTGCGTGTGCACGGGGCTTGCCGCCCGCCCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGccc                                                     <  5:33697/149‑1 (MQ=255)
                                                                                             acaGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTTCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCAc                                                                                                      >  4:186028/1‑100 (MQ=255)
                                                                                             acaGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCAc                                                                                                      <  3:186028/100‑1 (MQ=255)
                                                                                                gTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGACGCACTCGCGCGGCAAGACAGGTCTGAAGCCAGTCTGTCACTGGAGATGCACGTGAACGATGGCCTGCATAGATTCGGCGCCGCCGCCGCGGCGTACACCCCt                                                   >  1:221337/1‑148 (MQ=255)
                                                                                                       aTGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCggg                                            >  8:46794/1‑148 (MQ=255)
                                                                                                                 ccTTGTGTTGCGTGTGGAAGGGTGTTTGCCCCCGCGCGCCGCAAAACCCGTCAGAAGCCAATATGTGACTGGCGATGCCCCTTAAACCAGGGCGGCATAGCTTAGGCGCCCCCCCCGCGGCGTGACCCCCCGACCCGGCCGGcaacaat                                 >  6:38735/1‑148 (MQ=255)
                                                                                                                 ccTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTAc                                                         >  5:253834/1‑125 (MQ=255)
                                                                                                                 ccTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTAc                                                         <  6:253834/125‑1 (MQ=255)
                                                                                                                 ccTTGTGTTGCGTGTGGAAGGGTCATTCCGCCCTCGCGCGGCAAGACCGGTCTGAAGCCAGTAATTCACTGGCGAATGACCTGAAACCAGGCCTGCCACCATTCGGCGCCGCCCCCCCCGCCTACGCCCCTGGCCCggcaggcaaaa                                   >  2:141804/1‑144 (MQ=255)
                                                                                                                  cTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCCCGGCGTACGCCCCTGACCGGGCAGGCATCAAgg                                >  8:75540/1‑149 (MQ=255)
                                                                                                                              gtgGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGccc                    <  5:81873/149‑1 (MQ=255)
                                                                                                                                 gAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCCGGCAGGCATCAAGGTAGTCCAGGCCCAcc                 >  7:201486/1‑149 (MQ=255)
                                                                                                                                      gTCTTGCCGCGCGAGTGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAg            >  2:185707/1‑149 (MQ=255)
                                                                                                                                            ccGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAg      <  7:46794/149‑1 (MQ=255)
                                                                                                                                              gCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGt    <  1:141804/149‑1 (MQ=255)
                                                                                                                                                aCTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGTGa  <  5:38735/149‑1 (MQ=255)
                                                                                                                                                  tCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGt    <  6:148349/145‑1 (MQ=255)
                                                                                                                                                  tCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGt    >  5:148349/1‑145 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAACAACTTGGTGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGTGA  >  NC_002947/3075948‑3076240

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