Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,587,102 G→A 7.0% C211C (TGC→TGT ydeQ ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009131,587,1020G→A 7.0% 2.7 0.0 ‑0.258C211C (TGC→TGTydeQputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (25/26):  total (28/30)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.18e-01

TCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTA  >  NC_000913/1586814‑1587399
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                         
gaaatgtgtataagagacgatatgtacgctaaTGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGGCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGcc                                                                                                                                                                                                                                                                                                <  2:166633/270‑1 (MQ=255)
   caATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTTTTGCATCAGTTGCTGTGTTCGAAACCACTTGTCGGCTACTGTCAGTGGTGGCACCCGAACGATAAAATGAGAGCTTTGTTTTGCTTGAGCAATATACGCCAAgg                                                                                                                                                                                                                                                                                            >  1:412076/1‑301 (MQ=255)
     attttgcATAGATAAAATTGAACCCAATGACTGACTGGCCGGTGCCGGCTGAAACCTGTTTTCCCGTTTGTCCATAGGTTGCGCTTAACCCCAGAGTTACTTTGCTTTTATTTACTGGCCCCAGTGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAgggg                                                                                                                                                                                                                                                                                          <  2:377801/295‑1 (MQ=255)
               gATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAGAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTGCTgcg                                                                                                                                                                                                                                                                                >  1:91693/1‑301 (MQ=255)
               gATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTACCCAGTATTTTTCCATTACTCATTAGAGTGACACCCACTCCCGATCCTTTTGTTGCATCAGGTGCTGTGTTAGCTACCACTTCTCGCGTACTGTCGGTGTTGGACCCCGAAAGATAACATGAGATCTTTTGTTCGCTTGGTCATTATCCGCCCAGGGTTTTTTTTgc                                                                                                                                                                                                                                                                                 >  2:69329/1‑300 (MQ=255)
                    aaGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATTAACGcc                                                                                                                                                                                                                                                                                                >  1:166632/1‑280 (MQ=255)
                           acCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTTTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATGTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCAGTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAg                                                                                                                                                                                                                                                                    <  2:225603/301‑1 (MQ=255)
                           acCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAACCACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGAAGTTTCCGCCACGGGGTATTTATGCGCTTccctgaagg                                                                                                                                                                                                                                                                    >  1:34637/1‑294 (MQ=255)
                                     gactgGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGATGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTTTAGTATACTGTCAGTGGTTGCACCCGAAAGATAAAAAGAGCACTTTTGTTCGTGTGAGAAATATCAGCCACGGGGTAATTCTGCGCTTCGCGCAAAGTTCGGgagaa                                                                                                                                                                                                                                                          >  2:331478/1‑300 (MQ=255)
                                       ctgGACGGTTCTGGCTTAAACCGTGTTTCCGTGGTGTCCATAGGTTGCGCTTAAACCCAGAGTTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGGCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATcc                                                                                                                                                                                                                                                        <  2:85512/301‑1 (MQ=255)
                                              gTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGCTGCTTTTGTTGCATCAGGTGCTGTGTTAGAAACCCCTTGTCGGCTACTGTCAGTGGTGTCACCCGAAAGATAAGATGAGAGCTTTTGTTCGCTTTAGAACTATCCGCCCAGGGGTATTTCTGTGTTTTCAGGACAGTTCGAGGGTTACACGTTGa                                                                                                                                                                                                                                                 >  1:112470/1‑301 (MQ=255)
                                                        aatacTTGTTTTCTGGTTTGTCATAGGTGGCGCTTAAACCCAGTGGTACTTTGCTTTTGTTTACTGTCCCCAGGGAGACATGTTCTCCCGTAGCCAGTATTTTTGCATTACGCATTAGAGTGACACCCACTCCCTATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGAc                                                                                                                                                                                                                                        <  2:34637/296‑1 (MQ=255)
                                                                     tGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCAACTCCCGATGCTTTTGTTGCATCAGGTGATGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAACAATATACGCAACGGGGTATTTCTGCGCTTCCAGGAAAGTACGGGAGATCAACGGTGACATTTAGACTATCAACGGTGCaa                                                                                                                                                                                                                          >  2:280087/1‑300 (MQ=255)
                                                                         ttGTCCATAGGTTGCGCTTAAACCAAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTCTTTTTCCATTACGCATTGGAGGTACACCACCTCCCGGTGATTTTGTTGCATCAGGTGATGTGTGAGCAAACAATTGTCGGGTACTGTGAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTAAGCCATATACGCCACGGGGTATTTCGGCGCTTCCAGGAAAGACAGGAAGATACACGTTGCACGTTGGACTATCAACATTGCAgccg                                                                                                                                                                                                                       >  2:263642/1‑301 (MQ=255)
                                                                           gTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTACAGGAAATGCCGGGAGATCCACGGTGACATTTAGACTATCAACGGTgacgcctc                                                                                                                                                                                                                      >  1:259019/1‑293 (MQ=255)
                                                                                   gTTGCGCTTAAACCCAGAGGTACGTTGCTTTTATTTACTGTCCCCAGGGAGACATTGTCTCCCGTAGCTAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCGTTTGTTGCATCAGGTGCTGTGTTAGCCAACACTTGTCGGCTACTGTCAGTGTTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTgg                                                                                                                                                                                                               <  2:63244/299‑1 (MQ=255)
                                                                                           tAAACCCAGAGGTACTTCGCTTTTATTGACTTGCCCCAGGGAGACAGTTTCTCCCGTAGCCAGTATTTTTCGCTTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGCGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAacaaca                                                                                                                                                                                                     <  2:383241/301‑1 (MQ=255)
                                                                                              aCCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGAGCTTCCAGGAACGTCCGGGAGATCCACGGTGACCTTTCGACTATACACGGTGCAGCCGACTGTTGGCATACAAACATTg                                                                                                                                                                                                  >  2:86856/1‑301 (MQ=255)
                                                                                              aCCCAGAGGTACTTTGCTTTGATTTACTGTCCCCCGCGAGACATTTTCTCCCGTAGCCAGTAGTTTTCCAGTACGCATTAGAGTGACACCAACTCCCGATGCTGTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTg                                                                                                                                                                                                  <  2:91693/301‑1 (MQ=255)
                                                                                                      ttaatttGGTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGa                                                                                                                                                                                          <  1:392740/297‑1 (MQ=255)
                                                                                                      gTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTGGTCGGCTACTGTCAGTGGTGGCAGCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACACCAGGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATGTCGACTCTAAACGGGGCAGCCGCCTGTTGGCATACCAACATTGTTATTTGa                                                                                                                                                                                          >  2:289466/1‑301 (MQ=255)
                                                                                                                tttatttaAGGCACCCAAGGAGACAGTTTCGCCCGGAGCCAGTATTTTTCCATTACGCAGTAGAGTGACACCAACTCCCGATGCGTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTGTTGTGCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTc                                                                                                                                                                                <  1:189520/301‑1 (MQ=255)
                                                                                                                 ataagtgaaagccacAGGGAAGCATTTTTGCCGGTAGCCAGTATTTTTCCAATAAGCAGTAGTGTGACACACACTACTGATGCTTTTGTTGCATCAGGTGCTATGTTAGCAAACACTTGTTGGCTAATGTCAGTGGTGGCACCTGAAAGATAGAATGAGAGCTTTTGTACGCTTGAGCAATATACGCCAAGGGGTATTTATGCGCTTCCAGGAGAGACCGGGAGATCCACGGTGACATTTCTACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTc                                                                                                                                                                                <  2:200381/286‑1 (MQ=255)
                                                                                                                       aCTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCACCGGTTCAGCCGCCTGTTTGCATAACAACATTTTTATTTTAAATAATATTCCACTTaa                                                                                                                                                                         >  1:389053/1‑301 (MQ=255)
                                                                                                                       aCTGTCCCCAGGGAGACATTTTCTCCCGTAGACAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTTTTAGAAAACACTTGTCGGCTACTGTCAGTGGTGGCAACCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTACGCTTCCAGGAAAGTACGGGAGATCCACGGTGACCTTTCGACTATAGAGGGTGGAGCCGCCTGTCGGAATAACAACATTGTTATTTGAAATAATATTCCACGTaa                                                                                                                                                                         >  2:64564/1‑301 (MQ=255)
                                                                                                                          ttCCCCAGTGAGACATTTTCTCCCGTAGTCAGTATTTTTCCATTACTCATTATGTTGACACCCACTTCCGCTGCTTTTGGTGCATCAGTTGATGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTTGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCCGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAAt                                                                                                                                                                      <  2:389062/300‑1 (MQ=255)
                                                                                                                          gTCCCCAGGGAGACATTTTCTCCCGTAGCAAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTATGTTAGCAAAAACTTGTCGGCTACTGTCAGGGGTGGCACGCGAAAGATGAAATGAGAGCTTTTGGTCGCTTGAGCAAAATACGCCACGGGGTATTTCTGCCATTCCAGGAAAGTACGTGAAATCCACGGCGACATTTCGAATATCAAGGGTGCAGCCGCCTGTTCGCATAAACACATTCATAGTGGAAATAATATCACACGTCAAAt                                                                                                                                                                      >  2:280712/1‑301 (MQ=255)
                                                                                                                          ctCCCCAGGGCGACACTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCGTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTTACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAAt                                                                                                                                                                      <  2:259024/300‑1 (MQ=255)
                                                                                                                            ccccAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGCTAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGAGCTTCCAGGAAAGTCCGGGAGATACCCGGTGACATTTCGACTATCACCGTTCCAGCCGCCTGTTGTCATAGCAACATTGTTATTTGCAATAATATTCCACGTAAAATTa                                                                                                                                                                    >  2:371770/1‑301 (MQ=255)
                                                                                                                            ccccAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATACTATTc                                                                                                                                                                                >  1:200379/1‑289 (MQ=255)
                                                                                                                                               cccGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGCAAGATAAAATGAGAGCTGTTGTTCGCTTGAGCAATATACGCTCTGGGGTATTTCTGCGCTTCCCGGACAGTCCGGGAGATCCACGGTGACATTTCGACTAGCAACGGTGCAGCCGCCTGTTGGCATAACAACACTTTTATTTTAAATAATCTTCCCCGTAAACTTACGCGGATTTCCGCTCCCTa                                                                                                                                                 >  1:131727/1‑301 (MQ=255)
                                                                                                                                                        ttggagcggaaatgtttgtacgggacagcacACACTCCTGATGCTTTTGTGGCATTAGGTGCTGTGTTAGGAAAAACTTGTGGGTTACTGTCAGTGGTGGCAACCGAAAGATAAAATGAGAGCTTTTGTTGGCTTGAACAATATATGCCAAGGGGTATTTCGGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGGTATTTGAAATAATATTCCCCGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCaa                                                                                                                                        <  2:327717/273‑1 (MQ=255)
                                                                                                                                                           tttgttCCATTACGCCGTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCGGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAACAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAACTACGCGGATTTCCGCTTCCTAATGTTGCAAttt                                                                                                                                     <  1:280081/297‑1 (MQ=255)
                                                                                                                                                                ttcctgacgctgtattgTGACACCCACTCCCGCTGCTTTGTTTGCATCAGGTGCTGTGGTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATACAATGAGAGCTTTTGTCCGCTTGAGCAATATTCGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTTACATTTCGACTATCAACGGTGCATCCGCCTGTTGGCATAACAATATTGTTATTTGAAATAATATTCCCCTTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATAc                                                                                                                                <  2:394028/285‑1 (MQ=255)
                                                                                                                                                                         cattaggggtaacaccACTCCCGATGGTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCGACTGTAAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCGTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAAt                                                                                                                        <  1:371761/286‑1 (MQ=255)
                                                                                                                                                                         cATTAGAGTGACACCCACTCCCGATGCTGTTGTTGCATCAGGTGCTGTGTGAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGCGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAAt                                                                                                                        <  2:11456/300‑1 (MQ=255)
                                                                                                                                                                              tagTGACACCCACTCCCGATGCGTTGGTTGCATCAGGTGCGGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACggg                                                                                                                   <  1:86856/299‑1 (MQ=255)
                                                                                                                                                                                    cacCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAACAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACTGTGCAGCCGCCTGTTGGCAGAACAACATTTTTATTTGAACTAATATTCCACGTAAAATTAAGCTGAATTCCGCTTCCTCATTTTGcccatgttccttatacaactcttcccctgc                                                                                                             >  1:327709/1‑271 (MQ=255)
                                                                                                                                                                                     acCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTAtttca                                                                                                                                                                                          >  1:271304/1‑220 (MQ=255)
                                                                                                                                                                                     acCCACTCCCGATGCTTTTCGTGCACCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATTAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGa                                                                                                                                                                                          <  2:271310/222‑1 (MQ=255)
                                                                                                                                                                                          cTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCAGGGGGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATACTATTCCACGTACAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACCGGCAATGACTTccc                                                                                                       >  2:78951/1‑300 (MQ=255)
                                                                                                                                                                                                      ttgttgCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTTAGTGGTGGCACCCGAAAGATAACATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGATATTTGAAATCATATTACACGTAAAATTACTAGACTATCCGCTCACTAATGTTGCAATTTTATCCATATGACTACGATCAATGACTTACTCCGCTTTCATGa                                                                                           >  2:321914/1‑300 (MQ=255)
                                                                                                                                                                                                           gCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGCGATCCACGGTGACATTTCGACTATCGCCGGTGCAGCCGCCTGTTGGCATAACAACATTGTTCTTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCCCTTCCTAATGTTGCAATTTTATACATCTGAATACGGGCAATGAATTCCCCCGCTTTTATTACAACt                                                                                      >  1:159758/1‑300 (MQ=255)
                                                                                                                                                                                                             aTCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATCATATTCAACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGACTCCGGGCAATGACTTCCCCCGCTTTTATGACAACTCCa                                                                                   >  1:282641/1‑301 (MQ=255)
                                                                                                                                                                                                                  gTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGCAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGCAAGTCCGGGAGATCAACGGTGACATTTCGACTCGCAACGGTGCAGACGCCGCTTCGCATAACAACATTGTTATTTGAAATAATCTTCCACGTAAACTTACGCGTATTTCCTCTTCATAATGTTGCACTTTTATACATGTTCATCCGGGCCATCACTTCCCACCTTTTTTTTAGACCTCCAcccg                                                                               >  1:423505/1‑298 (MQ=255)
                                                                                                                                                                                                                  gTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCgg                                                                               <  1:289460/300‑1 (MQ=255)
                                                                                                                                                                                                                       gtgtTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGCCTTCCCCCGCTTTTATTACAACTCCACCGgccccc                                                                          >  1:303852/1‑300 (MQ=255)
                                                                                                                                                                                                                             caactccacgcctcgtgaCTGTCAGTGTCCGCTCCCGAAAGATACTACGCTTGCTTGTTTCCGCGTGAGAGATATACGCAAAGTGGTATTTCTGCTCTTCCAGTATACTCCGTTAGATCCACAGTGACATTTTGATTATCAACGGTGCAGCCCCCTTTTTGCATAACAACATTTCTATTTGATATAATTTTCCACGTAAAATTACGCGGATTTCCGCTCCTTACTTTTTCTATTTTATACATATTTATACGGGCAATGTCTTCCTCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTg                                                                   <  2:112470/284‑1 (MQ=255)
                                                                                                                                                                                                                                      tGTCGGTTACGGTCAGGGGTGGCACCCGAAAGATAAAATGAGGGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAatat                                                          <  1:78951/301‑1 (MQ=255)
                                                                                                                                                                                                                                      tGTCGGCGACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGTTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAatat                                                          <  2:303860/301‑1 (MQ=255)
                                                                                                                                                                                                                                           gCTACTGTCAGGTGGGGCAACCGAAAGATCAAATTAGGGCTTTTGTTCGCTTGAGCAATATACGACAAGTGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGAGCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAAAATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATagag                                                      <  1:331470/300‑1 (MQ=255)
                                                                                                                                                                                                                                               cTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGCCCGCCCCAACTGGGGTACTATAGAGTTTTa                                                 >  1:177544/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       tAAAATGATAGCGTTTGTTCGCTTGAGCAATATACGCCAAGGTGTATCTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTAt                         <  2:177545/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ttggTTCGCTTTAGAAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCaa            <  2:412085/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       tttggCTTGGGGAATATAAGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATAc         <  1:280706/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gtCGCTTGATCAATATACGCCAAGGGTTAGTGCTGGGCGTCCAGGAAAGGCTGGGATATCCACGGTGACATTTCGACTAACAACGGTGCAGCTGCCTGGTGGCATAACAACAGTGTTATTTGAAATAATATTCCACGTAAAATTACGCTGATTGCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACGCCACCGGCCGCCCCAACTGGTGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAAtcc         <  1:263637/299‑3 (MQ=255)
                                                                                                                                                                                                                                                                                          cGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGACGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTACACGTAAAATTACGCGGCTTTCCGCTTCCTCATGTTGCAATTTTATACATATGCATACGGGCAATGACTTCCACCGCTTTTATGACAACTCCACCGGCAGCCCCAACGGGGGTAAAGTAGAGTTTTAATGGTACTGGCTACGGTGTCGTATCCCCAATTTCCAATACGtt      >  2:225052/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                              tGATCAATTTACGTCAAGGGGTATTTCTGCGCTCCCAGGAAAGTCCGGGAGCTCCACGGTGACAGTTCGACTATCAACGTTGCATCCTCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCTGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAAGACGTGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTa  <  2:423514/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                         
TCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATTACGCATTAGAGTGACACCCACTCCCGATGCTTTTGTTGCATCAGGTGCTGTGTTAGCAAACACTTGTCGGCTACTGTCAGTGGTGGCACCCGAAAGATAAAATGAGAGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTA  >  NC_000913/1586814‑1587399

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ 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.