Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,960,142 T→A 29.2% intergenic (+50/‑78) bolA → / → PP_1758 DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,960,1420TA29.2% 34.0 / 15.4 24intergenic (+50/‑78)bolA/PP_1758DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base A (3/4);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTGCGCGA  >  NC_002947/1960006‑1960290
                                                                                                                                        |                                                                                                                                                    
cTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTAAGCGGcccccttt                                                                                                                                         <  1:307141/150‑5 (MQ=255)
 tGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGc                                                                                                                                                <  1:431159/142‑1 (MQ=255)
 tGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGc                                                                                                                                                >  2:431159/1‑142 (MQ=255)
            cccTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCaaa                                                                                                                             <  1:137786/149‑1 (MQ=255)
              cTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTAAGCGGccc                                                                                                                                              <  2:462641/131‑1 (MQ=255)
              cTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTAAGCGGccc                                                                                                                                              >  1:462641/1‑131 (MQ=255)
               tACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGt                                                                                                                          >  2:122006/1‑150 (MQ=255)
                         aaGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGt                                                                                                                <  2:226547/150‑1 (MQ=255)
                                      ggTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGa                                                                                                   <  1:122006/150‑1 (MQ=255)
                                        tCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAgg                                                                                                 <  1:193218/150‑1 (MQ=255)
                                                          cGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTc                                                                               <  1:366895/150‑1 (MQ=255)
                                                             ccGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAGCGATTGCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGt                                                                            >  1:6065/1‑150 (MQ=255)
                                                                    tgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGc                                                                                                                                                   <  2:449001/72‑1 (MQ=255)
                                                                    tgCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGc                                                                                                                                                   >  1:449001/1‑72 (MQ=255)
                                                                                          cGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGtt                                                                           <  2:559804/122‑1 (MQ=255)
                                                                                          cGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGtt                                                                           >  1:559804/1‑122 (MQ=255)
                                                                                               ttttGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGt                                          >  2:304288/1‑150 (MQ=255)
                                                                                                                  gCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCAtc                                                          <  1:487891/115‑1 (MQ=255)
                                                                                                                  gCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCAtc                                                          >  2:487891/1‑115 (MQ=255)
                                                                                                                   cATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                      <  1:304288/150‑1 (MQ=255)
                                                                                                                               ggggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                      <  2:335843/138‑1 (MQ=255)
                                                                                                                               ggggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                      <  2:335913/138‑1 (MQ=255)
                                                                                                                               ggggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                      >  1:335913/1‑138 (MQ=255)
                                                                                                                               ggggCCGCTAAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGaa                      >  1:335843/1‑138 (MQ=255)
                                                                                                                                       ttAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTGCGCGa  >  1:464490/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
CTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGGCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTGCGCGA  >  NC_002947/1960006‑1960290

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: