Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,616,570 C→A 15.5% Q102K (CAA→AAA)  ylyA → sporulation‑related RNA polymerase‑binding factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,616,5700CA15.5% 75.2 / 2.3 35Q102K (CAA→AAA) ylyAsporulation‑related RNA polymerase‑binding factor
Reads supporting (aligned to +/- strand):  ref base C (26/3);  new base A (2/4);  total (28/7)
Fisher's exact test for biased strand distribution p-value = 8.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCCTTGCATTGTCCAAAATCGAAAACGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCATA  >  BS168_reference_genome/1616440‑1616708
                                                                                                                                   |                                                                                                                                          
tctcCTTGCATTGTCCAAAATCGAAAACGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAAttt                                                                                                                                       >  5:275962/1‑136 (MQ=255)
tctcCTTGCATTGTCCAAAATCGAAAACGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCGCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATATTCTTTATTCTTTCCAAttt                                                                                                                                       >  5:23825/1‑136 (MQ=255)
tctcCTTGCATTGTCCAAAATCGAAAACGGCACATTCGGATACTGTGAACAAACTGGGGCCCCCAT‑TCCTCTAGCAAACCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATATTCTTTATTCTGTCAAAttt                                                                                                                                       >  3:67035/1‑136 (MQ=255)
      tGCATTGTCCAAAATCGAAAACGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGaaaaa                                                                                                                               >  6:110387/1‑138 (MQ=255)
            gTCCAAAATCGAAAACGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGCCGGCCCGAACAGCAAACGATTTTAGTTATTCTGTCAAATTTGAGAAAAAAACGCtt                                                                                                                        >  2:19357/1‑139 (MQ=255)
            gTCCAAAATCGAAAACGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCtt                                                                                                                        >  5:52364/1‑139 (MQ=255)
             tCCAAAATCGAAAACGTCACATACGGATACTGAGAAGAAACTGGGGCCCCCAT‑TACACTTTCAAAGCTCGCAGTGCTGCCGACGGCCCGAACAGAAAACGAATTTATATATTCTGGCCAAATTGAGAAAACCACGCTTc                                                                                                                       >  6:116412/1‑139 (MQ=255)
                       aaaaCGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGCCGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGa                                                                                                              >  6:246559/1‑138 (MQ=255)
                       aaaaCGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGa                                                                                                              >  1:59496/1‑138 (MQ=255)
                                         attgttAATAAAGTGGCCCCCCAAT‑TCCTCTTCCAAAGGGGGGTGGCGTGCCGAGGGCCCTAACCCCATAGGTTTTTATTTTTTCTGAAAATTTTTAGAAAATAACGCTTCCGATATGGAAATCAACGGATATCGAata                                                                                           <  2:59496/137‑1 (MQ=255)
                                             tgAAGAAACTGGGGCCCCCCT‑TGCGCTTGCAAAGCTGGCAGTGGTGGCGCCGGCCCCAGCAGCACACGCTTTTATTTATTCTGTCCCGTTTGAGGAAAAACAGCGTCCGTTATGGCAATCAACGGATATGGAGTATGGt                                                                                       >  2:142692/1‑139 (MQ=255)
                                             tgAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGt                                                                                       >  7:86441/1‑139 (MQ=255)
                                             tgAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGt                                                                                       >  3:106142/1‑139 (MQ=255)
                                             tgAAGAAACTGGGGCCCCCAT‑TCCGCTTGCAAAGCTGGCAGTGCTGCCGCCGGCCCGAACAGCCAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATTTTGCAATCAACGGATATCGAATATGGt                                                                                       >  6:177245/1‑139 (MQ=255)
                                                  caccgggggCCCCCATGTCCTCTTCCAACGCGGGCAGGGCTGCCGACGCCCCTAACACAAAAGGATTTTATTTTTTATGACAATTTTTAAAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAgg                                                                                   <  1:117446/134‑1 (MQ=255)
                                                     cTGGGGCCCCCAT‑TCCTGTTGCAACGCTGGAAGGGCGGCCTACGCCCCGAACACCAAACGATTTTCTTTATTCTGTCAATTTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACt                                                                               <  5:110387/139‑1 (MQ=255)
                                                                        ttGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAgg                                                                                   >  8:18876/1‑117 (MQ=255)
                                                                        ttGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAgg                                                                                   <  7:18876/117‑1 (MQ=255)
                                                                          gCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGtgtctgt                                                           >  3:300532/1‑139 (MQ=255)
                                                                          gCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGtgtctgt                                                           >  1:115899/1‑139 (MQ=255)
                                                                          gCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGtgtctgt                                                           >  6:238407/1‑139 (MQ=255)
                                                                                cTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTAATCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAAACAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGt                                                      >  6:277044/1‑138 (MQ=255)
                                                                                   gCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGATAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAat                                                  <  6:275962/139‑1 (MQ=255)
                                                                                          tgcCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATAAATAACAGTGTCTGTCCGGTTAATATg                                               >  5:301305/1‑135 (MQ=255)
                                                                                          tgcCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATAAATAACAGTGTCTGTCCGGTTAATATg                                               <  6:301305/135‑1 (MQ=255)
                                                                                             cGAGGGCCCGAACAGCAAACGTTTTTTTTTTTTCTGTAAAATTTCAAAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGata                                        <  5:246559/139‑1 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATgaa    >  1:142130/1‑137 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCa    >  8:80512/1‑139 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCa    >  8:76309/1‑139 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCa    >  6:58728/1‑139 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCa    >  1:89218/1‑139 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCa    >  3:56225/1‑139 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCa    >  4:196029/1‑139 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCa    >  4:252728/1‑139 (MQ=255)
                                                                                                                                 tCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAAAATGGTCAGGCACAGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGAAATATCTTTACTCCAGTGCTTTTTATTGATAAAATGc     >  3:133905/1‑138 (MQ=255)
                                                                                                                                   cAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCATa  <  7:80512/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                          
TCTCCTTGCATTGTCCAAAATCGAAAACGGCACATTCGGATACTGTGAAGAAACTGGGGCCCCCAT‑TCCTCTTGCAAAGCTGGCAGTGCTGCCGACGGCCCGAACAGCAAACGATTTTCTTTATTCTGTCCAATTTGAGAAAAAAACGCTTCCGATATGGAAATCAACGGATATCGAATATGGTCAGGCACTGTATGAATAACAGTGTCTGTCCGGTTAATATGCCTGATATTTCTTTACTCCAGTGCTTTTTATTGATAAAATGCATA  >  BS168_reference_genome/1616440‑1616708

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: ‑