Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,064,254 Δ1 bp coding (235/1218 nt) brpA ← LytR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,064,2530C.96.2% 116.7 / ‑1.7 26coding (236/1218 nt)brpALytR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (14/11);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCTTGAG  >  CP000730/1064108‑1064376
                                                                                                                                                 |                                                                                                                           
tgcatgGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCttt                                                                                                                          <  1:64263/149‑1 (MQ=255)
                      tAGCCAACTGTAGGAAGGTAGCTTAGAGTATCACGGGGTATACTAATCATTCTAATTTGATGTTTTAATTGATTCAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCTTCCTGTCTTTTTTTCTTCTTCCATCTTTatct                                                                                                   >  2:119581/1‑148 (MQ=255)
                       aGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTatcat                                                                                                  <  1:263162/149‑1 (MQ=255)
                       aGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTatcat                                                                                                  >  1:320413/1‑149 (MQ=255)
                          cAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGa                                                                                               >  1:292992/1‑149 (MQ=255)
                          cAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGa                                                                                               >  1:325939/1‑149 (MQ=255)
                             cTTTAGGACTGTAGCTGATAGTATCACGTGGAATACTCAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATTGCGTCTGATCTAGAGTGTTCCGCATCCTGTCTTATTTTCTTCTACCATCGTTATCATCGATAc                                                                                            <  2:220316/149‑1 (MQ=255)
                             cTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATAc                                                                                            <  1:307751/149‑1 (MQ=255)
                              tttAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATAcc                                                                                           <  1:247263/149‑1 (MQ=255)
                              tttAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATAcc                                                                                           <  2:85103/149‑1 (MQ=255)
                              tttAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATAcc                                                                                           <  1:58249/149‑1 (MQ=255)
                                   gAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTaaga                                                                                      >  2:208791/1‑149 (MQ=255)
                                    aaTGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTaagaa                                                                                     >  2:214431/1‑149 (MQ=255)
                                           cttaTAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATa                                                                              >  1:121684/4‑149 (MQ=255)
                                                                cTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCTATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGa                                                         >  2:127064/1‑149 (MQ=255)
                                                                cTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGa                                                         <  2:251940/149‑1 (MQ=255)
                                                                    gCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGa                                                     >  2:112343/1‑149 (MQ=255)
                                                                     cATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGAc                                                    >  1:123770/1‑149 (MQ=255)
                                                                          tAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTttat                                               <  1:107236/149‑1 (MQ=255)
                                                                                          ttgattAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGa                               >  1:163586/1‑149 (MQ=255)
                                                                                                aaaaGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACg                         <  2:309359/149‑1 (MQ=255)
                                                                                                             aTCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTTTTTTTTCATCGTTATCATCGACACATAAGAAAAGAATAGTGATAGGTTGTTCTTCGGGTTTTACTTTATTATCTCTTAATTTGGTTTGACGATTAGCATTTTTg            >  2:132201/1‑149 (MQ=255)
                                                                                                                   gCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGttt      >  1:12102/1‑147 (MQ=255)
                                                                                                                   gCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCt      <  1:294117/149‑1 (MQ=255)
                                                                                                                    cGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCtt     >  1:176406/1‑149 (MQ=255)
                                                                                                                       cTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCTTgag  <  1:214431/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                           
TGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCTTGAG  >  CP000730/1064108‑1064376

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