Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,081,045 T→A 16.6% intergenic (+36/+26) PP_5332 → / ← PP_5333 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,081,0450TA16.6% 74.4 / 9.8 36intergenic (+36/+26)PP_5332/PP_5333hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base A (4/2);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

GAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTG  >  NC_002947/6080906‑6081188
                                                                                                                                           |                                                                                                                                               
gAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGAag                                                                                                                                       >  1:462557/1‑150 (MQ=255)
 aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                      >  2:790842/1‑150 (MQ=255)
 aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                      <  2:462557/150‑1 (MQ=255)
 aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                      <  1:790842/150‑1 (MQ=255)
       gtgCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGTGGGCTTGCCCGCGAAGAGGCCgg                                                                                                                                <  2:1046661/150‑1 (MQ=255)
                gaCAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCCTCCAGAACGGTAGTTAAAATATCTGCGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTAt                                                                                                                       <  1:134665/150‑1 (MQ=255)
                   aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTTTTTTCTTCTTCGCGGGCTTTCCCGCTAAGAGGCCGGGGCAATTCTCaa                                                                                                                    >  2:109658/1‑150 (MQ=255)
                   aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCaa                                                                                                                    <  1:129015/150‑1 (MQ=255)
                    gCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCATATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTACCCGCGAAGAGGCCGGGGCAATTATCgaa                                                                                                                   >  1:433431/1‑150 (MQ=255)
                      gAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGGCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTCCCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCACTTATCAAAgg                                                                                                                 >  2:993239/1‑150 (MQ=255)
                      gAATGTGTGGGATGACCTGGTGGAGTTGAAGCATGGGGAGAGTGGCTTGCACCACGTCCAGATCAGTTTTGAAAAAATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAgg                                                                                                                 >  1:21674/1‑150 (MQ=255)
                                            gAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCa                                                                                           >  2:220790/1‑150 (MQ=255)
                                                 gAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTg                                                                                      >  1:735730/1‑150 (MQ=255)
                                                       tGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCa                                                                                >  1:704133/1‑150 (MQ=255)
                                                                   ttCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAg                                                                    >  2:219233/1‑150 (MQ=255)
                                                                   ttCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAg                                                                    <  2:21674/150‑1 (MQ=255)
                                                                                cAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGCAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCa                                                       <  2:605230/150‑1 (MQ=255)
                                                                                 aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAg                                                      <  2:730537/150‑1 (MQ=255)
                                                                                       gTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGTCCTCTTCGCGGGCTTGCCCGCGCAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCAtt                                                >  2:709740/1‑150 (MQ=255)
                                                                                              aaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAAGCCGGGGCAATTATCAAAGGTCTGCGTCTTACACCACCCTCCCCTTGTCTGCATCCAGCGCATAGTCCCCATCGGCCAAGTAATTGCtctcc                                         >  2:549651/1‑146 (MQ=255)
                                                                                               aaTAGCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTc                                                                                                      >  2:885325/1‑88 (MQ=255)
                                                                                               aaTACCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTc                                                                                                      <  1:885325/88‑1 (MQ=255)
                                                                                                   tCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTnctc                                    <  2:704133/150‑1 (MQ=255)
                                                                                                       cgCCCTGCTCTTTTTTTCCCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCCCCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCAcc                                <  2:690416/149‑1 (MQ=255)
                                                                                                           cTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                            >  1:625698/1‑70 (MQ=255)
                                                                                                           cTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                            <  2:625698/70‑1 (MQ=255)
                                                                                                                       tACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGa                <  1:219233/150‑1 (MQ=255)
                                                                                                                         ccGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCaccacc                                                                                               >  1:248954/1‑69 (MQ=255)
                                                                                                                         ccGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCaccacc                                                                                               <  2:248954/69‑1 (MQ=255)
                                                                                                                               tctTCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc        >  2:1048892/1‑150 (MQ=255)
                                                                                                                                 ttCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgg      >  2:393398/1‑150 (MQ=255)
                                                                                                                                 ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgg      >  1:101244/1‑150 (MQ=255)
                                                                                                                                  tCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg     <  1:393398/150‑1 (MQ=255)
                                                                                                                                   cgcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATTGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgggg    <  2:735730/150‑1 (MQ=255)
                                                                                                                                    gcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCACAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCAGCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTCGCCGCTGACCCACAGCGGGgt   <  2:280520/150‑1 (MQ=255)
                                                                                                                                     cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGgtg  >  1:791558/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                               
GAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTG  >  NC_002947/6080906‑6081188

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: