Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,419,548:1 +C 100% intergenic (‑37/+99) 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.0% 76.3 / ‑3.1 25intergenic (‑37/+99)PP_1242/PP_1244hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (9/15);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGA  >  NC_002947/1419403‑1419635
                                                                                                                                                  |                                                                                        
ggCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggcct                                                                                       >  8:195731/1‑147 (MQ=255)
      ggATGGCGGGTTGTTTGCGGAAGAGCGGGTGAGGGTTGTTTTGGGCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGGTCTGTTGTGTTGTCGGTGGGAATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTcgc                                                                                 <  2:132377/149‑1 (MQ=255)
          gCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCggg                                                                              >  1:330151/1‑148 (MQ=255)
                ttCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGccc                                                                       <  5:311181/149‑1 (MQ=255)
                  cGAGGCGGAACAGCGTGGTAGTCTTGTTCGCGCCTTGGTAGAAGCAGTTCTTGCCGAATGTGGTGGTTTCTTTTGTGTTGGCTGTGGGCATTGTGTTGTTTCCCTTGATTCTGGGTTGCCTGGGCTGGCATCTTCGCTTGCACGccct                                                                      <  2:230292/148‑2 (MQ=255)
                          aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                             <  7:195731/149‑1 (MQ=255)
                                gggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtg                                                       >  8:333773/1‑149 (MQ=255)
                                    gAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTATTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                     >  6:19333/1‑147 (MQ=255)
                                      gggttgtgTTCGCCTTGGTAGAAGCAGGTCTGGCCGACTGTGGTGGGTGCTGTTGTGTTGTCGGGGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGt                                                 <  2:330151/147‑1 (MQ=38)
                                        gttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                     >  3:93063/1‑143 (MQ=255)
                                        gttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                     <  4:93063/143‑1 (MQ=255)
                                         ttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGc                                               <  4:274464/148‑1 (MQ=255)
                                            ttttCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGa                                           >  6:243575/1‑149 (MQ=255)
                                                cGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTATTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACgg                                        <  5:19333/148‑1 (MQ=255)
                                                 gCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCt                                      <  7:333773/149‑1 (MQ=255)
                                                  ccTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtg                                     >  2:301522/1‑149 (MQ=255)
                                                          gtaGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAc                              <  3:166197/146‑1 (MQ=255)
                                                          gAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAc                              <  1:301522/148‑1 (MQ=255)
                                                                  tCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGa                     <  5:243575/149‑1 (MQ=255)
                                                                    ttGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCggg                                                                              <  4:129859/90‑1 (MQ=255)
                                                                    ttGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCggg                                                                              >  3:129859/1‑90 (MQ=255)
                                                                              tggtggGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGccc                                                                       <  7:162290/87‑1 (MQ=255)
                                                                              tggtggGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGccc                                                                       >  8:162290/1‑87 (MQ=255)
                                                                                     ttCTGTTGGGTTGCCGGTGGGCATTGTTTTTTTTCCCTGGATTCTGGGTGGCCTGGGCTGGCCGCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTTACCCTGGGGTAGCGGGGGTGCCCGCGAaga  <  2:305739/149‑1 (MQ=37)
                                                                                                           aTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGc           <  4:327938/118‑1 (MQ=255)
                                                                                                           aTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGc           >  3:327938/1‑118 (MQ=255)
                                                                                                                                                  |                                                                                        
GGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGA  >  NC_002947/1419403‑1419635

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: