Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 914,845 G→A 6.4% L5L (CTG→TTG)  ybjE ← putative transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_000913914,8450G→A 6.4% 4.3 0.0 ‑0.065L5L (CTG→TTG) ybjEputative transporter
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (28/32):  total (30/35)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.03e-01

CGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGT  >  NC_000913/914558‑915133
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                
cgACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAAc                                                                                                                                                                                                                                                                                        <  2:479429/298‑1 (MQ=255)
         tttttCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTTGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAAccc                                                                                                                                                                                                                                                                              <  2:145987/299‑1 (MQ=255)
                    gATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAg                                                                                                                                                                                                                                                                 >  2:420956/1‑301 (MQ=255)
                        gtggtTGCGCCACGGAAGGCCTCGCTCCAGCCACATCAGGGCGGCAAGATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAAGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAAGCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATg                                                                                                                                                                                                                                                              <  1:44977/300‑1 (MQ=255)
                                         ggCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATg                                                                                                                                                                                                                                            >  1:276155/1‑301 (MQ=255)
                                                cTCCAGCAACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGacacac                                                                                                                                                                                                                                      <  2:231695/300‑1 (MQ=255)
                                                         aaTCAGGGCGGAAATATTACACAGTAAAATAAAGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTAGTGAGGTTATAGAGAAACGCCAGACTGATACCCATAAAAAAGAGACTAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGATTGTTGGCGAAGCGGAATGAAGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTACTTGAAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGc                                                                                                                                                                                                                            <  1:428550/300‑1 (MQ=255)
                                                                 cggcAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGa                                                                                                                                                                                                                     <  2:336765/300‑1 (MQ=255)
                                                                 cggcAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGa                                                                                                                                                                                                                     <  2:890815/300‑1 (MQ=255)
                                                                      atatTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAACGATGAAGTCGAAATGCGCGATGACACACTACTg                                                                                                                                                                                                                                 >  2:705245/1‑283 (MQ=255)
                                                                      atatTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTg                                                                                                                                                                                                                                 <  1:705230/283‑1 (MQ=255)
                                                                      atatTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTa                                                                                                                                                                                                               >  1:331545/1‑301 (MQ=255)
                                                                      atatTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTa                                                                                                                                                                                                               >  1:28962/1‑301 (MQ=255)
                                                                       aaataaccaGTAAAATAACGGTAATACTGACGGCAGAATAAGGCAGGATCGCCAACAGGTTACGGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTaa                                                                                                                                                                                                              <  1:691248/294‑1 (MQ=255)
                                                                                gTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATACCTTTTAACGCAACTTGTTGGAGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCAAGAAAACATGAAAACCCAGCCCTTGCAAAGATGAAGTAGAAATGCGCGATAACACACTACTGAAAGCGGAAGGCAGAGTAAAGTTGCAAt                                                                                                                                                                                                     >  1:18963/1‑301 (MQ=255)
                                                                                gTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAAt                                                                                                                                                                                                     >  2:285314/1‑301 (MQ=255)
                                                                                gTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAAt                                                                                                                                                                                                     >  2:364695/1‑301 (MQ=255)
                                                                                gTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCCCGATGACACACTACTGAAAGCGGAAGGAAGAGTAAAGTTGCaaa                                                                                                                                                                                                     >  1:352393/1‑300 (MQ=255)
                                                                                gTAAAATAACGGTAATACTGACGGCAAAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCaaa                                                                                                                                                                                                     >  1:168330/1‑300 (MQ=255)
                                                                                         cGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGaa                                                                                                                                                                                            >  2:481277/1‑301 (MQ=255)
                                                                                         cGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGaa                                                                                                                                                                                            >  2:279501/1‑301 (MQ=255)
                                                                                                  tGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGc                                                                                                                                                                                   <  1:420952/301‑1 (MQ=255)
                                                                                                    aCGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCAt                                                                                                                                                                                 <  2:18963/301‑1 (MQ=255)
                                                                                                     cGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGc                                                                                                                                                                                                                            >  2:344610/1‑257 (MQ=255)
                                                                                                           aataatGCAGAATCGCCAACAGGTTACTTGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGa                                                                                                                                                                            <  1:285311/299‑1 (MQ=255)
                                                                                                           aataatGAAGAATCGCCAACAGGTTACTGGCGATGTTATCGAGCACCGCCAGACTGAGACCCATAAAAAAGAGAATAAGGTGAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTAGCTAAGCGGAATGGTGTAACCCACAATCAGGGGAACCAGAATGAATAACAGCCCAGAAAACATGAAAACCCCGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGc                                                                                                                                                                          <  2:747354/301‑1 (MQ=255)
                                                                                                              aatGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGt                                                                                                                                                                       <  2:331548/301‑1 (MQ=255)
                                                                                                                     aaTCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTagag                                                                                                                                                                   <  2:168332/298‑1 (MQ=255)
                                                                                                                     aaTCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTagag                                                                                                                                                                   <  2:276158/298‑1 (MQ=255)
                                                                                                                     aaTCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTAt                                                                                                                                                                <  1:481273/301‑1 (MQ=255)
                                                                                                                     aaTCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACAAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTAt                                                                                                                                                                <  1:394060/301‑1 (MQ=255)
                                                                                                                     aaTCGCCAACAGGTTACTGACGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTAt                                                                                                                                                                <  2:352396/301‑1 (MQ=255)
                                                                                                                     aaTCGCCAACAGGGTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTagag                                                                                                                                                                   <  1:364691/298‑1 (MQ=255)
                                                                                                                                 gTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTagag                                                                                                                                                                   >  2:847417/1‑286 (MQ=255)
                                                                                                                                 gTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTagag                                                                                                                                                                   <  1:847391/286‑1 (MQ=255)
                                                                                                                                  ttACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTa                                                                                                                                                   >  2:576172/1‑301 (MQ=255)
                                                                                                                                       ggCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCACAGAAAACATGAAAACCCAGTCCTTGCAAAGATAAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGACATGTTATGCATAAGGAGCAGTAGAGTATTCGTGTTCATTTAACGaa                                                                                                                                              >  2:861372/1‑300 (MQ=255)
                                                                                                                                                      aaaCGCCAGACTGATAACCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGGAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTaa                                                                                                                               <  1:342766/301‑1 (MQ=255)
                                                                                                                                                                cTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGGAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTCCAAACtt                                                                                                                     >  1:548167/1‑300 (MQ=255)
                                                                                                                                                                        cATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATg                                                                                                             <  1:576166/301‑1 (MQ=255)
                                                                                                                                                                                   gaATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGtt                                                                                                                                                           <  1:83753/244‑1 (MQ=255)
                                                                                                                                                                                   gaATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGtt                                                                                                                                                           >  2:83754/1‑244 (MQ=255)
                                                                                                                                                                                                 aTCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAgg                                                                                     >  2:556504/1‑300 (MQ=255)
                                                                                                                                                                                                      gCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTTGCGTCAAAATCGCCGAGGCGTTc                                                                                >  2:494923/1‑300 (MQ=255)
                                                                                                                                                                                                                  aTTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAgg                                                                                                                                                                             <  2:144860/195‑1 (MQ=255)
                                                                                                                                                                                                                  aTTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAgg                                                                                                                                                                             >  1:144858/1‑195 (MQ=255)
                                                                                                                                                                                                                   ttAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAATATCGCCGAGGCGTTCCCTGACGCCCgggg                                                                  >  1:128218/1‑301 (MQ=255)
                                                                                                                                                                                                                    taataaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGc                                                                 <  1:556499/301‑1 (MQ=255)
                                                                                                                                                                                                                    taaaaaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAAGTAACCCACAATCAAGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAAACCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGc                                                                 <  1:861344/301‑1 (MQ=255)
                                                                                                                                                                                                                          cTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGccc                                                            >  1:416109/1‑300 (MQ=255)
                                                                                                                                                                                                                                aCGGAGCCTGTTGGCGAAGCTGAAAGAGGTAACCCAACATCAGGGGAATCAGAATGATTAACACTCCAGAACACATGAAAACCCAGGCCTTGCAAAGATGAAGTCGAAAGGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGGAATTAAAAGGAAATGTTGTGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGGGTCAAGATCGCCGAGGCGTTCCCTGAAGGCCTGGGCAGCCCACAGGGa                                                     <  1:742189/301‑1 (MQ=255)
                                                                                                                                                                                                                                                       aTGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTAc                              <  2:416113/301‑1 (MQ=255)
                                                                                                                                                                                                                                                         gAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                            >  2:607423/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                cccACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGtt                     >  2:247301/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                cccACAATCAGGGGAACCAGAATGATTAACAACCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCATATGGa                                                     >  1:326133/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                cccACAATCAGGGGAACCAGAATGATTAACAACCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCATATGGa                                                     <  2:326136/269‑1 (MQ=255)
                                                                                                                                                                                                                                                                cccACAATCAGGGGAACCAGAATGATTAACAACCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGa                                                     <  2:326127/269‑1 (MQ=255)
                                                                                                                                                                                                                                                                cccACAATCAGGGGAACCAGAATGATTAACAACCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGa                                                     >  1:326124/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                  aacaATCAGGGGAACCAGAATGATTAACAGCCCAGAAATCATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAc                   <  1:279498/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aTCAGGGGAACAAGAATGATTAAAAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTcg               <  1:494918/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTagag                                                                                                                                                                   <  1:440187/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTagag                                                                                                                                                                   >  2:440191/1‑146 (MQ=255)
                                                                                                                                                                                                                                                                               aCCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTaa                                                                                                                                                                                                              <  1:634562/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aCCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTaa                                                                                                                                                                                                              >  2:634573/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                                   gAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCATGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAATATTCGTTTTCATTTAACGATATTCTTGCGCTTTAATTACAAACTGCACCGATTTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAATGCCGGGGCCGCCCACATGACTGTGGCCTGATGCTGCTTTTTACAGGGATGTTACCTCTCGCCCGACGCgtg  >  1:451424/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                
CGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCAACAGGTTACTGGCGAGGTTATCGAGAAACGCCAGACTGATACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACAGCCCAGAAAACATGAAAACCCAGTCCTTGCAAAGATGAAGTCGAAATGCGCGATGACACACTACTGAAAGCGGAAGGACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGT  >  NC_000913/914558‑915133

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.