Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,323,881 C→T 100% P27L (CCG→CTG)  xkdD → phage PBSX; conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,323,8810CT100.0% 141.8 / NA 41P27L (CCG→CTG) xkdDphage PBSX; conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/18);  total (23/18)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCGGGCAAAGCGAACAGAAAAAAACATTGAACAATGGTATAAGGACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACCGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCAATGCCTCTTTGATGGCCGCTGCGC  >  BS168_reference_genome/1323748‑1324015
                                                                                                                                     |                                                                                                                                      
tCGGGCAAAGCGAACAGAAAAAAACATTGAACAATGGTATAAGGACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATc                                                                                                                                    <  5:284244/138‑1 (MQ=255)
    gCAAAGCGAACAGAAAAAAACATTGAACAATGGTATAAGGACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTatat                                                                                                                               <  1:128495/139‑1 (MQ=255)
                 aaaaaaaCATTGAACAATGGTATAAGGACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAATCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATc                                                                                                                                    <  3:138582/121‑1 (MQ=255)
                 aaaaaaaCATTGAACAATGGTATAAGGACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATc                                                                                                                                    >  4:138582/1‑121 (MQ=255)
                                       taAGGACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGAt                                                                                            >  3:233783/1‑139 (MQ=255)
                                          ggACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTa                                                                                                                   >  8:71092/1‑113 (MQ=255)
                                          ggACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTa                                                                                                                   <  7:71092/113‑1 (MQ=255)
                                          ggACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTc                                                                                         >  3:254055/1‑139 (MQ=255)
                                          ggACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTc                                                                                         >  7:294165/1‑139 (MQ=255)
                                           gccgacgGTAAGTGAAAGTGTTGCCAAAGAAAAAACCGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCGGCCCCGGGCGGGGATACGAAACTTATATATATTTATTAAGCTTAAAGTGAAGAAGAAAAAGAAAGATTTcc                                                                                        <  6:208699/137‑1 (MQ=255)
                                           gacgagGTGAATTGAAAGTGTTGGCAAAGACAAACCAGGCAGATAAACGCCCTGCGCGGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACTGTGAAGAAGAAAAAGAAAGATTTcc                                                                                        <  6:204845/139‑1 (MQ=255)
                                              gacgGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTa                                                                                     >  3:199736/1‑139 (MQ=255)
                                              gacgGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTa                                                                                     >  3:51679/1‑139 (MQ=255)
                                              gacgGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTa                                                                                     >  4:164637/1‑139 (MQ=255)
                                               acgGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGCATTCCGTAc                                                                                    >  5:111965/1‑139 (MQ=255)
                                                     aGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAAc                                                                              >  3:45699/1‑139 (MQ=255)
                                                                  gCAAAGACAAAACACGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGGGCGGGGATACGAAAGTGATCTATATTTATTCCGCTTACCGGGAAGAACAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTg                                                                 <  6:92893/139‑1 (MQ=255)
                                                                   cAAATACAAAACAGGCAGAGATAAGCCCTGCGCCGGGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCTGCTTACAGTGAAGAATAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGc                                                                <  2:100333/139‑1 (MQ=255)
                                                                   aaaaGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACCTTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGc                                                                <  2:104904/138‑1 (MQ=255)
                                                                     aaGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTTTTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAg                                                              <  4:304125/139‑1 (MQ=255)
                                                                     aaGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAg                                                              <  8:164858/139‑1 (MQ=255)
                                                                                  cAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCt                                                 <  1:65127/139‑1 (MQ=255)
                                                                                  cAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCt                                                 >  6:255151/1‑139 (MQ=255)
                                                                                  cAGAGAAAAGCCCTGCGCCGTGGCGGGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAATAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCt                                                 <  4:67373/139‑1 (MQ=255)
                                                                                     agaAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAg                                              <  4:45699/139‑1 (MQ=255)
                                                                                               tGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAACCGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCg                                    >  1:243298/1‑139 (MQ=255)
                                                                                                 cgcCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGc                                   >  5:24791/1‑138 (MQ=255)
                                                                                                 cgcCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGc                                   >  7:308424/1‑138 (MQ=255)
                                                                                                      tGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTAtatttat                                                                                                                           >  4:269221/1‑45 (MQ=255)
                                                                                                      tGGCGTGCTGTCCCGTGCGGGGATACGAAACTGATCTAtatttat                                                                                                                           <  3:269221/45‑1 (MQ=255)
                                                                                                           tGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAATAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACATCAGCTCTTATTCGCAAAAAAGCAAt                        >  4:138293/1‑139 (MQ=255)
                                                                                                           tGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCACACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCGCGGTTTCGCAAAAAAGCAAt                        >  7:24666/1‑139 (MQ=255)
                                                                                                           tGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCAAt                        >  6:46532/1‑139 (MQ=255)
                                                                                                           tGCTGTCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCAAt                        >  7:24664/1‑139 (MQ=255)
                                                                                                                tCCCGTGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGCCCTCAGCTCTTATTCGCAAAAAAGCAATGCctc                   >  1:133558/1‑139 (MQ=255)
                                                                                                                     tGCGGGGATACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGC‑AAAAAGCAATGCCTCTTTGAt             >  2:224327/1‑139 (MQ=255)
                                                                                                                            aTACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCAATGCCTCTTTGAt             <  6:179941/133‑1 (MQ=255)
                                                                                                                            aTACGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCAATGCCTCTTTGAt             >  5:179941/1‑133 (MQ=255)
                                                                                                                               cGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCaa                         >  5:34746/1‑118 (MQ=255)
                                                                                                                               cGAAACTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCaa                         <  6:34746/118‑1 (MQ=255)
                                                                                                                                 aaaCTGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCAATGCCTCTTTGATGGCCGCTgcgc  <  4:199736/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TCGGGCAAAGCGAACAGAAAAAAACATTGAACAATGGTATAAGGACGACGGGAAGTGAAAGTGTTGGCAAAGACAAAACAGGCAGAGAAAAGCCCTGCGCCGTGGCGTGCTGTCCCGTGCGGGGATACGAAACCGATCTATATTTATTCAGCTTACAGTGAAGAAGAAAAAGAAAGATTTCCGTACTCAAACGGGCGGCTGATTGCAGCTGTATTTGACCTCAGCTCTTATTCGCAAAAAAGCAATGCCTCTTTGATGGCCGCTGCGC  >  BS168_reference_genome/1323748‑1324015

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