Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,419,549 +C intergenic (‑38/+98) PP_1242 ← / ← PP_1244 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,419,5481.C96.2% 82.0 / ‑2.9 26intergenic (‑37/+99)PP_1242/PP_1244hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (12/13);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAA  >  NC_002947/1419414‑1419676
                                                                                                                                       |                                                                                                                                 
ccGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGc                                                                                                                      <  3:171335/148‑1 (MQ=255)
 cGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGCCTGTGGTGGGTTCTGGTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCAc                                                                                                                    >  5:264451/1‑149 (MQ=255)
            cGGAACAGCGGGTGTGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTGGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATCCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTccc                                                                                                          <  1:74118/148‑1 (MQ=255)
               aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                      <  1:91960/149‑1 (MQ=255)
                            tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgt                                                                                               <  6:176879/143‑1 (MQ=255)
                            tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgt                                                                                               >  5:176879/1‑143 (MQ=255)
                              ttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGTTGGGCATTGTTTTGTTTCAATGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCg                                                                                       >  4:122753/1‑149 (MQ=39)
                                gttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGt                                                                                      <  2:193244/148‑1 (MQ=255)
                                     cGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                      <  7:178753/127‑1 (MQ=255)
                                     cGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                      >  8:178753/1‑127 (MQ=255)
                                           ggTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGTCtgtg                                                                            <  2:292128/147‑1 (MQ=255)
                                           ggTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtggt                                                                          >  8:302039/1‑149 (MQ=255)
                                              aGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGc                                                                                                              >  3:301335/1‑110 (MQ=255)
                                              aGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGc                                                                                                              <  4:301335/110‑1 (MQ=255)
                                               gAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAc                                                                       <  2:76576/148‑1 (MQ=255)
                                               gAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAc                                                                       >  1:76576/1‑148 (MQ=255)
                                                  gCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCAGCTCCCACAGGTACTGTTCAGTGCGTTGAGGCTGGGGTGCCCCtg                                                                   >  3:242570/1‑149 (MQ=38)
                                                     ggTGTTGCCGACTGTGGTGGGTTCTGTTGTTTTGTCGGTGGGCATGGTTTTGTTTCCCTGGATTCTGGGTTTCCTTGGCTGGCCTCTTCTCGGGCACGCCCTCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg                                                                <  8:215273/149‑1 (MQ=255)
                                                             cGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTAGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                      >  8:49280/1‑103 (MQ=255)
                                                             cGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTAGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                      <  7:49280/103‑1 (MQ=255)
                                                               aCTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGt                                                      <  7:302039/149‑1 (MQ=255)
                                                                                                 ttgttttgttTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGGAGTGCGAGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGGAAGgt                    >  3:113557/1‑149 (MQ=37)
                                                                                                         tttCCCTGTATTCGGGGTTGCGTGGGCTGGCCTGTGCGCGGGCACGCCGGCTCCCACCGTTAGTGTGCTGTGGGTGACGGCTGTGGTGACCCTGTGGGTGGGGGGGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAg               <  2:214924/146‑1 (MQ=255)
                                                                                                         tttCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAg               >  1:214924/1‑146 (MQ=255)
                                                                                                               tGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTgg      >  8:258029/1‑149 (MQ=255)
                                                                                                                    tcggggTGGCCTGGGCGGGCATCTTCTCGGGCACGCCCGCTCCCACCGGTCCTGTGCAGTGCGTGACGGCGGTGGTGACCCTGTGGGCGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTaa  <  4:242570/145‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
CCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAA  >  NC_002947/1419414‑1419676

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: