Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 27,584 A→G 19.3% intergenic (‑431/+15) PP_0022 ← / ← PP_0023 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294727,5840AG19.3% 45.7 / 6.5 26intergenic (‑431/+15)PP_0022/PP_0023hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base G (3/2);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCA  >  NC_002947/27450‑27730
                                                                                                                                      |                                                                                                                                                   
cACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCt                                                                                                                                      >  1:465305/1‑150 (MQ=255)
               cGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCGGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCtttttt                                                                                                                       >  2:327479/1‑149 (MQ=255)
                        ggTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTAACGACGCGGTACGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCGCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTgg                                                                                                              <  1:15574/150‑1 (MQ=255)
                         gTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTggg                                                                                                             >  1:142543/1‑150 (MQ=255)
                              ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTc                                                                                                        >  2:44698/1‑150 (MQ=255)
                              ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTc                                                                                                        <  2:142543/150‑1 (MQ=255)
                                   aaTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTc                                                                                                   <  2:317725/150‑1 (MQ=255)
                                                aagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTc                                                                                     <  1:44698/150‑1 (MQ=255)
                                                           tCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGCCCTTCTGCTTTTCAgg                                                                          <  1:523258/150‑1 (MQ=255)
                                                                   cgGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTc                                                                  <  2:429585/150‑1 (MQ=255)
                                                                               tCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGt                                                      >  2:609159/1‑150 (MQ=255)
                                                                                  gCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCg                                                   >  1:290134/1‑150 (MQ=255)
                                                                                    gagaCGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGa                                                 >  1:376746/1‑150 (MQ=255)
                                                                                      gaCGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGt                                               >  2:318868/1‑150 (MQ=255)
                                                                                                              gCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg                       <  2:465305/150‑1 (MQ=255)
                                                                                                                       aaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCg                                                   >  1:511090/1‑113 (MQ=255)
                                                                                                                       aaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCg                                                   <  2:511090/113‑1 (MQ=255)
                                                                                                                       aaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCtt              <  1:425770/150‑1 (MQ=255)
                                                                                                                         aaaGCAGCTCTCGGGAGCTGCTTTTTTTTTACTCCCCCTTTTTGTGCTTGTGAGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGGGACTCGCTGCCTTTTGTCCGGAGTGCTATCTTGATTGGTGCGCTTTGGGCCTTGCTTGa            >  2:276201/1‑150 (MQ=255)
                                                                                                                          aaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTTTTTTTTTGTGTTTTCAGGGCCGGATGATCCTTTTGCTTTTCAGGGGGATCTCGCTGCCTTTTGTCCGCACTGCCATCTTCATTAGTGCGCTTTGGCCCTGTCTTTAt           >  2:186901/1‑150 (MQ=255)
                                                                                                                          aaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGAt           <  1:403944/150‑1 (MQ=255)
                                                                                                                            gcagcTCCCTAGACCTGCTTTTTTCTTACTCCCCCTTTCTGTTCTTGTGGGGTTTCATGTTCGGATGATCCTTCTGCGTTTCAGTGGGGAGTCGCTGGCTTCTGTCCGGAGTGACAGCTTCTTTGGTGAGCTTTTGGCCTGGCTCGATGc         >  2:411694/1‑150 (MQ=255)
                                                                                                                                 tCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGa            >  1:503951/1‑142 (MQ=255)
                                                                                                                                 tCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGa            <  2:503951/142‑1 (MQ=255)
                                                                                                                                 tCCCGAGAGCTGCTTATTTCATACTTCTGCTTTTTTGTCTTGGGGGGTTTCAGGTCCGGATGGTCCTTTATGCTGTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATg    <  1:318868/150‑1 (MQ=255)
                                                                                                                                   ccGAGAGCTGCTTTTTTCTTACTTCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCa  >  2:590760/1‑150 (MQ=255)
                                                                                                                                      |                                                                                                                                                   
CACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCA  >  NC_002947/27450‑27730

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