Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,504 A→T 14.3% intergenic (‑54/+102) PP_1662 ← / ← PP_1663 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,5040AT14.3% 68.3 / 5.6 28intergenic (‑54/+102)PP_1662/PP_1663hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/8);  new base T (2/2);  total (18/10)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGT  >  NC_002947/1859356‑1859638
                                                                                                                                                    |                                                                                                                                      
ccTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACaa                                                                                                                                       <  2:581938/150‑1 (MQ=255)
       gggACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCAt                                                                                                                                >  1:273978/1‑150 (MQ=255)
       gggACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCAt                                                                                                                                >  1:97832/1‑150 (MQ=255)
                cagagtcgtcgCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGa                                                                                                                       >  1:21716/6‑150 (MQ=255)
                                      ttGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCt                                                                                                 >  1:151408/1‑150 (MQ=255)
                                                    gTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGTTTCAt                                                                                                                                <  2:123717/105‑1 (MQ=255)
                                                    gTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACTTGTTTCAt                                                                                                                                >  1:123717/1‑105 (MQ=255)
                                                     tAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggtgg                                                                                                    >  2:96023/1‑132 (MQ=255)
                                                     tAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCtgggtgg                                                                                                    <  1:96023/132‑1 (MQ=255)
                                                         cTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCt                                                                                                          <  1:327613/122‑1 (MQ=255)
                                                         cTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCt                                                                                                          >  2:327613/1‑122 (MQ=255)
                                                          tCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGc                                                                             >  1:179450/1‑150 (MQ=255)
                                                                     aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTggg                                                                  >  1:218118/1‑150 (MQ=255)
                                                                        aTGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                               <  2:273978/150‑1 (MQ=255)
                                                                               cGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCa                                                        <  2:97832/150‑1 (MQ=255)
                                                                                          cGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCTTAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcagg                                             <  2:21716/150‑1 (MQ=255)
                                                                                                tcttctTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgc                                       <  1:351926/150‑1 (MQ=255)
                                                                                                 cttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgcg                                      >  2:365644/1‑150 (MQ=255)
                                                                                                    cttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCGGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGa                                   <  2:151408/150‑1 (MQ=255)
                                                                                                            agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgggg                                                             <  1:309667/116‑1 (MQ=255)
                                                                                                            agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgggg                                                             >  2:309667/1‑116 (MQ=255)
                                                                                                            agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                           >  1:276320/1‑150 (MQ=255)
                                                                                                            agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTGTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCATGTATTGCAGGCAGTTCGCGATCAGTTga                           >  2:129976/1‑149 (MQ=255)
                                                                                                                  ttAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCa                     >  1:430438/1‑150 (MQ=255)
                                                                                                                     gCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCg                  >  2:371154/1‑150 (MQ=255)
                                                                                                                           taagtaagTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGc            >  2:602965/1‑150 (MQ=255)
                                                                                                                                     aaaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGt  >  2:302502/1‑150 (MQ=255)
                                                                                                                                         cttcATGAAACTTGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcag                                              >  2:344044/4‑102 (MQ=255)
                                                                                                                                         cttcATGAAACTTGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcag                                              <  1:344044/99‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
CCTTCGTGGGACATTCTCGATTCGTCGCGTTCGGTTACTTGGACAGTCATCGGTAAACTCCACTTCTTCAAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGT  >  NC_002947/1859356‑1859638

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