Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,595,353 T→C 9.9% T50A (ACG→GCG)  cshB ← ATP‑dependent RNA helicase; cold shock

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,595,3530TC 9.9% 108.0 / 4.6 42T50A (ACG→GCG) cshBATP‑dependent RNA helicase; cold shock
Reads supporting (aligned to +/- strand):  ref base T (20/16);  new base C (2/2);  total (23/19)
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.19e-01

TTCCTGATAAATCTGGTTGGCAAGCTCTCTTGTCGGAGCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGTTTCAATTCATAAAGTTCAAATTTCG  >  BS168_reference_genome/2595228‑2595490
                                                                                                                             |                                                                                                                                         
ttCCTGATAAATCTGGTTGGCAAGCTCTCTTGTCGGAGCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCg                                                                                                                              >  2:177669/1‑139 (MQ=255)
               gTTGGCAAGCTCTCTTGTCGGAGCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGctttct                                                                                                                  <  6:80487/136‑1 (MQ=255)
               gTTGGCAAGCTCTCTTGTCGGAGCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGctttct                                                                                                                  <  3:134130/136‑1 (MQ=255)
                    cAAGCTCTCTTGTCGGAGCCGTAATGACTTCCTGAACGTCATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTaa                                                                                                          <  4:33069/139‑1 (MQ=255)
                    cAAGCTCTCTTGTCGGAGCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTaa                                                                                                          <  1:177669/139‑1 (MQ=255)
                    cAAGCTCTCTTGTCGGAGCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTaa                                                                                                          <  4:232063/139‑1 (MQ=255)
                                     gCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGcc                                                                                         >  1:191469/1‑139 (MQ=255)
                                     gCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGcc                                                                                         >  7:244845/1‑139 (MQ=255)
                                     gCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGc                                                                                          >  7:245931/1‑138 (MQ=255)
                                                tACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGata                                                                              >  7:32435/1‑139 (MQ=255)
                                                   cTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTgg                                                                                 <  4:252167/133‑1 (MQ=255)
                                                   cTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTgg                                                                                 >  3:252167/1‑133 (MQ=255)
                                                                  aGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTaaaaa                                                            >  3:236416/1‑139 (MQ=255)
                                                                     aGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCa                                                         <  2:191469/139‑1 (MQ=255)
                                                                     aGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCa                                                         <  8:244845/139‑1 (MQ=255)
                                                                     aGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGAAATCAGTCCGTTCGTAAAAACCa                                                         <  8:245931/139‑1 (MQ=255)
                                                                     aGGAGCTATTTTGTTTAATAAAGGCAGTAAATTATCGTGTTGCTTTCCAGTACCCGTCTGAGACTGGCCGATGTCGCTTTCTTTTTTTAACACTCAATGAATCAGCCTTTTCTGTATATCAGTCGGTTAGTAAAAACCa                                                         <  4:263734/139‑1 (MQ=255)
                                                                       gatcttttttttttAATAAACGAAGTAAAAAATCGTGTGTATTTCCAGAACCAGGCTTACACTGGCCGTTTACGCTTTCTTTTTTTAACCCTGCAGGAATCAGCGTTTTCTGGATATCTGTGGGTTCGTAAAAACCa                                                         <  6:194024/131‑1 (MQ=255)
                                                                       gATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAg                                                       >  3:98644/1‑139 (MQ=255)
                                                                        aTCGATTTTTTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAgg                                                      >  8:142411/1‑139 (MQ=255)
                                                                        aTCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGCCTGAGACTGGCCGATGACGcttt                                                                                                                    <  1:133264/77‑1 (MQ=255)
                                                                        aTCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGCCTGAGACTGGCCGATGACGcttt                                                                                                                    <  1:129247/77‑1 (MQ=255)
                                                                        aTCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGCCTGAGACTGGCCGATGACGcttt                                                                                                                    >  2:129247/1‑77 (MQ=255)
                                                                        aTCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGCCTGAGACTGGCCGATGACGcttt                                                                                                                    >  2:133264/1‑77 (MQ=255)
                                                                         tCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGc                                                     >  3:11246/1‑139 (MQ=255)
                                                                         tCGAGTTTGTTTAATAAATTAAGTAAATAAGCGTGAGTCTTTCCAGTACCCGGCTGAGACTGGACGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCGTTTCTGGATATCAGGCGGTTCGTAAAAACCAAGGc                                                     >  1:91314/1‑139 (MQ=255)
                                                                           gATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCg                                                    <  5:93017/138‑1 (MQ=255)
                                                                            aTTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTGTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAGCCAGGGCGag                                                  >  8:295811/1‑138 (MQ=255)
                                                                                  tttAATAAAGGCAGTAAATAAGCGTGAGGCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACt                                             <  7:142411/138‑1 (MQ=255)
                                                                                    taataaAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGc                                           <  5:5189/138‑1 (MQ=255)
                                                                                      ataaAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCAt                                         <  6:242392/138‑1 (MQ=255)
                                                                                            gCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTataa                                   >  4:119707/1‑138 (MQ=255)
                                                                                                      aGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACAATGGCGATGTACTGCATCTATAATAAATGGttt                         >  8:169822/1‑138 (MQ=255)
                                                                                                          tGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGtt                          >  5:284278/1‑133 (MQ=255)
                                                                                                          tGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGtt                          <  6:284278/133‑1 (MQ=255)
                                                                                                          tGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGTTTCAAt                     >  8:79384/1‑138 (MQ=255)
                                                                                                            aGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGcttt                                                                                                                    >  6:185334/1‑41 (MQ=255)
                                                                                                            aGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGcttt                                                                                                                    <  5:185334/41‑1 (MQ=255)
                                                                                                            aGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGTTTCAATTCa                  >  5:269309/1‑139 (MQ=255)
                                                                                                            aGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGTTTCAATTCa                  >  1:189193/1‑139 (MQ=255)
                                                                                                                  tCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGTTTCAATTCATAAAGt            >  5:106953/1‑139 (MQ=255)
                                                                                                                            gTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGTTTCAATTCATAAAGTTCAAATTTCg  >  3:218450/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                                         
TTCCTGATAAATCTGGTTGGCAAGCTCTCTTGTCGGAGCCGTAATGACTACCTGAACGACATCCTTAGCAGGATCGATTTTGTTTAATAAAGGCAGTAAATAAGCGTGAGTCTTTCCAGTACCCGTCTGAGACTGGCCGATGACGCTTTCTTTTTTTAACACTGCAGGAATCAGCCTTTTCTGGATATCAGTCGGTTCGTAAAAACCAAGGCGATGTACTGCATCTATAATAAATGGTTTCAATTCATAAAGTTCAAATTTCG  >  BS168_reference_genome/2595228‑2595490

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