Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,076,097 C→T 9.7% intergenic (‑45/+29) 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,0970CT 9.7% 105.3 / 3.2 42intergenic (‑45/+29)PP_2685/PP_2686hypothetical protein/transglutaminase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (18/19);  new base T (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 = 6.06e-01

TGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGTGA  >  NC_002947/3075959‑3076240
                                                                                                                                          |                                                                                                                                               
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                                                                                                >  5:166696/1‑123 (MQ=255)
                                                                 tccAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAAc                                                                                                <  6:166696/123‑1 (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                                                         <  6:253834/125‑1 (MQ=255)
                                                                                                      ccTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTAc                                                         >  5:253834/1‑125 (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)
                                                                                                                                          |                                                                                                                                               
TGAACGTGGAGATCTGGTCGATGCCGGCATTAGTGCGTGAGTCGCTGATGAAGACAAGCCCGTCCTCCAGGTGCATGGCGACACAGTAAGTCATGCGGTTTCCCTTGTGTTGCGTGTGGAAGGGTCTTGCCGCACTCGCGCGGCAAGACCGGTCTGAAGCCAGTCTGTCACTGGCGATGCACGTGAACGCTGGCCTGCATCGATTCGGCGCCGCCGCCGCGGCGTACGCCCCTGACCGGGCAGGCATCAAGGTAGTCCAGGCCCACCGCCAGCTTCAGGTGA  >  NC_002947/3075959‑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: ‑