Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7350AT32.3% 57.9 / 13.3 31Y2* (TAT→TAA) 
*65K (TAA→AAA) 		rcbA
ralR		double‑strand break reduction protein
endodeoxyribonuclease toxin RalR
Reads supporting (aligned to +/- strand):  ref base A (0/21);  new base T (10/0);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 2.25e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGT  >  NC_000913/1413499‑1413962
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
gatgatgatTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTATTTTTTTATATTAACACCTc                                                                                                                                                                                                                        >  2:128730/1‑251 (MQ=255)
    ttGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACa                                                                                                                                                                                                                                  <  2:342510/237‑1 (MQ=255)
     tGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTCATTTTTTACATTAacaaccccctc                                                                                                                                                                                                                    >  2:549307/1‑242 (MQ=255)
              ataatCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCa                                                                                                                                                                                                          <  1:150465/251‑1 (MQ=255)
                    aGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTTGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑tt                                                                                                                                                                                                   <  2:336612/251‑1 (MQ=255)
                       tATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAATAACACCCCCATTATTTTTTTCCGTGG‑tttgt                                                                                                                                                                                                >  1:552772/1‑248 (MQ=255)
                                      tGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACACCACCCTTATTATTTTCCGTGG‑TTCGtttttttcttcttttt                                                                                                                                                                                 >  2:10904/1‑238 (MQ=255)
                                             aaTGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAAGTACACCCCCCTTATTATTTTCAGTAG‑tttttttttttctcttttttaggtgtt                                                                                                                                                                          >  1:410077/1‑226 (MQ=255)
                                                                         cGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTc                                                                                                                                               <  1:10904/250‑1 (MQ=255)
                                                                            aaTTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGa                                                                                                                                           <  1:197217/251‑1 (MQ=255)
                                                                                   aTTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAACACCCCCCTTTTTATTTTCTATGG‑TTTGGTTTTTCCATCTTTTGAGTGCTTTTTTTTCTCTTCCCCCCTAAACGGTTCggggttttcat                                                                                                                                    >  1:430275/1‑243 (MQ=255)
                                                                                                aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCAGTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTAGTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCGCTTCCTCCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                                                                                                        <  2:512915/250‑1 (MQ=255)
                                                                                                    tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTc                                                                                                                    <  1:128730/250‑1 (MQ=255)
                                                                                                    tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCAATAACACCTCCCTTTTTATTTTCAGTGG‑GTTGTTTTTTTCCTTTTTTTAGTGCTTTTTTTTTATTTTCACCCTAACCCGTTCggggtttgcctccctttcaacttcctcc                                                                                                                   >  2:189048/1‑226 (MQ=255)
                                                                                                       gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCgg                                                                                                                  <  1:90119/249‑1 (MQ=255)
                                                                                                                 aaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTTGCTCGGTAATATCCTGaa                                                                                                      <  2:552772/251‑1 (MQ=255)
                                                                                                                             ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATTAACAACCCCCTTTTTTTTTTTATTGGGTTGTTTTTTTCCTTCTTT‑TAGTGCCTCTTTTTTTCTTCCACCCCTACCGGTTTGGGGTTCGCATCCGTTACACCTCTCTCgcgaaaaacacgaaaatgcgttctcc                                                                                          >  2:552116/1‑226 (MQ=255)
                                                                                                                              tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAc                                                                                           <  2:397879/249‑1 (MQ=255)
                                                                                                                                        agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAACACCCCCCTTTTTATTTTCACTGG‑GTGTTTTTTTTCATCTTTTTAGTGCTTCTTTTTTATTTCTCCCCTAAACGGGTCGGGGTTTGCCTCCGTTATACTTCTCCCCGGAATATCCTGAAATggggtttaccgctcttgcttg                                                                               >  2:551619/1‑232 (MQ=255)
                                                                                                                                             aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAAc                                                                          <  1:390565/251‑1 (MQ=255)
                                                                                                                                             aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAAc                                                                          <  2:262715/251‑1 (MQ=255)
                                                                                                                                                   tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACaa                                                                     <  1:405955/250‑1 (MQ=255)
                                                                                                                                                     acaatCTGAATTCGTTAGTTTAGATTTTGAATATCTTGTCGAGTTAGTAAGAGTTCCACCTTCCTTTACAAGAGTCGCGGTAATTTTATACATGAACACCTCCATTCTTATTTCCAGTGA‑TTCGTTTATTTCATCTTTCGAGTGCTTCTTTTTCACTTCCACCGTAACCGGTTCGCGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAAGGCTTACACCGTTACTGATGGACAACCACAAAAt                                                                  <  2:289229/247‑1 (MQ=255)
                                                                                                                                                            gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACAt                                                            <  2:17803/250‑1 (MQ=255)
                                                                                                                                                                      gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACAtt                                                 <  2:538620/251‑1 (MQ=255)
                                                                                                                                                                         tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGt                                               <  1:229095/250‑1 (MQ=255)
                                                                                                                                                                           gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATAAAAACCACCCTTATTATTTTCAGCGG‑GTCGTTTATTCCCTTTTTTTACTAATTTTTTTTTTCTTCCACCCCAAACGGTTCGGGATTTGCATCCCATACACTTCGCTCCGGAATATCCTTAAAAGGCTTTTACCGCTTCTCATGGGAAACCACaaaaaagaaaatggttaaaattgttat                                            >  1:252183/1‑228 (MQ=255)
                                                                                                                                                                                  aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCAt                                     <  1:549307/251‑1 (MQ=255)
                                                                                                                                                                                            ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACTCCTCCATTATTATTTCCAGTTG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTtctc                           <  1:52473/251‑1 (MQ=255)
                                                                                                                                                                                               gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATa                        <  1:345774/251‑1 (MQ=255)
                                                                                                                                                                                                                      tttCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGt  <  1:270107/250‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
GATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGG‑TTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGT  >  NC_000913/1413499‑1413962

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