Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,015,961 A→T 11.1% intergenic (‑176/+58) PP_5679 ← / ← PP_4421 hypothetical protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,015,9610AT11.1% 100.5 / 5.2 36intergenic (‑176/+58)PP_5679/PP_4421hypothetical protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (15/17);  new base T (2/2);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGA  >  NC_002947/5015815‑5016106
                                                                                                                                                   |                                                                                                                                                 
gatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGTTAGCACTACAAAAAAAAAAAGCCCGGCTAAAgg                                                                                                                                                >  1:394529/1‑150 (MQ=255)
 atgaATAAGCCCGCACCGGTCTCCGGCGGGGTCTTTTACGTTCGGTGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGc                                                                                                                                               >  1:19264/1‑150 (MQ=255)
           ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa                                                                                                                                     >  1:156232/1‑150 (MQ=255)
           ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa                                                                                                                                     >  1:160349/1‑150 (MQ=255)
               aCCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa                                                                                                                                 <  1:565746/150‑1 (MQ=255)
                  ggtacccGGCGCGGGCTTTAACTTTCCGGGCCTTGGACAGAATCTGACCCGAAAGTACGTCCTTGTGAAATGATGACAGCGTCTTGGCTTTGATTGTGATAGAAATACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAagag                                                                                                                              <  2:393293/145‑1 (MQ=255)
                             cgGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGc                                                                                                                   <  1:605154/150‑1 (MQ=255)
                                        tttCGGGGCCTTGGACAGTATCTGAACGGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTAATAGCAATACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTg                                                                                                        <  2:730349/150‑1 (MQ=255)
                                               gCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCCTAAAACGAAGAGGGAAATTCAGCTGGCATCATTGCCaaaga                                                                                                 >  1:819726/1‑147 (MQ=255)
                                                    ggACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTTTTGGGTTTCATCGGAAGAGGACGACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATg                                                                                            <  1:101520/150‑1 (MQ=255)
                                                     gACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                             >  2:237524/1‑116 (MQ=255)
                                                     gACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                             <  1:237524/116‑1 (MQ=255)
                                                     gACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATgg                                                                                           <  1:136595/150‑1 (MQ=255)
                                                      aCAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGc                                                                                          <  2:19264/150‑1 (MQ=255)
                                                                  ccGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAAAGGGTAATTCAGCTGGTATCATTGCCAACGTTTATGGTTCAGACCTTCaa                                                                              >  1:309585/1‑149 (MQ=255)
                                                                         aaCGTCCTTGTGAACTAAAGAAATCGTCTTGGCTTTGATTGTAATAAAAAAAAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGATGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtcatca                                                                       <  2:160349/149‑1 (MQ=255)
                                                                                ttGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGc                                                                >  1:767278/1‑150 (MQ=255)
                                                                                            gaCAGCGTCTTGGCTTTGATTGTGCTAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCtt                                                    >  1:240733/1‑150 (MQ=255)
                                                                                               aGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtt                                                                                                                                         <  1:439267/63‑1 (MQ=25)
                                                                                               aGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtt                                                                                                                                         >  2:439267/1‑63 (MQ=25)
                                                                                                 cGTCTTGGCTTAGAAAGTGCTCGACATAGAACACAACAAAGCCCGCCTAAATGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgtgt                                               <  2:309585/150‑1 (MQ=255)
                                                                                                    cTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGc                                             <  1:3062/150‑1 (MQ=255)
                                                                                                    cTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGc                                             >  2:3062/1‑150 (MQ=255)
                                                                                                          tttGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                             >  1:182199/1‑63 (MQ=255)
                                                                                                          tttCATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                             <  2:182199/63‑1 (MQ=255)
                                                                                                            tGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtgat                                    <  2:240733/150‑1 (MQ=255)
                                                                                                                    aTAGCACTACAAAAAAAAAAAAGCCCCCCTAAAGGCGGGCTTAAAACAAAAAGGGTAATTCACCTGGCACCATTGCCAACGTTGATGGCTCAGACCTTCACCCCATCAAGAACCCCCACCCAACCTTCGTGTCGGTTGTTGATGAGGCGa                             >  2:680696/1‑150 (MQ=255)
                                                                                                                        cACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCaa                                                                      >  1:633663/1‑104 (MQ=255)
                                                                                                                        cACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCaa                                                                      <  2:633663/104‑1 (MQ=255)
                                                                                                                                 aaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGtt                <  2:767278/150‑1 (MQ=255)
                                                                                                                                   aaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCg              <  1:680696/150‑1 (MQ=255)
                                                                                                                                          gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTc           <  2:403999/146‑1 (MQ=255)
                                                                                                                                          gCCCGCCTTTAGGCGGGCTTAAAACGAAGAGGGCAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTc           >  1:403999/1‑146 (MQ=255)
                                                                                                                                              gCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTc                     <  2:582221/132‑1 (MQ=255)
                                                                                                                                              gCCTAAAGGCGGGCTTAAAACGAAGAGGCTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTc                     >  1:582221/1‑132 (MQ=255)
                                                                                                                                               ccTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAaga  >  2:369810/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGA  >  NC_002947/5015815‑5016106

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