Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,166,687 T→A 15.6% intergenic (‑488/‑293) PP_5483 ← / → PP_1921 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,166,6870TA15.6% 77.4 / 6.4 32intergenic (‑488/‑293)PP_5483/PP_1921hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base A (3/2);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TTTTTTCCCCGGGTGAGCATTCAGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGA  >  NC_002947/2166542‑2166826
                                                                                                                                                 |                                                                                                                                             
ttttttCCCCGGGTGAGCATTCAGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAg                                                                                                                                           >  1:176277/1‑150 (MQ=255)
  ttttCCCCGGGTGAGCATTCAGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCg                                                                                                                                         <  2:486383/150‑1 (MQ=255)
        ccGGGTGAGCATTCAGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGgatga                                                                                                                                   >  2:198850/1‑150 (MQ=255)
                gCATTCAGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAggg                                                                                                                           <  1:558710/150‑1 (MQ=255)
                  aTTCAGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCt                                                                                                                         <  2:150820/150‑1 (MQ=255)
                      aGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCAGGAGCgggt                                                                                                                                      <  2:553908/133‑2 (MQ=255)
                      aGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCAGGAGCgggt                                                                                                                                      >  1:553908/1‑132 (MQ=255)
                                 ggCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTg                                                                                                          <  1:198850/150‑1 (MQ=255)
                                            gTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAg                                                                                               <  1:76755/150‑1 (MQ=255)
                                                     cGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCtcaatcaa                                                                                      >  1:17421/1‑150 (MQ=255)
                                                      gTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAAc                                                                                     >  1:271210/1‑150 (MQ=255)
                                                              ccTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTc                                                                           >  2:415243/1‑150 (MQ=255)
                                                                    tCAGAAAGTTACTCCCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGtt                                                                     <  1:221049/150‑1 (MQ=255)
                                                                         aaGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACTACCCGCTCCTGGAGCCGAATGATCCTAGCCCTGTTGGGGTTTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGCTTACAgc                                                                <  2:57675/150‑2 (MQ=255)
                                                                                tctcAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagg                                                         <  1:415243/150‑1 (MQ=255)
                                                                                           tATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtagggagg                                              >  1:384449/1‑149 (MQ=255)
                                                                                           tATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtag                                                   <  1:212910/145‑1 (MQ=255)
                                                                                           tATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtag                                                   <  1:212875/145‑1 (MQ=255)
                                                                                           tATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtag                                                   >  2:212910/1‑145 (MQ=255)
                                                                                           tATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtag                                                   >  2:212875/1‑145 (MQ=255)
                                                                                             ttGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtagggagtag                                            <  2:100690/150‑1 (MQ=255)
                                                                                                      gTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGACCTGCCGGTGCGGAGGGCTGTGGGGGTGTGCCTGGTTGAGTGAAGC‑‑‑ATATACCTGGCTGTTCCAGGTTGCAGGTCgtagggagtagggagtagggagtaggga                                  <  1:417784/150‑1 (MQ=255)
                                                                                                                 gCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTg                                                                                                          <  2:128939/70‑1 (MQ=255)
                                                                                                                 gCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTg                                                                                                          >  1:128939/1‑70 (MQ=255)
                                                                                                                         ccATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTGGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTGGGGAGTAGGGAGTAGGGGCAGGGGGAGTc                >  1:130097/1‑150 (MQ=255)
                                                                                                                                gtgAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCgg         <  2:17421/150‑1 (MQ=255)
                                                                                                                                gtgAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCgg         >  2:144746/1‑150 (MQ=255)
                                                                                                                                  gAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtag                                                   <  1:302982/106‑1 (MQ=255)
                                                                                                                                  gAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtagggagtag                                                   >  2:302982/1‑106 (MQ=255)
                                                                                                                                       atcccGCTCCAGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGa  >  2:258334/3‑150 (MQ=255)
                                                                                                                                        tcccGCTCCAGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtag                                                          >  2:165108/2‑93 (MQ=255)
                                                                                                                                        tcccGCTCCAGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCgtag                                                          <  1:165108/92‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TTTTTTCCCCGGGTGAGCATTCAGTGATTAGCTGGCAATGGCGGGTGAGTTGACGTCATCTCCCTCGATCAGAAAGTTACTCTCAGCAACCTATTGCAGAGGGTTTAAGGATCGCCGTTGGCCATTATGTGAAACGACCCGCTCCTGGAGCGGGATGATCCGAGGGCTGTTGGGGTGTCGCTGGTTGAGTGAAGCTCAATCAACCTGGC‑‑TTCCAGGTTGCAGGTCGTAGGGAGTAGGGAGTAGGGAGTAGGGAGTAGGGGCAGGGAGAGTCGGCGCGGGTGCAGA  >  NC_002947/2166542‑2166826

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