Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,322,225 G→A 9.0% intergenic (‑191/‑1232) PP_3792 ← / → PP_5631 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,322,2250GA 9.0% 116.3 / 3.5 46intergenic (‑191/‑1232)PP_3792/PP_5631hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/23);  new base A (2/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

GCTGCGCCGGGATGAGGTCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAGGGAATCTCCGAACAA  >  NC_002947/4322081‑4322361
                                                                                                                                                |                                                                                                                                        
gCTGCGCCGGGATGAGGTCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCag                                                                                                                                     >  2:17666/1‑150 (MQ=255)
gCTGCGCCGGGATGAGGTCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACAAAGCCGCGATAAGGCAAAGACGGAATTGAGAACAGCCGGTCTCTGGGCGCACTggcaagggcaag                                                                                                                                     >  2:738442/1‑150 (MQ=255)
 cTGCGCCGGGATGAGGTCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCaga                                                                                                                                    <  2:825261/150‑1 (MQ=255)
     gcCGGGATGAGGTCGACGTCGACAAAAACACGGCTCACATGTGCCGGAACGTGGTTGAAAAGCAGGCAGTCGCCCCGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAgagaga                                                                                                                                <  1:916177/150‑1 (MQ=255)
       cGGGATGAGGTCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATc                                                                                                                              <  1:663356/150‑1 (MQ=255)
                 tCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCagaga                                                                                                                                  <  2:333930/136‑1 (MQ=255)
                 tCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCagaga                                                                                                                                  >  1:333930/1‑136 (MQ=255)
                                   cGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGcgcc                                                                                                  <  1:1094389/150‑1 (MQ=255)
                                       tCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGt                                                                                              >  1:592284/1‑150 (MQ=255)
                                       tCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATATTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCTCCGCGt                                                                                              <  1:624067/150‑1 (MQ=255)
                                           aTCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCTCTCGAAGCAGAAGCGCCGCGTCTCg                                                                                          >  1:731001/1‑150 (MQ=255)
                                             cTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAAGGCCAGAGAGACCGCGCACTCCCGCTCGAAGCAGAAGCGCCGCGTCTCGAc                                                                                        <  1:304902/150‑1 (MQ=255)
                                              tGCCGGACCGCGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTAGAACCGCCGGGCCCTTGGCGCACCGGCAAGCGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCCCGACg                                                                                       >  2:324166/1‑150 (MQ=255)
                                                 cGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAg                                                                                    <  2:592284/150‑1 (MQ=255)
                                                   gATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACAGTATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTc                                                                                  <  1:738442/150‑1 (MQ=255)
                                                                caGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTc                                                                     >  1:777493/2‑150 (MQ=255)
                                                                 aGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGGCGCCATTCATCTTcc                                                                    >  2:890820/1‑150 (MQ=255)
                                                                      gcagTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCt                                                               >  1:439374/1‑150 (MQ=255)
                                                                         gTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGAc                                                                                        <  1:320522/122‑1 (MQ=255)
                                                                         gTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGAc                                                                                        >  2:320522/1‑122 (MQ=255)
                                                                            gAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGc                                                         >  1:798103/1‑150 (MQ=255)
                                                                            gAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGc                                                         >  1:1056173/1‑150 (MQ=255)
                                                                                   ccAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTcac                                                  <  1:17666/150‑1 (MQ=255)
                                                                                   ccAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTcac                                                  >  2:1016386/1‑150 (MQ=255)
                                                                                       gACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATc                                                                        >  1:444357/1‑124 (MQ=255)
                                                                                       gACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATc                                                                        <  2:444357/124‑1 (MQ=255)
                                                                                          gAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAAt                                           >  1:849433/1‑150 (MQ=255)
                                                                                            aaTACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGAc                                                                                        >  2:353652/1‑103 (MQ=255)
                                                                                            aaTACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGAc                                                                                        <  1:353652/103‑1 (MQ=255)
                                                                                                  cccGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCa                                        <  1:841198/145‑1 (MQ=255)
                                                                                                  cccGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCa                                        >  2:841198/1‑145 (MQ=255)
                                                                                                     gCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGAGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTc                                                       <  1:93079/127‑1 (MQ=255)
                                                                                                     gCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGAGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTc                                                       <  1:92766/127‑1 (MQ=255)
                                                                                                     gCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGAGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTc                                                       >  2:92766/1‑127 (MQ=255)
                                                                                                     gCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGAGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTc                                                       >  2:93079/1‑127 (MQ=255)
                                                                                                        ccGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCACAGGCCAGCGAGCTCGAGCAATACTGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGc                             <  2:969834/150‑1 (MQ=255)
                                                                                                                ttGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGc                     >  1:363207/1‑150 (MQ=255)
                                                                                                                  gAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCtt                   >  1:2873/1‑150 (MQ=255)
                                                                                                                   aaCAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTa                  <  2:39946/150‑1 (MQ=255)
                                                                                                                    aCAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAg                 <  2:1012372/150‑1 (MQ=255)
                                                                                                                     cAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGACCCCGCAATACCGCTCAACGCATAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAgg                <  1:1030851/150‑1 (MQ=255)
                                                                                                                        ccGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAGGGaa             <  2:849433/150‑1 (MQ=255)
                                                                                                                            gCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAGGGAAtctc         <  1:1016386/150‑1 (MQ=255)
                                                                                                                             ctctTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAGGGAATCTcc        <  1:285362/150‑1 (MQ=255)
                                                                                                                                   gcgcACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAGGGAATCTCCGAACaa  <  1:357896/150‑1 (MQ=255)
                                                                                                                                           gCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCt                            >  2:224391/1‑116 (MQ=255)
                                                                                                                                           gCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCt                            <  1:224391/116‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
GCTGCGCCGGGATGAGGTCGACGTCGAGAAAAAGACGGCTCACATCTGCCGGATCGTGGTTGAAAAGCAGGCAGTCGAGCGGACCAAGACGAAATACACCCGCACCGTAATGTTGAACAGCCGGGCTCTTGGCGCACTGGCAAGGGCCAGAGAGATCGCGCAATACCGCTCGAAGCAGAAGCGCCGCGTCTCGACGGAGTCGCCATTCATCTTCCAGCCTTCTGGCTCGTCACCACGAATCCAGCGTTCATTGCTGGGCAGCTTAGGGAATCTCCGAACAA  >  NC_002947/4322081‑4322361

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