Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,466,276	T→C	9.9%	pseudogene (885/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_000913	1,466,276	0	T→C	9.9%	3.4 0.0 0.0	52	pseudogene (885/2513 nt)	ydbA	pseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (21/21):  total (25/27)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGA  >  NC_000913/1465984‑1466566                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                    tggaggggaggagagaaaaaggagcAAAGGAGTGGAAAAATAAACAAAACGGGGGGACTGTGATTACTGGGGATAATGCCACAGTAAACAACAAGGGCAAAACCACCGGTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGc                                                                                                                                                                                                                                                                                             <  2:132076/278‑1 (MQ=11)     agggggtggggcAGATGGGGAAACAAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCAAAATAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAAGAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAt                                                                                                                                                                                                                                                                                          <  1:270538/290‑1 (MQ=255)                acaaTGGGGCAGCGAAGTGGAAAAATAAACACAACGTCGGGCTTGATATTACGGGGGAGAAAGACACAGTAAACAACAATGGCAAAACCAGCGTTGACGGTAAAGATGCCACCGGTACGGAAATCAACGGGAATAACGGGAAGGTGATGCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGAGATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAGAGGTATACAGATTGACGGCGACCAGGCCATCGTCAATAACg                                                                                                                                                                                                                                                                               <  2:238601/297‑1 (MQ=14)                                  ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGCCAATCTGGATGTCAACGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACACGGGCACCATGACCTTCACCGAACCGACGTCCATGGTTATCCCGATCGACGGGAACAAGGCCATCTTCACTAACGCAGACGAAGGCACCAGCa                                                                                                                                                                                                                                                             >  1:288355/1‑300 (MQ=18)                                      atatAATCAAAAAGGTGAGCTTGTTATTACGGGTGATAAAGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCAACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAc                                                                                                                                                                                                                                                        <  1:587733/301‑1 (MQ=255)                                       tataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGT‑AACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGACATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAACTGGAAGACCGCGGCGGCGGTCACGGTATTGATATCACCGGTGGCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTAGACGATCCAGAGCACATGGGTCAGCAGAGCGACGGAGCCAAGGCCAGCGTGACACACGCAGGGAAGAGCACCGTCACCAACgg                                                                                                                                                                                                                                                      >  2:336201/1‑301 (MQ=255)                                        ttaagtaaaaCGGTGAGATTGTTATTACGGTGGATAAAGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAGTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACgg                                                                                                                                                                                                                                                      <  2:509097/293‑1 (MQ=255)                                              aaaaCGGGGAGCTGGTTATTACGGGGGAAAAAGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAAGAAACGGTAATAACGGGAAAGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcac                                                                                                                                                                                                                                                <  2:616198/299‑1 (MQ=255)                                                 aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGcac                                                                                                                                                                                                                                                                                           >  1:400538/1‑253 (MQ=255)                                                                  tcGGGGGAAAATGCAAAAGTAAACACAAAGGGGAATACAACCGTTGACCGGAAAGATGCAACCGGTGCGGAAATCAAGGGTAAAAACGGGAAAGTGATTCAGGACGGCGATCTGGAGGTGAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                                                                                                            <  2:426026/300‑1 (MQ=17)                                                                            acggaatgtcaaaCTACGATCGCAAAACCAACGGCAACGGAAATGATACCAACTGGACGGAAATCGGAGGGAAGGACGAGAGACATATGCAGGCCGGCGCTCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAGCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAc                                                                                                                                                                                                                  <  1:193290/291‑1 (MQ=14)                                                                               aaacacgaaaCAACAAGTGAAAAAACAACTTGGACGTGAAGGTTTCTACCGGTACGGAAATTAACGTGAAGAACGCGTAATGGCGCCAGGAGTGCGATCGGTGTGTCTGGTGCGGCTGTCATGTGATCGATATCACCGGTGAAAGCTCGACGGTTGATAACAATGGCACCATTACCGTCACCGATCCGGAGTCCATTGGTATCCAGATCGATGGTGACAAGGCCATCGTCAATAACGAAGGCTAGAGCATCATCACCAACGGTGGCACTGGCACCCAGATTAATGGTGATGACGCCACGGc                                                                                                                                                                                                               <  1:663264/294‑1 (MQ=255)                                                                                  caGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAGGTGATTCAGGACTGTCACTTAGACACAACTGACGCCGCCGACGAAAGCTTAGTCATAGGTGACAGCGCAACGCTTGATAACAAGGGCACCATGACCGTCACGCATACGGAGTCCATGGGGATCCAGATCGACGGTGACGAGGCCCTCGTCACTAACGAAGGCGAGAGCACCATCACCAACGGTGGAACCGGCCAGCAGAGTACAGGTGAGGCGACCACGGCaca                                                                                                                                                                                                            >  2:140347/1‑299 (MQ=9)                                                                                    gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGCCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGCTCCGGAGTCCATGGGTATACAGATCGACGGTGACAAGGCCATCGTCACTCACGAAGGCGAGAGCACCATCACCCACGGTGGCCCCGGCACCCAGATTAATGGTGATCACGCCACGGCAaaca                                                                                                                                                                                                          >  1:195059/1‑301 (MQ=255)                                                                                     aaaaAAAAAAGGGAAAAACCAAAGGGGACGGTAAAAACTCAACCGGGACGGAAATCAACGGTAATAACGGGAAAGGGATTCAGGACGGAGATCTGGATGTCAGCGGCGGAGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCaa                                                                                                                                                                                                             <  2:628904/296‑1 (MQ=17)                                                                                      aaaAAACAAAGGCAAAAGCACCGTTGAGGGAAAAGATTCCACCGGTACGGAAAGCAAGGGTAATAACGGGAAAGTGATTCAGGACGGCGGTCTGGATGTCAGCGGGGGCGGTAACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaac                                                                                                                                                                                                        <  2:195052/301‑1 (MQ=255)                                                                                         caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCCTCGTCAATAACGAAGGCGACAGCACCATCCCCACCGGTGGCACCGGCACCCAGATTACTGGagctgaacacccggcaacaaaaaaa                                                                                                                                                                                                     >  2:32191/1‑282 (MQ=17)                                                                                         caaatagaacacacccaccGGTGACGGGAAAGCTGTCAGCGGTACGGAAATCAACGTTAATAACGGGTAATTGATTCAGGACGGCGAGCTTGAGGTCGCCGGCGGCCGTCACGGTATTTATAGCACCGGTGACAGGGCGACGGTGGATAACAAGTGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                                                                                                     <  2:516320/287‑1 (MQ=17)                                                                                                         ctcgtttgcGGTAAAGACTCAACCGGGAAGCAAGTCAACGGGAATAACGGCAAAGTGATTCAGGACGGCGAACTGGATGGCAGCGGCGGCGGTCATGGTATTGATATCACCGGTGACAGCGGGACGGTGGATAAAAAGGGCACCATGACCGTCACCGATCCGGAGGCCAGGGGTATCCAGATCGACGGTGCCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTg                                                                                                                                                                                    <  1:73220/293‑1 (MQ=17)                                                                                                                aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGTGACCAGGCAATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACCCAGCTTAACTGTGATGACGCCACGGCAAACAACAACGGCAAAACCCCC‑GTTGACGGCaa                                                                                                                                                                             >  1:383368/1‑301 (MQ=255)                                                                                                                   gTAAAAATTCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACACGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTggtgatggtg                                                                                                                                                                           <  1:32193/300‑10 (MQ=255)                                                                                                                        gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCGGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGCTGACGCCACGGCAAACAACCACGCAAAAACCACC‑GTTGAcagacgtccccaaa                                                                                                                                                                     >  1:67330/1‑288 (MQ=255)                                                                                                                        gATCCCACCGGTGCGGAGAGCAAAGGTAAACACGGGAACGTGAGTAAGGAAGGCGATCTGGATGTCAGCGGGGGCGGGAACGGCATTGAGATCACCGGTGACAGCGCGACGGTGGATACCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACgg                                                                                                                                                                                <  2:67326/290‑1 (MQ=17)                                                                                                                              cccGGTACGCAAATCAACTGTAATAACGGGAAAGTGAGTCAGGACGGCGATCTGGATGTCAGCGGCGGCTGGCACGGTATTGGTATCACCGGAGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCAATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAACATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTa                                                                                                                                                               <  1:186186/300‑1 (MQ=255)                                                                                                                                   gacttaaaTCAACTGTCATCACCGGCAAGCGATTCCGGACGGCCTACGGGATGTCAGCTGCGGTTGTCAGGTAACGGATATCCCCGTTGACAGCGCTACGGTGGCTACCAAGCGCACCATGACCGTCACCCATCCGGACTCTATGCGTAACCCGATCGCCTGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTCGCACCGGCACCCAGATTAATGGTGATGCCGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCTAGGATTCCACCGGTACGGa                                                                                                                                                           <  1:336207/295‑1 (MQ=17)                                                                                                                                                ggTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACAGtaata                                                                                                                                             >  2:112581/1‑301 (MQ=255)                                                                                                                                                 gTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAACAGGCGTCCATGGGTATCCAGATCGACGGTGACACGCCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACACAGATTAATGGTGAGTACGCCACGGCAAACAACAACGGCCAAACCACC‑GTTCACGGCACGGATTCCACCGGAACGTACATCAACCGtaataa                                                                                                                                            >  1:635111/1‑301 (MQ=17)                                                                                                                                                                  atgcAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAgg                                                                                                                           <  1:112586/298‑1 (MQ=255)                                                                                                                                                                            gCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTAACGACAAGGATTCCACCGGTACGGAAATAAACGGTACAAACGAAAAAGATATTAAGAACgggtatcg                                                                                                                 >  1:720653/1‑295 (MQ=17)                                                                                                                                                                               aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGCAGCCACGGCCACCAACACCGGCAAAACCACC‑GTTAACGGCAAGGATTCCCCCGGACCGGAAAACAACGGTAAAAACAGAAAAGTTATTCCGGCCGGcaaacggga                                                                                                              >  2:741975/1‑293 (MQ=17)                                                                                                                                                                                  aggAGGTCAGCGGCGGCGGTCACGGTATTGATAACAGCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAAGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                                                                                            <  1:223114/299‑1 (MQ=255)                                                                                                                                                                                                        cGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCAACCAGATTAATGGTGATGACGCCACGGCAACACACAAAGGCAAACCACCC‑GTTGACGGAAATGCTTCCCCCGGTACGGAAATCACAGGTAAACACGAAAACGTGATTCCGGACGGCGACCTGAATGTCGGCGGCGGCGGTACCGGTa                                                                                       >  2:474083/1‑299 (MQ=17)                                                                                                                                                                                                          gTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGTGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGACGGCGCACGCACCGTCACCAACGGTGGCACCGGCCCCCCGATGAATGGTGAGACCGCCACGGCAAACAACAACGCCACACCCACC‑GTTTCATGCACGTGTTCACCCGGCACTGAACACAACGGTACTCCCCGAAAACTTATTCCTGACACCACTCTGTTTGGCCGCGGCCGCGGTcgctttatttg                                                                                   >  2:501379/1‑291 (MQ=255)                                                                                                                                                                                                          gTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGATAATAACGGAAAAGTTATTCAGGACGGCCATCTGGATGTCAGCGGCGGCGGTCACGGTCTTgg                                                                                   >  1:613197/1‑300 (MQ=255)                                                                                                                                                                                                            atggATATCACAGAGTACAGCGCGCAGGTGGAAAACAAGGGCACCCTGCACGACACCGATCCGGAGTCCATGGCTATGCAGCACGAAGGGGAAAAGGCCATCGTAAATACCGAAGGCGACAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATTGTTATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GGTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                 <  1:742002/298‑1 (MQ=17)                                                                                                                                                                                                            atagaaagcacccgggaaacggcgATGGTGGATAACAAGGTGACCAGGACAGTCATCGATCCGGAGGCCATGGGTATCCAGATCGACGGTGACAAGGTCATCGTCAAGAACGAAGGCGAGAGCAACACCACCAACGGTCGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGCCAAAACCACC‑GTGGACGGCAAGGATTCCACCGGTACGTAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                 <  1:474097/280‑1 (MQ=17)                                                                                                                                                                                                            aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACCAGGCCTTCCTCTTTTACGCATCTGACACCACCATCACCAGCGTAGGCACCGGCACCCAGGTTACTCGTGATGACGCCACGGCAACCAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTCCGCAAATAACCGGTAATAACGGCACAGTTCTTCAGCACGCTGCTATGGATGTCATCGCCGGCGCACACAGTCTTGcta                                                                                 >  2:208891/1‑301 (MQ=14)                                                                                                                                                                                                              tGAAATCACCGGTGAACTGGCGACGTGTGAGAACAAGGGCACCATGACCGTCATCTATCAGAAGTCCATGGGTATCAAGATTGACGGTGACAAGGCCATGGTCAATAACGAAGGCAAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                               <  1:509167/301‑1 (MQ=17)                                                                                                                                                                                                               gATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCCCCGATTCGGCGTCCGTCGCTGTCCACATCGACGGTGACATGGCTCTCGTCGATCGCGCAGTCACGAGCACCATCACCAACAGGTGCACCGGCACCCAGATTGATGGTGAGGCCGCCCCGGCCAACGACAAAGGCAAAACAACC‑GTAGACGGACACGATTCCACCGGTCCGGACATCAACGATAATACCGGCAAAGTTCTTCAGCGAGCAGAGCAGGAGGTCAGCGGCGGCGGCAACTGTATTGatacc                                                                               >  2:577337/1‑298 (MQ=255)                                                                                                                                                                                                                atacCACCGGGGACAGCGGGAAGGTGGAGAACAAGGTCACCAGGAACGTACACGATCCGGAGGCCATGGGTATCGAGAAAGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGTCAAAACCACC‑GTTGACGGCAAGGATTCCAACGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAc                                                                             <  1:501402/301‑1 (MQ=17)                                                                                                                                                                                                                                       gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCGAACAACACCGGAAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTCCGACAATAGACGGCAATACCAGAAAAGTTAATCAGGACGGCGACCTGGACGTCAGCGGCGGCGGACACGGTACTgcagtcacaggctagaccgacagggtggcta                                                      >  1:672610/1‑271 (MQ=14)                                                                                                                                                                                                                                       gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGAGGTGACGCCTAGTCATATAACAGCGGC‑CAACCACCAGTTGGAGGAAAGAAGCCCACAGGAACGGAAATCAACAGGACAAACGAAAAAGTGATACAGGACGGCGAGATGGACGGCACAAGAGGCGGAAACGGTAGTGATATCGAGGGCGACAGCGCGACGGTGGaga                                                      >  1:257601/1‑299 (MQ=9)                                                                                                                                                                                                                                                                        ttgtggtatatggagaccaagatatgaggTGACAAGGACATCGCAAAGAAAGAAGGCGTGTGTCACATCAACAACGGCGGCACCGGCACACAGAGTACTGTGGATGACGCCACCGGAAAAAACAACGGCAAAACCACC‑GTTGACGGGAAGGATGCCACCGGTACGGAAATCAACGGTAATAACGGAAAATTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                     <  1:393205/274‑1 (MQ=14)                                                                                                                                                                                                                                                                             gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGGTGACGCTACGGCATACAACAACGGCAAAACCACC‑GTTGACGGCAACGATTCCAGCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGAACGGCGGTCTGGAAGTCAGCGGCGGCGGACACGGTATACATAACACCGGCGACAACGCCAGGGTGGATAACAGGGGCACAACGActggcgcaaacccggaggccct                 >  1:396560/1‑279 (MQ=255)                                                                                                                                                                                                                                                                             gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACCCAGATTAATGGTGATGACGCCACCGCAAATAACAACGCCAAAACCCCC‑GTTGACGACCCGGAGTCACTCAGTACGGACAACAACGCTACTAACTGAAAGGTTATTCAGAACGGCGATCTGGATGTCAGCGGAGGCGGTCACGGTATTGATATCCCCGGCGACAGCGCCACGGTGGAAAACAAGGGCACCATACCCGTCACCGATCCGGAggcgat                 >  1:396555/1‑295 (MQ=255)                                                                                                                                                                                                                                                                                  tGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGACAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAACGTCACCGATCCGGAGTctggtgctta           >  1:577101/1‑292 (MQ=255)                                                                                                                                                                                                                                                                                  tGGGTATCAAAACCGACGGGGACAAGGCCAGCGTAAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGAACCCGGATTAATGGTGAGGACACCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                     <  2:577078/291‑1 (MQ=255)                                                                                                                                                                                                                                                                                    ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATCACACGGGCACCATGACCGCCACCGAGCCGGAGTCCCTCGGTATc          >  2:176631/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                     ggaTCCAGAGCGAAGGGGACAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGAACCCAGATTAATGGGGATGACGCAACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc         <  1:176638/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGGACGGAAATCACCGGTAATAACGGAAAAGTTATTCACGAAGGCGAGCTGGATGTCAGCGGCGGCGTTCACGGTATTGATATCACCGGCGACCACGCAACGGTGGATAACCAGGGCCCCATGCCCAGCACCACACCGCAGTCCCTCGGTATCCAgtatga  >  1:388295/1‑295 (MQ=17)                                                                                                                                                                                                                                                                                                    gTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg                                                                                                               >  2:584194/1‑183 (MQ=255)                                                                                                                                                                                                                                                                                                    gTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg                                                                                                               <  1:584218/183‑1 (MQ=255)                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                    TTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACC‑GTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGA  >  NC_000913/1465984‑1466566

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