Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,893,643 A→T 15.4% intergenic (+118/‑90) pyk → / → PP_4302 pyruvate kinase/urea transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,893,6430AT15.4% 60.8 / 4.7 26intergenic (+118/‑90)pyk/PP_4302pyruvate kinase/urea transporter
Reads supporting (aligned to +/- strand):  ref base A (14/8);  new base T (2/2);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAAC  >  NC_002947/4893506‑4893789
                                                                                                                                         |                                                                                                                                                  
ggATCGAGACTTTGATCTGCTTTTGCCGGGGGCGCGTTGCGCCCCTTGCGCCACACACGGCCGCTCCCACCGGGCCACGCTGCCACAAGGCAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCTTTGccctccaatc                                                                                                                                         <  4:225436/149‑5 (MQ=255)
         cTTTGATCTTAGTTTGCGGGGGGGGGGTTGCGCCCCTTTCGCCACACATGCCCGCTCCCACCGGGCCACGCTGGCACATGGTATTGCAGCCCTTGGGGTTGGGGCTTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctc                                                                                                                                <  8:45789/149‑1 (MQ=255)
                                                      acaAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTccc                                                                                   >  6:329996/1‑149 (MQ=255)
                                                                tCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCt                                                                         >  4:111785/1‑149 (MQ=255)
                                                                  ccACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCTTTGccctccaat                                                                                                                                          >  3:270389/1‑79 (MQ=37)
                                                                  ccACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCTTTGccctccaat                                                                                                                                          <  4:270389/82‑4 (MQ=37)
                                                                  ccACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                                                           >  3:73161/1‑129 (MQ=255)
                                                                  ccACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                                                           <  4:73161/129‑1 (MQ=255)
                                                                                       aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTacag                                                   >  1:117670/1‑147 (MQ=255)
                                                                                       aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                      <  2:117670/145‑1 (MQ=255)
                                                                                       aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAcc                                                  <  5:329996/149‑1 (MQ=255)
                                                                                                 gTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCg                                        >  8:40148/1‑149 (MQ=255)
                                                                                                  tCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGt                                       >  6:206053/1‑149 (MQ=255)
                                                                                                            aGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                      <  5:203309/124‑1 (MQ=255)
                                                                                                            aGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                      >  6:240211/1‑124 (MQ=255)
                                                                                                            aGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                      >  6:203309/1‑124 (MQ=255)
                                                                                                            aGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                      <  5:240211/124‑1 (MQ=255)
                                                                                                                  cTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGAc                       >  5:79396/1‑149 (MQ=255)
                                                                                                                     gtgtCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                    >  5:46976/1‑149 (MQ=255)
                                                                                                                       gtCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGc                  >  4:256555/1‑149 (MQ=255)
                                                                                                                       gtCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGc                  >  2:90419/1‑149 (MQ=255)
                                                                                                                        tCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTgg                    <  3:188163/146‑1 (MQ=255)
                                                                                                                          gcgATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCAc               >  8:160778/1‑149 (MQ=255)
                                                                                                                               gggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGACGACCACCGATCCCCCCTCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt          >  5:228294/2‑149 (MQ=255)
                                                                                                                                 gAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTa        <  6:79396/149‑1 (MQ=255)
                                                                                                                                    ggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACt      >  1:78367/1‑148 (MQ=255)
                                                                                                                                       aaaGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAAc  <  7:40148/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAAC  >  NC_002947/4893506‑4893789

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