Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,288,410 T→A 11.1% intergenic (‑91/+93) PP_4660 ← / ← mdaB LysR family transcriptional regulator/drug activity modulator B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,288,4100TA11.1% 91.4 / 3.7 36intergenic (‑91/+93)PP_4660/mdaBLysR family transcriptional regulator/drug activity modulator B
Reads supporting (aligned to +/- strand):  ref base T (14/18);  new base A (2/2);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

TCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCC  >  NC_002947/5288262‑5288557
                                                                                                                                                    |                                                                                                                                                   
tCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTT                                                                                                                                                     <  1:96135‑M1/149‑1 (MQ=255)
  gAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGC                                                                                                                                                   >  8:142521‑M1/1‑149 (MQ=255)
   aGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCC                                                                                                                                                  <  4:86043‑M1/149‑1 (MQ=255)
   aGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGGCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCC                                                                                                                                                  <  6:290209‑M1/149‑1 (MQ=255)
        cAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                                                                             <  8:236941‑M1/149‑1 (MQ=255)
         aGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCC                                                                                                                                            <  6:219751‑M1/149‑1 (MQ=255)
           tCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCA                                                                                                                                          <  2:349169‑M1/149‑1 (MQ=255)
              aTGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACaa                                                                                                                                       <  5:72167/149‑1 (MQ=255)
               tGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCT                                                                                                                                              >  8:96781‑M1/1‑141 (MQ=255)
               tGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCT                                                                                                                                              <  7:96781‑M1/141‑1 (MQ=255)
                 aCGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTa                                                                                                                                    <  4:239074/149‑1 (MQ=255)
                    gCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtatta                                                                                                                                 >  1:306034‑M1/1‑142 (MQ=255)
                     cgacgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTg                                                                                                                                 >  1:171396/1‑148 (MQ=255)
                     cgacgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTg                                                                                                                                 >  1:165804/1‑148 (MQ=255)
                            tACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTcc                                                                                                                         <  2:160429/149‑1 (MQ=255)
                              ccTGGAGTTCTTCGGATCGGGTTTTCGCATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                       <  2:153297/149‑1 (MQ=255)
                              ccTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCC                                                                                                                                            >  7:151048‑M1/1‑128 (MQ=255)
                              ccTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCC                                                                                                                                            <  8:151048‑M1/128‑1 (MQ=255)
                               cTGGAGTTGTTCGGTTGGGGTTTTGTCATCAGTGCCTGTTGTTGGTCTGTGTAAGGATTAGATATCCGGTTGTCGAGTCCAAGTCCAGCGGAGCCTGTATGTGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                       <  1:14402/148‑1 (MQ=255)
                                                  gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                   >  2:165804‑M2/1‑146 (MQ=255)
                                                  gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                   >  2:171396‑M2/1‑146 (MQ=255)
                                                  gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                   >  7:212589‑M2/1‑146 (MQ=255)
                                                    tttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaaga                                                                                                                                          <  1:139235‑M1/108‑6 (MQ=255)
                                                    tttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaaga                                                                                                                                          >  2:139235‑M1/1‑103 (MQ=255)
                                                     ttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGTCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgg                                                                                                 >  4:111362‑M2/1‑143 (MQ=255)
                                                       cacaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtggga                                                                                               <  4:342810‑M1/148‑42 (MQ=255)
                                                        acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCGGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagc                                                                                             <  3:309816‑M1/149‑44 (MQ=255)
                                                           tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcggg                                                                                          <  2:306034‑M1/149‑47 (MQ=255)
                                                                                      gATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCCCAACCACAACTAAATCAAAGTACATGAAAACAGCGTAAAAACAGTTGGAGCGGGCCATCCCGCCAAGAAGCCCGTAAAGa                                                                >  2:277791/1‑148 (MQ=255)
                                                                                                       tCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCACGCCCGCGAAGACGCCGGTAAAGACAGAAAATTACTGTTCCt                                              >  7:183298/1‑149 (MQ=38)
                                                                                                                                       gtgggagcgGGCAAGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGc                     <  6:30488/136‑1 (MQ=255)
                                                                                                                                       gtgggagcgGGCAAGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGc                     >  5:30488/7‑142 (MQ=255)
                                                                                                                                       ctgttCGGGGGCTTGCCCGCTCCCACAATTACTGAAAAGTCATTTAAACCTGCGTAATACCTGTGGTAGCGGGAAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATaa              <  1:277791/146‑1 (MQ=38)
                                                                                                                                              cgGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTAcc        >  7:116338/1‑148 (MQ=255)
                                                                                                                                                gtcTTGCCCGGTCCTTCACCTTCTGCAAGTTCGTTGAAACCTGCGTAATACCTTTGGTAGGGGGAAAGCGCTCTAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTTCCAAc     <  8:183298/147‑1 (MQ=255)
                                                                                                                                                 gCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAcg    >  2:311109/1‑149 (MQ=255)
                                                                                                                                                   ttGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAcgcc  >  4:188795/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCC  >  NC_002947/5288262‑5288557

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