Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,025,078 A→C 19.2% T72P (ACG→CCG)  yhdH → putative sodium‑dependent transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,025,0780AC19.2% 83.5 / 2.2 47T72P (ACG→CCG) yhdHputative sodium‑dependent transporter
Reads supporting (aligned to +/- strand):  ref base A (17/19);  new base C (2/7);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGCGATTTGGAAGTTCCCTTACGTAGCGGGGACCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACAGTGTAATCGGGGGATGGATT  >  BS168_reference_genome/1024942‑1025200
                                                                                                                                         |                                                                                                                          
gggCGCGATTTGGAAGTTCCCTTACGTAGCGGGGACCAATGGGGGAGGAGCCGTTTTTCTTATCTTCGTGTTTTTTACCATCCTCTTAGGCGATCCGCTTTTGGTGGGGGAATTTTTTTTTGG‑TAGACGGAACCAGAc                                                                                                                           >  3:74908/1‑138 (MQ=255)
  gctcgATTTGGAAGTTCCCTTACGTCGCGGGGACCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTATTTAACATCCTCTTAGGCTATCCGCTTTTGGTTGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGaa                                                                                                                        <  8:263286/136‑1 (MQ=255)
   cgcgATTTGGAAGTTCCCTTACGTAGCGGGGACCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTTTTTACCATACTATTACGCTATCCGCGTTTGGTGGGGGTATTTATATTATG‑AAGGCGGAACCAGGCGa                                                                                                                         >  7:214815/1‑137 (MQ=255)
                  cccTTACGTAGCGGGGACCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCata                                                                                                         <  2:353113/138‑1 (MQ=255)
                      tACGTAGCGGGGACCATTGGGGGAGGAGCCTGTTTTCTTATCTTCGGGTTATTTACCATCCTCTTAGGCTATCCGCGTTTGGGGGGGGAATTTATTTTTTG‑AAGACGGAACCAGCCGAATGCCATTTATGCATATaaa                                                                                                     >  1:106615/1‑138 (MQ=255)
                             cGGGGACCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTGTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGCCGAATGCCATTGATGCATATAAAAAAGAAgc                                                                                             >  7:148303/1‑139 (MQ=255)
                             cGGGGACCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAgc                                                                                             >  5:423746/1‑139 (MQ=255)
                                  aCCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTAGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGcgccgc                                                                                        >  1:132958/1‑139 (MQ=255)
                                          gggAGGAGCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGCCTATCCGCTTTTGGTGGGGTAATTTATTTTTGC‑ACCACGGACCCAGACGAAGGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCt                                                                                <  3:426387/139‑1 (MQ=255)
                                           ggaggaGCCTTTTTTCTTATATTAGTGTTATTTACCATACTCTTAGGCTATACGCATTTGGTGGGTTAAGTTATTTTTTG‑AAGACTGAACCAGACTAATTCCATTGCTGCTTATAAAAAAGAATATCCTCGATCATCtt                                                                               >  6:425104/1‑139 (MQ=255)
                                            gaggagCCTTTTTTCTTATCTTCGTTTTATTAACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGCAATTTATTTTTGC‑AAGACGGAACCAGACTAATGCCATTGATCCAAATAAAAAAGACGCGCCGCGATCAGCTTg                                                                              <  4:310103/139‑1 (MQ=255)
                                            gaggagCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTg                                                                              >  3:201925/1‑139 (MQ=255)
                                            gaggagCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCGTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTg                                                                              >  5:82365/1‑139 (MQ=255)
                                              ggagCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTTGTGGGGGAATTTATTTTTTG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAATAAATAAGCGCCGCGCTAAGCGTGGt                                                                            >  1:225861/1‑139 (MQ=255)
                                              ggagCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCa                                                                                                                   <  5:115325/100‑1 (MQ=255)
                                              ggagCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCa                                                                                                                   >  6:115325/1‑100 (MQ=255)
                                                     tttttATTATCTTCGTGTTTTTTACCATCCTCTTAGGCTATCCGCTTTTGGGGGGGTATTTTATTTTTGT‑AAGACGCAACCAGACTAATGCCATTGATGAATATAAAAAACACGCGCCGCGATCAGCTTGGTTTCTTa                                                                      <  8:109903/138‑1 (MQ=255)
                                                       tttCTTAGATTGTTGTTATTAACCATCCTATTCGGCTACCCGGTTGGGGGGGGGGAATTATTTTTTGC‑ACCACGCAACCACCCGAATCCAATTTATCAAAAAAAAAAACACGCGCCGCGATCAGCTTGGTTTCTTACg                                                                    <  2:90920/138‑1 (MQ=255)
                                                            tATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTg                                                               >  7:109903/1‑138 (MQ=255)
                                                             aTCTTCGTGTTATTTACCATCCTCTTATGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGa                                                             >  3:35182/1‑139 (MQ=255)
                                                                tgcGTGTTATTCCCCACCCTCTCCGGCTCGCCGCTTTTGGTGGGGCAATTTATTTTGGG‑ACCCCGTCCCCACCCGAAGCCCATGGAGGCAAAAAAAACCGCCGCGCCGCGATCAGCTTGGTGTCTTACGGGGTGGAt                                                            <  2:17266/135‑1 (MQ=255)
                                                                 tCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AACCCGCAACCAGACTAATGCCATTGATGAATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCg                                                          <  6:82240/138‑1 (MQ=255)
                                                                  cGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGc                                                        >  7:169785/1‑139 (MQ=255)
                                                                      ttATTTACCCTCATATTAGGCTTGCGGCGTTGGGTGGGGTAATTTATTTTTGG‑ATGCCGGAACCACTCCATCGACTTTGTCGCATATAAAATAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAg                                                    <  8:360614/139‑1 (MQ=255)
                                                                         ataacaaTCCGCTTCGGCTCTCCGCTTGTGGTGGGCGAATTTATTTTTGG‑ACTCCGGAACCACACGAATGCCATTGATGCATATAAAACACACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAgcg                                                  <  6:171274/132‑1 (MQ=255)
                                                                           tcccATCCTCTTAGGCTAGCCGCTTTTGGGGGGGTAATTTATTTTGGC‑AAGACGCACCCAGCCGAAGGCCGGTGAGCCATATAAAAAACACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAgcggc                                                <  7:242241/136‑1 (MQ=255)
                                                                                tCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCgtgggt                                          >  3:230362/1‑136 (MQ=255)
                                                                                   tctTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCgtgggtgtg                                       >  3:73194/1‑133 (MQ=255)
                                                                                        ggCTAGCGGGTTTGGGGGGGGGATTTTATTTTTGC‑ACGACGCAACCCCCCGAATGCCATTGCTGCATATAAAAACGACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGc                                  <  8:144109/139‑1 (MQ=255)
                                                                                          ttttCCGGTTTTGGTGGGGTAATTTATTTTTTT‑ACTACGCAACCATAGAAATGCCATTGTTTAATATAAAATAGACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTg                                <  8:552/136‑1 (MQ=255)
                                                                                          cTATCGGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGCAACCAGACGAATGCCATTTATGCATATAAAAAAGACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTg                                <  6:82365/139‑1 (MQ=255)
                                                                                          cTATCGGCTTTTGGGGGGGCATTTTATTTTTGC‑AACCCGCAACCCGACGAATCCCTTTGAGGCATAAAAAAACGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTg                                <  3:426611/139‑1 (MQ=255)
                                                                                          cTATCCGCTTTTTGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTg                                <  2:132958/139‑1 (MQ=255)
                                                                                          cTATCCGCTTTTGGGGGGGCAATTTATTTTTGC‑AAGACGCAACCATCCCAATGCCATTGATGCATAAAAAAACGACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTg                                <  3:377463/139‑1 (MQ=255)
                                                                                           tATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGGTTGTTTctgctgg                               >  4:426611/1‑138 (MQ=255)
                                                                                             tCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTGTTAGTGGTTTCg                             >  6:300696/1‑139 (MQ=255)
                                                                                             tCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCg                             >  5:197953/1‑139 (MQ=255)
                                                                                             tCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTGTGTTTGTGGGGTCg                             >  4:377463/1‑139 (MQ=255)
                                                                                                   ttgggtggGGTATTTTTTTTTTGT‑ACCACCGAACCAGCCTAATGCCTTTGATCCATAAAAAAACCACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTa                        <  2:127830/135‑1 (MQ=255)
                                                                                                   ttgggtggGGTATTTTATTTTTGG‑AAGACGGAACCAGACTAATGCAATTTATGAATATAAAAAAGACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTa                        <  6:423746/135‑1 (MQ=255)
                                                                                                   ttggggggggtATTTTATTTTTGG‑AAGAGGGAACCTGACTAAGGCCATTTATGAATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTa                        <  8:122515/132‑1 (MQ=255)
                                                                                                     gggtggGGGAATTTATTTTTGG‑AAGCCGCACCCAGACGAAGGCCATTGATGCAAAAAAAAAACCCGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACa                      <  7:18767/137‑1 (MQ=255)
                                                                                                                  tatttTTGG‑ACGACGGACCCAGACGAATGCCATTGATGCATATAAAAAAGACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACAGTGTAATCggggg         <  2:19366/138‑1 (MQ=255)
                                                                                                                   ttttttGGCACGACGCAACCAGACTAATGCCATTGATCCATAAAAAAACCACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACAGTGTAATCggggg         <  8:169785/137‑1 (MQ=255)
                                                                                                                       ttGG‑AAGACGGAACCATCCGATGGCCATTGATGATTATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACAGTGTAATCGGGggatgga    <  3:199482/138‑1 (MQ=255)
                                                                                                                         gg‑ATGACGCAACCCTACTATTGCAATTTTTGCATATAAAAAAGACGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACAGTGTAATCGGGGGATGGAtt  <  2:225861/138‑1 (MQ=255)
                                                                                                                          gcaaGACGGAACCAGACGAATGCCATTGATGCATAAAAAAAACAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACAGTGTAATCGGGGGATGGAtt  <  4:201925/136‑1 (MQ=255)
                                                                                                                                         |                                                                                                                          
GGGCGCGATTTGGAAGTTCCCTTACGTAGCGGGGACCAATGGGGGAGGAGCCTTTTTTCTTATCTTCGTGTTATTTACCATCCTCTTAGGCTATCCGCTTTTGGTGGGGGAATTTATTTTTGG‑AAGACGGAACCAGACGAATGCCATTGATGCATATAAAAAAGAAGCGCCGCGATCAGCTTGGTTTCTTACGGGGTGGATCGGCGTAGCGGCGTGTTTCTTAGTGCTGTCGTTTTACAGTGTAATCGGGGGATGGATT  >  BS168_reference_genome/1024942‑1025200

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