Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,856 A→G 5.6% pseudogene (1465/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009131,466,8560A→G 5.6% 4.6 ‑1.8 ‑0.0141pseudogene (1465/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  new base (6/3):  ref base (51/76):  total (62/79)
Fisher's exact test strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.03e-01

Only 100 of 142 total aligned reads displayed.

ATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAA  >  NC_000913/1466556‑1467154
                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                           
ataccggtgcacgtagACCAGGCGGTTGTGAACAACGAGGGCGAGCGCGACATCACCAACGGTGGCACCGGCCATTAGATGAACTGTGATGGAGCCATGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGTAAGTGCTGTCAATAACGGTAAGGTGCTGAAGGCCGGCGATCTGGCTGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGa                                                                                                                                                                                                                                                                                                             <  1:500869/286‑1 (MQ=14)
       tagACGGCGACCAGGCCGTAGATAACAACGAAGGCGAGAGCGGCATCACCAACGGTGGTACCGGCACGAAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGAAAGGTGATTCAGTACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                                                                                                                                                                                                                                                      <  1:1065844/299‑1 (MQ=255)
         gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGCCCGCGCACCGGTGGCTACCGAGGACACCATGACCGTCACCGCTCCCCAGTGCACCGGCAACCAGGTTTACGCCg                                                                                                                                                                                                                                                                                                    >  1:252485/1‑301 (MQ=11)
            ggCGAGCAGGGCGTGGTTAAAAACGAAGGCAGGAGCGCCATCCCCAACGGGGGCACCGGCACGCATATGAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAAGGAAATGGCTGGCAATAACGGGAAGGTGAGTCAGGACGGCGATCGCGATGTAAGCGGAGACGGTAACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAc                                                                                                                                                                                                                                                                                                  <  2:549959/300‑1 (MQ=14)
                        ttgacaaacaactagtggcatggcgcCATCACCAACGGTGGCACCGCGACTCAGATTAACGGTGATGACGCAACTGCAAACACCAACGGCAAAACCACCGTTAACGGCAAGTATTCCACCGGTACGGAAATTGCTGGCAATAACGTGAAGGTGATTCAGGCCGGCGATCTCGATGTCAGCGGCGGCGGTTACGGTATTGTTTTCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                                                                                                                                                                                                                                                     <  1:218266/278‑1 (MQ=14)
                          tgttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGCAAGGCGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCCCGGTATTGATATCAACGGAGACAGCGCAACGGGGGCAGACAAGGGAACCATGGCCGTCCCCGATCCGGAGTCCATCGGTATCCAGATTCACGCCGACCAGGCCATCGTCGAt                                                                                                                                                                                                                                                                                  >  1:293939/1‑301 (MQ=14)
                            ttaaacatgaagcgtaagtccccAGCACCAACTGTGTCACCTGCACGCAGAGTAACGGTGAAGACGTAACAGCCAACAACAACCTCAAAACAACCGGTGAACGACACGATTACACCGGTATGGGACTGGCTGGCAATAACAGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGACATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATACAGATTGACGGCGACCAGGCCATCGTCAATaa                                                                                                                                                                                                                                                                                <  1:970004/280‑1 (MQ=14)
                             aaacaacAAAGGCGAGGAAACCGACACCAACGGGGGCACCGGCACCAAGATTAACGGGGATGACGCAAAGGCAAACAACAGCGGCAAAACCACCGTTGACGGCAAGGATTACACCGGTACGGAAATTGACGGAAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGAGGCGGCGGTCACGGTATTGAGATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAAc                                                                                                                                                                                                                                                                               <  2:797092/300‑1 (MQ=11)
                              aaaaacGAGGGCGAGAGAGCAATCAACGAGGGTGGCACCGGCAGGCAGAATAACGGTGAGGAAGCCACGGGAAACAACAACGGCAAAACCAACGTGGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                                                                                                                                                                                                                                                              <  1:960977/298‑1 (MQ=17)
                                    gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCg                                                                                                                                                                                                                                                                        >  1:472410/1‑301 (MQ=255)
                                     aaGGGGAGGGCGCCCTGAACAACGGTGGCAGCGGCAAAAAGATGACAGGGGATGACCGCACGGCAAACAACAACGGCAAAACAAGCGTTGACGGAAAGGATTACAACGGTACGGAACTGGCTGGAAATAAAGGGAAGGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCGCGGTATTGATATCATCGGCGACAGCGAAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCg                                                                                                                                                                                                                                                                        <  2:1047952/300‑1 (MQ=14)
                                       ggCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATGGACGGCGACCAGCCCCTCGTCAATAACGAAGGCgaga                                                                                                                                                                                                                                                                     >  1:554723/1‑301 (MQ=255)
                                            agcgcccACCACCACCGGTGGCACCGGACAGCAGAGTAACGGTGATTACACCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGAAACGGTGGATAACAAGGGCACCATGACCGGCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                                                                                                                                                                                                                                                 <  2:533965/295‑1 (MQ=17)
                                                 ccagcgcaaaCGGTGGAACCGGAACGCAGATTAAGGTGGATGAGGCAACGGTAAAAAAGACGGGGAAAACCACCGTTGACGTCAAGGATTCCGCGGGGGCGGAATTGGATGGGAATAAGGGAAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                                                                                                                                                                                                                                                            <  2:389098/292‑1 (MQ=14)
                                                  caccacccgctgtgcgaCCGTCACACAGATTAAGGCTGAAGATGCCACTGAAACCTAAAACGTCAAATCCATCGTGGCGGGCAAGGATTCAACCGCTAATGAACTTGGCATCAATATCGGGAGTTGTATGCAGGACGGCTATCTGGCTGTCGGCGGCGGCTGTCACGGTATTGATACCACCGGCGACAGCGCAACGGTGTATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCAATCGGTATCCAGATTGACGGCGATCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                                                                                                                                                                                                                          <  2:752704/285‑1 (MQ=9)
                                                  cATCACAAACGGTGGCAGCGGCACGCTGATTAACGGGGATGAAGCCACGGGCAACAACAAGGGCAAAACCCACGTTAACGGCAAGGATTCGACCGGGACGGAAATTGCGGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAcc                                                                                                                                                                                                                                                          <  1:990378/301‑1 (MQ=17)
                                                   aTCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACCAAAATTGCTGGCAATAACGGGAATGTGATTCAGGACGGCGATCTGGCTGTCAGCGGCGGCGGTCCCGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAAGGCACCATGACCTTCACCGCTCCGGACCCCATCGGTATCCAGCTTGACGCCGACCAGCCAATCGCCCATAACCAACGcagaggcccccacccca                                                                                                                                                                                                                                                         >  1:369507/1‑285 (MQ=14)
                                                    tCAACAACGGTGGCACCGGCAACCAGATTAATGTTTATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAt                                                                                                                                                                                                                                                                                  <  1:806134/275‑1 (MQ=255)
                                                     caaccaccggtgaacatctaCGCAGACTAAAGGGGATGACGCAACCGCAAAAAATAACGCAAAAACCACCGGTGACGGCAAGGATTCGACCGGTACGGAAATTGCTGGAATTAACGGGAAGGTGAGTCAGGCCGGCGATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTAACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAt                                                                                                                                                                                                                                                       <  2:725652/282‑1 (MQ=17)
                                                     cACCAACGGTGGCACCGGCACGCAGATTAACAGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGAAATAACGGGAAGGTGATTCAGGACGGCGAGCTGGATGTCAGAGCCGGAGGTCACGGTCTTGATATCCCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACAGTCACCGATCCGGAGTCCATCGGTATACAGATTGACGGCGACATGTCAATCGTCAATAACGAAGGCCAGTGCACTATAACAaaa                                                                                                                                                                                                                                                       >  1:66016/1‑300 (MQ=14)
                                                       ccAACGGTGTGCCAGGTACGAGGATAACCGGTGAGGAGGCAACGGCAAACAACAACGGCAAAGCCAGCGGGTAGGGCAAGTATTGTACCGGTAAGGAAAGTGCTGGCAAAAACGGGAAGGTGAGTCAGGACGGCGGCAGTGACGGATGGGGATAAGGGAAAGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATgg                                                                                                                                                                                                                                                     <  2:319194/301‑1 (MQ=12)
                                                           cGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCCATAACGAAGGCGAGAGCCCTATCACCAATggcgg                                                                                                                                                                                                                                                  >  1:862181/1‑300 (MQ=255)
                                                                      gacatgaGATTAACGGAGATGCAGCCAAGGAAAAAAAAAACGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGGACGGAAATTGATGGTAATAACTGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTAACGGTATTGATCTCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTc                                                                                                                                                                                                                                       <  1:932052/294‑1 (MQ=17)
                                                                             gATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGCAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGAACaa                                                                                                                                                                                                                                >  1:1007699/1‑300 (MQ=255)
                                                                             gATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCCCCAATGGCGGCACCGGCACTCAGATCaa                                                                                                                                                                                                                                >  1:146996/1‑300 (MQ=255)
                                                                             gATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAAAGCACTATCACCAATGGCGGCACCGTCTCTCAGATCaa                                                                                                                                                                                                                                >  1:443059/1‑300 (MQ=255)
                                                                             gATTAACGGTGATGACGACACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGACGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCCCCAATGGCGGCACCGGCCCTCAGATCaa                                                                                                                                                                                                                                >  1:866243/1‑300 (MQ=255)
                                                                                      ggaggaagCCACGGGAAACAAACACGGCAAAACCACGGTTGACGGGAAGGGTGCCACCGGTACGGAAATTGACGGTAAAAACGGGAAAGGGATTCAGGACGGCGATCGGGATGTCAGCGGCGGCGGTCAGGGTATGGAGATCACCGGCGACAGCGAAACGGTGGATAACAAGGGCACCATGACCGTCACCGAGCCGGAGTCCATCGGGATCCAGATTGACGGCGACCAGGCCAGCGTCAAGAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAacga                                                                                                                                                                                                                       <  2:1057506/293‑1 (MQ=14)
                                                                                       gatgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGAGCACGGTATTGATAGCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACAGTCAACGATCCGGAGTACATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACAAATGGCGGCACCGGCACTCAGAGCAACAGTAacgacg                                                                                                                                                                                                                     >  1:335514/1‑301 (MQ=255)
                                                                                                   gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACAATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGCGGGAACCGGCACACAGACAAACGGAAacgcggcaaccccacaca                                                                                                                                                                                                         >  1:352483/1‑286 (MQ=17)
                                                                                                   gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCCTCACCGATCCGGAGTCCAACGGTATCCAGATTGACGGCACCCAGGCCAACGTCAATAACGAACTCGAGAGCACTATCACCAATGGCGGCACCGGCACCCAGATCAACGTTAACGACCCCCCcgccacaa                                                                                                                                                                                                         >  2:362962/1‑296 (MQ=11)
                                                                                                     aaaCAACAACGGGCAAACCACCGGTGACCGCAAGGATTCCCGCGGCACGGAAATGGCTTGCAATAACGGGAAGGTGATTGAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                                                                                                                                                                                                       <  2:443053/301‑1 (MQ=255)
                                                                                                     aaaAAACAACGCGAAAACCACGCGGTACGCCAAGGACTGCACGGGGCGGGGAACGGTTGGAAATAACGGGGAGGTGAGTCAGGCGGGGGGACTGGATGACAGCGGCGGGGTGAAGGTGAGTGATACCCCCAGCGATAGCGCAACGGGGGAGAACGAGGGGTCCAGGGTCGATAACGGACCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATGGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                                                                                                                                                                                                       <  2:411642/301‑1 (MQ=12)
                                                                                                        caacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACAATGACCGTCACCGATCCGGAGTCCATCGGTATCACGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGCAGGCACCGTCACTCAGCTCAACGGTAACGACCCCCCCGCCAATACCacg                                                                                                                                                                                                    >  1:918942/1‑299 (MQ=17)
                                                                                                                 cgaaaCCAGCGGGGGCGGGAAGGATACGCACGGTACTGAAATGGCTGGCAATAACGTGAAGGTGATTCAGGCCGGTGGTCTGGGTGTCAGCGGCGGCGGTCACGGTATGGATATCGCCGGCGACAGCGCAACTGTGGATAACAGGGGCACCATGACCGTCACCGTTCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAc                                                                                                                                                                                            <  2:352478/298‑1 (MQ=255)
                                                                                                                   acaatgtagggcgaggggaaGTATTCCACAGGTACGGAAAGAGCTGGTAGTAACGGGAAGGGTAGTCAGGCTGGAGAGCGGGGTGACAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAGCGCAACGGTGGATAAAAAGGGCACCATGACCGTAACCGATCCGGAGTCTATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc                                                                                                                                                                                           <  1:131016/282‑1 (MQ=17)
                                                                                                                   aaaCCACCGTTGACGACAAAGATTCCGCCGGCACGGAAAGCGACGGCAATAACGGGAAAGTGATTCAGGACGGCGGTCTGGATGTCAGCGGCGGCAGTCACGGTATTGATATCACCGGAGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGGGTACATCGGTATCCAGATTGACGGAGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAGCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCAc                                                                                                                                                                                         <  1:421627/301‑1 (MQ=17)
                                                                                                                       ccccgggggACGGCAAGGAGTTGACCGGTACTGAAATTGCTGGCAATAACGGGAAGGGGATTCAGGACGGAGGTCTGGATGTCCGGGGTGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCTTAACCGATCCGGAGTCCAGCGGTATCCAGATTGACGTCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGt                                                                                                                                                                                      <  2:938584/293‑1 (MQ=18)
                                                                                                                            gTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGACATCGTCAATCACGACGGCGAGAGCACTATCACCAATGGCTGCACCGGCACACAGACAAACGGTACCGACGCCACCGCGCATAACAGTGGAAAAAACACCCTTGaag                                                                                                                                                                                 >  2:844281/1‑299 (MQ=18)
                                                                                                                                 cGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCGTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGCCCAGGCCATCGTCAATACCGACGGCGAGAGCACTAACCCCCATCGCGGCACCGGCACACAGATCCACGGTAACGACGCCACCACCAATACCAGTGGAAACACCACTGTGGAGGGaaaa                                                                                                                                                                            >  2:824467/1‑301 (MQ=18)
                                                                                                                                       ggATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGACACACGGTCACGACGCCACCGCGAATTACCGTGGAAAAACAACTGTTGATGGAAAAGATTc                                                                                                                                                                       >  2:846274/1‑300 (MQ=255)
                                                                                                                                       ggATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCCGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGAAAAAACCACTGTTGATGGAAAAGGtta                                                                                                                                                                       >  1:458662/1‑299 (MQ=255)
                                                                                                                                                             agggggaagaaagggAAGGTGATTCAGGACGGTGAGCTGGAGGTGCGCGGCGGGGGTCAAGGTATTGATAGAACCGGGGAAAGAGAAACGGGGGATAAAAAGGGCACCATGACCGTAACCGGTCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTTAATAACGAAGGCGAGAGCACGATCACCAATGGCGGCACCGGCACTCAGATTAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGc                                                                                                                                                 <  2:474481/288‑1 (MQ=255)
                                                                                                                                                                ggCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCCCCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAATATCACCAATGGCGGAACCGGCCCTCAGCTCAACGGTAACGAAGCCACCGCGACTAAACGTGGAAAAACCAATGTTGATGGAAAAGATTACCCGGGTACCAAACGCGCGGGCaaaa                                                                                                                                             >  2:715688/1‑299 (MQ=18)
                                                                                                                                                                ggCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATACCGAAGGCGAGACCACTATCACCAATGGCGGCCCCGGCACTCAGATCAACGGTACCGACGCCACCGCGAATAACCGTGGAAAAACCAATCTTGATGGAAACGAATCCCCGGGTCCAAAAAGCGCgggagata                                                                                                                                             >  2:544395/1‑296 (MQ=18)
                                                                                                                                                                 aagtaagcGGAAAGGGGATCAGGAGGGAGAGATGGAGGTTAACGGCGGCGGTCACGGTATTGATATCAGCGACGAAAGCGCAACGGTGGAGAACAAGGGTACCATGACCGAAAACGATCCGGAGTACATCGGGAGCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAGCACCAATGGAGGCACCGGCACACAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAatat                                                                                                                                            <  1:64191/293‑1 (MQ=255)
                                                                                                                                                                  cAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACACGGGCACCATGACCGTCCCCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCCTCGTCAATAACGAAGGCGCGAGCACTATCACCAATGGCGGCACCGGCACTCAGCTCAACGGTACCGACGCCACCGCGAATAAc                                                                                                                                                                                                       >  1:513901/1‑241 (MQ=255)
                                                                                                                                                                  cAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGCGCACTATCACCa                                                                                                                                                                                                                                                         >  2:755194/1‑191 (MQ=255)
                                                                                                                                                                  aaaaaaCGGGAAGGTGATGCAGGACGGAGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATACCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATc                                                                                                                                           <  1:1025842/297‑1 (MQ=255)
                                                                                                                                                                              ggTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                                                                                                                                                                                                                                                 <  2:105644/171‑1 (MQ=255)
                                                                                                                                                                              aagttaTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAc                                                                                                                                                                                                                                                           <  2:54604/173‑1 (MQ=255)
                                                                                                                                                                                tGAGTCTGGACGCGGATGTGTATGGCAGCGGCGGCGGTAACGGGAATGAAATCACCGGCGACAGCGCAACAGGTGAGAACAAGGGAACCAGTACCGTCACCGAGCCGGAGTCAAGCGGGAGCCAGAGTGGCGGCGAACAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATc                                                                                                                              <  2:918914/300‑1 (MQ=255)
                                                                                                                                                                                   aTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAAGGTATTGATATCACCGGCGAAAGCGCAACGGTGGATAAAAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGTGAGAGCACGATCACCAATGGCGGCACCGGCACTCAGGTCAACGGTAACGACGCCACCGAGAATAACAGGGGAAAAAACACAGCAGATGGAACAGATAACCCGGGAACCAAGATAGCGGGCAACGACGGCACTGGaaaatcgg                                                                                                                           >  1:794824/1‑296 (MQ=17)
                                                                                                                                                                                             ggggtggggcggTCAGCGGCGGCGGTAAAGTTATGGATAGAACCGGCAACAGCGGAACGGGGGATAACAAGGGCAACGGGACCGTCACCGAGCCGGAGTCCAGCGGTATACAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGAGCAACGGTAACGACGCCACCGCGAAGAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGcc                                                                                                                  <  2:823135/289‑1 (MQ=255)
                                                                                                                                                                                               gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTACATCTGGATGGTCGACTGa                                                                                                               >  1:739735/1‑301 (MQ=255)
                                                                                                                                                                                                ctggtagttgaaggggcggcggTCACGGTATGGTTCTCAACGGGTAAAGCGAAACGGTGGATAAAAAAGGAGCCATGCCCGGCACGCATCCGGAGTTCATCGGTATCCAGATTGACGGCGACCAGGCCAGCTTCAATAACGAAGGCGAGAGCACTATAACCAATGGCGGCACCGGCCCTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGAc                                                                                                              <  2:866216/287‑1 (MQ=255)
                                                                                                                                                                                                    ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGCAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACGCAGATTAACGGTGATGACGCCACCGCAAACAACAACGGCAAAACCGCCCTTGACGGAAAGGATTCCACGGGTCCGAAAATTGCCGACAATACCGGCaggggaacacgagacgagcgccagacttat                                                                                                          >  2:1052571/1‑272 (MQ=14)
                                                                                                                                                                                                     gATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGGTCAACGGTAACGACGCCACCGCGAATACCAGTGGAAAACCCACTGTTGATCGAAAAGATTCCACGCGTACCAAAATCGCGGGCAGTATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTa                                                                                                         >  1:408025/1‑301 (MQ=255)
                                                                                                                                                                                                       tgggccagggcggggtgcacgtgaatGATATCAACGGGGACAGCGCAACGGGGGATAACAAGGGCACCAGGACCGTGAGCGATCCGGAGTCCATCAGGATACAGATTGACGGCTACCAGGCCATCGTCAAAAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTAc                                                                                                        <  1:839139/276‑1 (MQ=255)
                                                                                                                                                                                                       ggTCAGCGGCGGCGGTAACGTAATTGAAATCCCCGGCGACAGCGCAACGGTGGGTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTACCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGAGGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTAc                                                                                                        <  2:554712/299‑1 (MQ=255)
                                                                                                                                                                                                          aaGCGGCGGCGGGCCCGGTATTGATATCACCGGCGAAAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGc                                                                                                    <  1:844307/300‑1 (MQ=255)
                                                                                                                                                                                                                   cggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAACTCGCGGGCAATATCGGAATTGTaa                                                                                                                                 >  1:743134/1‑263 (MQ=255)
                                                                                                                                                                                                                    cgTCACGGTATTGATAAACCGGGCGACCGCGGAACGGGGGATAACAAGGGCACCATGACCGTAACCGAGCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAGGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATg                                                                                          <  2:739709/300‑1 (MQ=255)
                                                                                                                                                                                                                     gTCACGGGTTTGATATCACCGGGGCCAGCGCAACGGTGGATAACAAGGGCACCCTTACCGTCACCGTACCGGCGTCCATCGGTAGCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATg                                                                                          <  2:136776/300‑1 (MQ=255)
                                                                                                                                                                                                                      acACGGTATTGATATCAACGGCAACAGGGAAACGGTGGATAACAAGGGCATCATGACAGTAACCGAGCCGGAGCCCATCGGTATCCAGATGGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAATATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGAGAATAACAGTGGAAAAACCACTTTTGATGGAAAAGAGTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgt                                                                                        <  1:273643/300‑1 (MQ=255)
                                                                                                                                                                                                                       gccagggaTTGATATCAACGGCGACAGGGCAACGGGGGATAACAGGGGAACCATGACCGGCACCGACCCGAAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACGAGCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATg                                                                                          <  2:1019016/292‑1 (MQ=255)
                                                                                                                                                                                                                        aCGGTTATGATATCACAGGCGACAGGGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAGCCAGAGTGACGGCGAGCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgtgt                                                                                      <  2:335509/301‑1 (MQ=255)
                                                                                                                                                                                                                         cGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTCGTCAACAACGAAGGCGAGAGCACTATCACCAATGGCGACACCGGCACGCAGATCAAAGGTGACGACGCCAACGCAAACAACAACGGAAAAACCAATGTTGATGGAAAAGATTCCACGGGTACCAACACCGCGGGCCATATAGGCAGTGTAAATAAGGACAGTAGcaggcacgtcaccggcggcgcgcagggcgta                                                                                     >  2:214263/1‑271 (MQ=14)
                                                                                                                                                                                                                         cGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGAGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGAAAATCGGGCTGGTAACCTGAACGTTACAGGGGGTGCACATGGTGtt                                                                                     >  2:490204/1‑301 (MQ=255)
                                                                                                                                                                                                                          agaacgtGAGATCACCGGCGACAGCGCACCGGCGGATAACAAGGGCACCAAGACCGTCACCGATCCGGCGTACATCGGTATCCAGATTTACGGCGACCAGGCCATCGGCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCATATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgtgt                                                                                      <  2:987584/294‑1 (MQ=255)
                                                                                                                                                                                                                           ggtatgtatacaaCCGGCGAAAGCGCAAAGGTGAAGAACAAGGGCACCAGGACCGGAACCGATGCGGAGGCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgtgt                                                                                      <  1:544404/287‑1 (MQ=255)
                                                                                                                                                                                                                           gTATTGATATCACCGGCGACAGCGCACCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAAGCCGTCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCAACCGCGAATCACAGTGGAAAAACCACTGGTGATGGAAAAGATTCCAAGGGAAACAAAAt                                                                                                                                                        >  1:1015010/1‑232 (MQ=18)
                                                                                                                                                                                                                           gTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAATGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAc                                                                                                                                                                                                                                                           <  1:1014933/133‑1 (MQ=255)
                                                                                                                                                                                                                              tgttatagCACCGGGGAAAGAGCAACGGTGGATAACAAAGGCACCATGACCGTCACCGAGCCGGAGGCCATCGGTATCCAGATTGAATTCGACCAGGCCATCGTCAAGAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATTGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTgaga                                                                                 <  2:600228/297‑1 (MQ=255)
                                                                                                                                                                                                                              tctcataacagctgataCAGCGTAACGGCAGATAACAAGGGCACCATGACCGTCACCGATCCGTAGTCCAGTGGTATCCATATTGACGGCGACCATGCCATCGTCAATAACGAAGGCGAGAGCCCTAACACCAATGGCGGCACCTGCACGCAGATCAACGGGAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTgaga                                                                                 <  1:422621/285‑1 (MQ=25)
                                                                                                                                                                                                                                taatatccccggcacacccgcACCGGTGGAAAACAAGGGCACCCTGACCGTCACCGATCCGGAGACCATCCGTATCCAGATTGCCGTCGACCAGGCCATCGTCAATAACCAAGGCGATAGCACTACCACCAAGGTCGGCACCGGCACTCCGAGCCACGGTAACGACGCCACCGCAAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTGCCAAAATCGCGGGCAATATCGGCATTTTAAATCTGGATGGTAGCCTTACTGTTACAGGCGGTGCGCATGGTGTTGAGAAc                                                                               <  2:794798/283‑1 (MQ=255)
                                                                                                                                                                                                                                tGATATAACCGGCGCAAGCGCAACGGTGGATAACAAGGGCACCAGGAACGTAACCGATGCGGAGTCCATCGGTATCAAGATGGACGGCGACCAGGCCATTGGCAAGCACGAAGGCGAGAGCATTATAACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCTACCGCGAAGAACAGTTGAAAACCCACTGTCGATGGAAAAGATGCCACGGTTACCCAAATCGCGGGCAATATCGGCAATGTAAATCTGGATGGTAGCCTCACTGTTACAGGCGGTGCGCATGGTGTTgaga                                                                                 <  1:1104344/299‑1 (MQ=255)
                                                                                                                                                                                                                                 gATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAGCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAACCCACTGTTGGTGGACAAGATTCCCAGGTTACCCAAATCGCGGGCGAAAGAGGCATTGTGAATCAGGATGGGACCCTGACGGTAACAGGCGGGGCGCCTGGGGTTGAGaaaa                                                                              >  2:1004886/1‑299 (MQ=21)
                                                                                                                                                                                                                                   tatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATc                                                                                                                              >  1:594753/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tcaatcgcAACGGTGGATAACCAGGGCATCCTGACCGTTACCGATACGGAGTACATCGGTTGCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGATAGCACTATCACCAATGGCGGCACCGTCACTCAGATCAACTGTAACGACGCCACCGCGAATACCAGTGGAAAAACCACTTTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAAc                                                                   <  1:887777/295‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAACGGCAAAACCACTGTTGATGGCAAAGATTCCACCGGTACGAAAATCGCGGGCAATAACGGCAGTGTGAATCTGGACGGTAGCCTGACTGTTACAGGCGGTGACCATGGTGTTgagactgcatcttcacaccacacgacacgg                                                       >  2:152858/1‑274 (MQ=14)
                                                                                                                                                                                                                                                         gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATCACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCCTTGTCACTCTGGATGGTAGCATGACTGTTACAGGCGGTGCGCCTGGTGTTGAGAACATTGGTGACACCGGCACGGTTaaca                                                      >  1:150874/1‑301 (MQ=255)
                                                                                                                                                                                                                                                           ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAACCCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATAACGGCATGGTAGAAATGGATGGTAGACTGAAGGTTACAGGGGGGGCGCTGGTTGTAGAACACATTGTTGatacaggagagactagagaa                                                    >  2:604055/1‑282 (MQ=21)
                                                                                                                                                                                                                                                            gagaaCAAGGGCACCATGACCGTGACCGATCCGGAGTCAAGCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACAAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAACACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaacaac                                                    <  2:1020596/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTAt                                                                                                                                                                                                                                                              <  2:921160/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGACTAACAGTGGAACAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCTCGGGCAATATCGTCATTGTACATCTGGATGGTAGCATGACTGTTACAGGCGGTGCGCAGTGTGTTGAGAAAATTGGTGACAACGGCCCGGTTAACAACAAAGga                                              >  2:661999/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                        cccATGACCGTCAACGATCCTGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCAGCCCAAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCAACGGCAATCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttg                                       <  1:824493/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aCCATGACCGTCACCGGTACGGGGTACATCGGGAGCCAAAATGACGGCGACCAGGGCAGCGTCAATAAAGAAGGCGAGAGCAATAGCACCAATGGGGGCACGGGCACTCCGATCAACGGGAACGAAGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttg                                       <  2:1007669/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aCCATGACCCTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTTAATAACGAAGGCGAGAGCACGATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttg                                       <  1:117956/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ataaCCTTCACCGAGCCGGAGTCCATCGGAAACCAGATTGACGGCGACCTGGCCACCGTCAATAACGAAGGCGAGAGCAATATCACCAAGTGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACCGCGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttgttg                                    <  1:1052602/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aTGAACGTCACCGGCCCTGATGCAATGTGTATCCCTATTGACGGCTACCACGCCATCGTCAATAACGAAGGCGATAGCACGACCACCAAAGGCGGCACTGGCACTCAGATCAACGGTAACGACGCCTCCGCGAATAACAGTGGAAAAACCACTGTTGACTGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATTGTAGCCTTACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTTGTGACAACGGCACGGTTAACAACAAATGAGATAttgttg                                    <  1:214265/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGGATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATc                                                                                                                              >  1:1056692/1‑209 (MQ=255)
                                                                                                                                                                                                                                                                                     ccGCTCAGGAGTCTATCGGGAGCCAGATTGACGGCGAGCAAGACATCGTGAATAACGAAGGCGGGAGCACTAGCACCAAGGGCGGCACCGGCACGCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAGGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATAc                          <  2:293935/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       gATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                                                                                                                                                                                                                         >  2:236601/1‑76 (MQ=255)
                                                                                                                                                                                                                                                                                       gATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                                                                                                                                                                                                                         <  1:236603/76‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          cacgAGTCCATCGTGATCCATACTGACGGCGACCAGGGCAGAGTCTATAACGAAGGCGAGAGCACTATCAACAATGGCGGCACCGGCATTCAGATCAACTGTAACGACGCTCCCTCGAAGAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACTGGTACCAAAATCGCGGGCAATATCTGCCTTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGTTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGAt                     <  1:823812/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                tacacagatATCCAGAAGCACGGCGACACAGCGCACGTCAATAACGAAGGCGAGGGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACATGGGAAAAACCACTGTTGATGGAAAAGATGCCCCGGGTACCAAAATCGCGGGGAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGAtt                <  2:45142/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   ccagtgaTCCAGATGTACTGTGACCCGGCCATCGTCAATAGAGAAGGCGAGAGCACTATCACCAACGGCGGCCCAGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGGTGATGGAAACGATTCCACGGGTAACAAAATCGCTGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTGACATGCGGTGCGTATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACATCAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCg            <  2:186442/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCCCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGCTTCCACGGGTACCCAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGTTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTTCGAACATTGATGACCACGGCACGGTTAACACCAAAGTCGATCTTGTTGTTTCAGATACTGACTAGATTGGCGTGCtc       >  1:760149/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                              gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGCTGGCAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGCTGGTAGCCTGACTGTTACAGGCGGTGAGCATTGTTTTGAAACCATTAGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGTTACTGGATCGGTTGGctgtcacatcaa  >  2:1116610/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                           
ATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAA  >  NC_000913/1466556‑1467154

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