Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,357,070	A→C	85.0%	L14R (CTT→CGT)	sapA ←	antimicrobial peptide transport ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	1,357,070	0	A→C	85.0%	18.0 ‑0.6 ‑0.2	74	L14R (CTT→CGT)	sapA	antimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  new base (22/39):  ref base (3/9):  total (26/48)
Fisher's exact test strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.93e-01

GTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCC  >  NC_000913/1356773‑1357361                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                    gTAAAACAGTCGGTGTGTTGAAACGGAACAACGCGAGGCAGGTGGAAGGGATAGGGCGCGCCGTCGTCGAGTACTCCCCAGGTTTCGGCAAGGGCCGGCATCAGGCGATAGGGATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTTACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGt                                                                                                                                                                                                                                                                                                   <  2:332350/300‑1 (MQ=255)  taaaacatacgggtgttgtgaaCGGAGAATCGCGACGCAGTTGTATGCGATAGTGCGCGCGGGGGTCGATTACTTCCCAGCTTTCGGAAAGGGCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCAGAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGt                                                                                                                                                                                                                                                                                                   <  1:61603/279‑1 (MQ=255)    aaCCCGTCGGTTTTTTGAAAAGGAATATCGCGGCGAAGGGGGAAGCGTGAGGGCGCGCCGTTGTCGAGTAAGTCCCAGCTTTCGGCAAGTGCCGGCATCGGGTGATAGGTATAGGGATCGACATCCAGCAGTCGATCGTAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCGTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCTCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAg                                                                                                                                                                                                                                                                                               <  1:689899/301‑1 (MQ=255)       cggTCGGTTTTTTGAAAAGGAAAATGGCGATGCAGGTGGAAGCGAGAGGCCGGGCCGTTGTAGAGTACTTCCCAGCCTTCCGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATAGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCaa                                                                                                                                                                                                                                                                                            <  2:622066/299‑1 (MQ=255)                  gtgagatgggtaaaagagacagggggggAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACGGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAag                                                                                                                                                                                                                                                                                  <  2:349101/274‑1 (MQ=255)                         ggAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATCAAACTGGGCGGCAAGGGTATCG‑CACTTCACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTTCCCGCTGCCGCCCTAGACAAAACCGCTTTCGCGGATATCAGCATGCGTGGGCGATTCAGTCGCGGCGATTTCCTTCCCCCTCCCACGTCCATCCATCACCacactacacgctca                                                                                                                                                                                                                                                                         >  1:661045/1‑288 (MQ=255)                               ttggggaggaGGTGGAAGCGATAGGTGGCGCCGTGGTGGAGTACTTCCCAGCGTGCGGCAAGGTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGAGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGAGATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCt                                                                                                                                                                                                                                                                    <  1:1100260/292‑1 (MQ=255)                               tctcgAAGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTAATTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCt                                                                                                                                                                                                                                                                    <  1:549174/298‑1 (MQ=255)                               tcggcgcgacgggggAAGCGAAAGGTGGCGCCGGTTGCGAGGACTTCCCAGCTGTCGGCAAGTTCCGGCATCAGGGGATAGGTATAGGGATCGACATCCAGCAGGCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACGGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCt                                                                                                                                                                                                                                                                    <  2:522145/288‑1 (MQ=255)                                        aGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACCAGTCCAGCAATCACCaaacg                                                                                                                                                                                                                                                                                  >  1:349106/1‑276 (MQ=255)                                          ttggaagggagtgttggtgcCGTTGTCGATTGGTTACCAGATTTCGGAAAGTTCGTGCAACAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGTCGTCAAGGGTATCGACAATTAACCCACTGCTCGCTTTTGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCataat                                                                                                                                                                                                                                                          <  2:932510/282‑1 (MQ=255)                                           ttgAAGCGAGAGGGCGCGCCGGTGTCGAGTACTTCCCAGCTTTCGGCAAGTTGCGGAATCAGGCGGTAGGTATAGGGACAGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCGGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAatat                                                                                                                                                                                                                                                        <  2:389019/299‑1 (MQ=255)                                                 taatgggtcggctcGTTGTCGAGGACTTCCCAGCTTTCGGCAAGGTCCGGCAGCAGGGGATAGGTATAGTGATCGACAACCAGCAGTCGATCATAAAACTGGGCGGCAAGGGGATCGACAATTAACCCACTGATCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAggg                                                                                                                                                                                                                                                  <  2:129733/288‑1 (MQ=255)                                                 ccatgggTCGCTCCGTGTTCGAGTACTCCACAGCTTGCGTCAAGTTCCGGCATCAGTCAATTGGCATAGGCATCGACACCCAGCAGTCTACCCTAAAACTGGGCGGATATGGTATCGACAATTAACCCACTGCTCGCTTTTGCTGGGTTAAAGGTTTTGCCCTGCCCGCTTACGTAATAGACAAAACCTCTGTCGCGGCTATCACCATGCGGGGGAGATTCAGTCGCGTCGATTGCCTGACCACTCACACGTCCAGCAATCCCCAAAAGAGACGATAATACCTGGCGCATAATATTAAgg                                                                                                                                                                                                                                                   <  1:509099/295‑1 (MQ=255)                                                  gttatgTCGGGCCGTTGGCGAGTACGTCCCAGCTTTCGGCAAGTGTCGGAAGCAGGGGATAGGTAGAGGCAGCGACATCCATCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCGTCTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAggg                                                                                                                                                                                                                                                  <  2:1269/295‑1 (MQ=255)                                                       gTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGCTCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTCAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATCTCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGCCCCATCACACGTACAGCAATCCCCAAAAGCGACGATAATACCTGGCCCATAACATTAAGGGAttttt                                                                                                                                                                                                                                            >  1:890950/1‑300 (MQ=255)                                                       gTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCCATCACACGTCCAGCAATCACCAAAAGAGACGATACTACCTGGCGGATAATATTAAGGGCTTTTa                                                                                                                                                                                                                                            >  1:820808/1‑301 (MQ=255)                                                        ttgagcCGTTGTCTAGTACTTCCCAGATGTCGGCAAGTTGCGGCATCAGGCGGTAGGTATAGTGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTa                                                                                                                                                                                                                                            <  1:509926/296‑1 (MQ=255)                                                                       gTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCGGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGACTGACCACTACCACGTCCAGCAATCACCAAAAGAGAAGATAATACCTGGCGCATAATATTAAGGGAGTTTATGTAACGCGACTAt                                                                                                                                                                                                                              >  1:885675/1‑299 (MQ=255)                                                                            tCCCAGCTTGCGGCAATTTCCGGCATAAGGCGATAGGTAGAGGGATCGACATCCAGAAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGCCAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTa                                                                                                                                                                                                                       <  1:668253/301‑1 (MQ=255)                                                                            cgcgagcagtcggaaagggcaggCATAAGGCAATAGGTAGAGGGATCGACCGACATAAGGAGAGCAGAAGCCTGGGCGGTAAGGGGATCGACAATTCACCCCCGGCACGCTGTGGATGGGTTAAAGGTGTTGCCTGGCCCGCTGAAGCAATAGACAAAGCCGCTGTCGCGGATATCAGCATGAGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGGCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTa                                                                                                                                                                                                                       <  1:350530/280‑1 (MQ=255)                                                                                ttctgtCGGCAAGTTCCGGAATCAGGGGAGTGATATAGGGAGCGAAATGCAGCAGTCGATCATACAACTGGGAGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGAAATAGACAAAACCGCTGTCGCGGATATCAGAATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAatat                                                                                                                                                                                                                   <  2:359274/296‑1 (MQ=255)                                                                                         taaGTTCCGGCAAAAGGCGATAGGTATAGGGATCGAGATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGAc                                                                                                                                                                                                           <  1:825758/299‑1 (MQ=255)                                                                                         aaaGTGCCGGGATGAGGAGATATGTAGAGGGAGCGACAGCCAGCAGTCGATCATAAAACGGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGAGATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACa                                                                                                                                                                                                          <  2:527518/300‑1 (MQ=255)                                                                                            gTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTTGAGGGATTAAAGGTGGTGACTTGCCCGCTGACGCAATAGCAAAAACCGCTGTCGCGGATATTAGCATTCGTCGGCAGTTACGCCGCGGCGCTTGCCTTACCACTCACCATTCCAGCACTCACCAACAGAGACGATAATAGCTGGCGCATAATACTACGGGATTTTATTTAGCGAGGCTCACTTACTCTTATTTAATTACAttt                                                                                                                                                                                                       >  2:949907/1‑301 (MQ=255)                                                                                                    tacaggagtagggTATAGGGATAGATAACAAGAAGTTGAGACTAAACGCGGGCGGGAAGGGTAACGAAAATTAACCCACGGCTCGTTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGATGCTATCTTACTAATATTTAATGACATTTGCCATTAc                                                                                                                                                                                               <  2:987054/290‑1 (MQ=255)                                                                                                                gTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATCCCTGGCGCATAATATTAAGGGATTTTATGTAACGAGGCTATCTTACTAATATTTAATGACATTTGACATTACCGTTTGTTTTCa                                                                                                                                                                                   >  1:947109/1‑301 (MQ=255)                                                                                                                gTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGCAGATTCAGGCGCGGCGCTTGCCTTACCACTCACACGTCCAGCAATCACCCAAAGAGCCGATAATACCTTGCGCATCCTATTAAGGGTTTTTATGTAAAGATGCTATCTTTCTAAGATTTAATGACATTTACCCTTACCGTTTTTGTTc                                                                                                                                                                                    >  2:24817/1‑300 (MQ=255)                                                                                                                  ataGGGAACGAAATCCAGCAGGCGATCATAAACCTGGGCGGCAAGGGTATCGACAATGAACCCGCTGCTCGGTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAgg                                                                                                                                                                                 <  2:143683/301‑1 (MQ=255)                                                                                                                     gggATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCCCACGTCCCGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTCCTAATATTTACTGACATTTGCCATTACCGTTTGTGTTCAGGGGt                                                                                                                                                                              >  1:234438/1‑301 (MQ=255)                                                                                                                     gggATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGgcg                                                                                                                                                                                                                                                                <  2:391779/219‑1 (MQ=255)                                                                                                                     gggATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGgcg                                                                                                                                                                                                                                                                >  1:391784/1‑219 (MQ=255)                                                                                                                      ggATCGAAATCCAGCAGTCGATCATAAAACTGGGCGGAAAGGGGATCGACAAGTAACCCACGGCTCGCTTTGGATGGGTGAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAgggg                                                                                                                                                                               <  2:548362/299‑1 (MQ=255)                                                                                                                       gATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTACACCACTAGTCGCTGTAGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCAGTCGCGGATATGAGCATGCGGGGGAGATTAAGGAGCGGCTCTTGCCTGCCCGCTCACAGATCAAGCAATAGCCACACGAACCGATAATACCTGCCCCCGAATATTCAGGGAGTTTATGTCAAGAGGATATCTAGCACAAATTTAATACCATTTTCACTTACTGTTTGTGTAACGGGGTCg                                                                                                                                                                            >  2:93004/1‑301 (MQ=255)                                                                                                                         tCGTCATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGCTTGCCTGCCCACTCACACGTCCAGCAATCACCAAAAGAGCAGATCATACCTGGCGCATAATATTACGGGCTTTTATGTAAAGAGGCAATCTTCCTAATATTTAATGACCTTTGCCATTACCGTTTTTTTTCAGGGGTAGTa                                                                                                                                                                          >  2:251189/1‑301 (MQ=255)                                                                                                                            aCATCCAGCAGTCGATCATAAAACGGGTCGGCAAGGGAATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGTTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGTATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATg                                                                                                                                                                       <  1:199932/301‑1 (MQ=255)                                                                                                                               tacagcatgagATCATAAAACTGGGGGGCAAGGGTAGCGACAATTAACCCACCGCACGCTTTGGATGGGTTAAAGGGGTTGACTTGCCCGCTGACGAAATAGACAAAACCGCTGTCGAGGATATCAGCATGAGGGGGAGAATCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAgg                                                                                                                                                                    <  2:890924/291‑1 (MQ=255)                                                                                                                               cccAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTGTGGATGGGGTAAAGGGGTTGACTTGCCCGCTGACGAAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGGATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAgg                                                                                                                                                                    <  2:820782/300‑1 (MQ=255)                                                                                                                                ccAGAAGTGGACCAGAAAACTGGTCGGCAAGGTGATCGACAATTAACCTACTGCTGGCTTCGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGAGTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAAGACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATg                                                                                                                                                                       <  2:868860/297‑1 (MQ=255)                                                                                                                                          aTCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGCTAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTTTGTTCAGGGGTCGTAATGAGGGCCGTGtggtt                                                                                                                                                         >  1:452584/1‑301 (MQ=255)                                                                                                                                          aTCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATAc                                                                                                                                                                                                                                                                      <  2:758751/192‑1 (MQ=255)                                                                                                                                          aTCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATAc                                                                                                                                                                                                                                                                      >  1:758777/1‑192 (MQ=255)                                                                                                                                           tCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATAGTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCTTAATGAGGGCCGTGtgtttg                                                                                                                                                        >  2:238224/1‑301 (MQ=255)                                                                                                                                           tCAGAAAACCGGGCGGCAAGGGTATCGACAAGTAACCCACTGTAGCCGTTGGAGGGGTTAAAGCTGTTGACTTGCCCGCTGACACAATAGACAAAACCGCTGTCGCGGATATCAGCATTCTGGGGATATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTCCCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGtggttg                                                                                                                                                        <  2:169896/301‑1 (MQ=255)                                                                                                                                                      ggggggaaGGGTATCGAAAATTAACCCACGGCTGGATTTGGATGGGTGAAAGGGGTTGACGTGCCCGCGGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAAc                                                                                                                                              <  2:234436/294‑1 (MQ=255)                                                                                                                                                                      acaagttaaccACTGCTCGGTGTGGATGGGTTAAAGGTGTTGACTTGGCCGCTGACTCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                                                              <  2:822628/291‑1 (MQ=255)                                                                                                                                                                             aCCCAGTGCTAGCTTTGGATGGGTGAAAGGGGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCATAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                                                      <  1:906085/301‑1 (MQ=255)                                                                                                                                                                                  cTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATAc                                                                                                                                                                                                                                                                      >  2:843302/1‑152 (MQ=255)                                                                                                                                                                                  cTGCTCGCTTTGGATGGGTTAAAGGTCTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGTGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATAc                                                                                                                                                                                                                                                                      <  1:843328/152‑1 (MQ=255)                                                                                                                                                                                         cTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAACGAGACGATAATACCTGGCGCATAATCTTAAGGGATTTTATGTAAAGATGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGTCCGTTTGGTTGTTCGTTGACCAGGTGACGTCCCCATCTTTTCAACATTCGTACCCgg                                                                                                          >  2:611828/1‑300 (MQ=255)                                                                                                                                                                                         cTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCgg                                                                                                          >  1:510291/1‑300 (MQ=255)                                                                                                                                                                                         cTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGCCGATTGCCTGACCACTCACACGTCCCGCCATCCCCAAAAGAGACGATAATACCTGGCGCGTAATATTAAGGGAGTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTTCCATTACCGTTTGTGTTCCGGGGTCGTAATGAGGGCCGTGTGTTTGGTCGTTGCACAGGTGACGTCGCCCTCTGTTCAACCTTCGTAc                                                                                                              >  2:236202/1‑297 (MQ=255)                                                                                                                                                                                                  ggTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGCTTGCCTGACCACTCACACGTCCAGCAATCACCCACAGCGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGGGGCTATCTTACTAATATTTAATTACATTTTCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCTTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCTTTATTtct                                                                                                   >  2:975490/1‑299 (MQ=255)                                                                                                                                                                                                  ggTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGCTAATACCTGGCGCATAATATTAAGGGGTTTTATGTAAAGAGGATATCTTACTAATAATTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTGATAAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCCTTTGTTCAACATTCATACCCGTGTttttgat                                                                                                 >  1:184045/1‑297 (MQ=255)                                                                                                                                                                                                  ggTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCg                                                                                                                                                                                             >  1:137189/1‑209 (MQ=255)                                                                                                                                                                                                  ggTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCg                                                                                                                                                                                             <  2:137188/209‑1 (MQ=255)                                                                                                                                                                                                    tttaaGGGGTTGAGTTGTCCGAGGACGCAATAGCGCAAACTGGTGTCGTGGAGATCAGAATGCGGGGCAGTTTCATGCGCGGCTATTGCCTGACCACTAACACGTCCATCACTAACCAAATGAGACGTTAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCGATCTTACTAATATTTAATGACATTTGCCATTACCGTTTTTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctat                                                                                               <  2:772975/298‑1 (MQ=255)                                                                                                                                                                                                       aaggggtgtacgtGCCCGCGGAGGCAATAGAAAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                            <  2:510284/289‑1 (MQ=255)                                                                                                                                                                                                        aGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCCGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCCTCTGTTCAACATTCGTACCCGTGATTTCTctataa                                                                                             >  1:56535/1‑299 (MQ=255)                                                                                                                                                                                                         ggTGGTGACTGGCCCGGTGACGTAAAAGACAAAACCGCTGTCGCGGATAGCAGCATGCGGGGGAAATTCAGGCGAGGCGATTGCCTGACCACTCACACGTCCAGCAATCAGCAAAAGAGACGATAATAACTGGCGCATAATATTAAGGGATTTTAGGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                            <  1:238226/299‑1 (MQ=255)                                                                                                                                                                                                          gttggggctggCCCGCTGACGAAATAGACAAAACCGCTGGCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACAAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGTGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                                                          <  2:56534/290‑1 (MQ=255)                                                                                                                                                                                                               cgcgtggccGCTGACGCAATAGAAAAAACCGCTGTCGGGGAGATCAGCAGGCGGGGGACATACAGGCGCGGAGATTGCCTGATCACTCACACGTCCAGCAATCACCAAAAGAGTCGATAATACCGGTCGCAGAATATTAAGGTATTGTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCTTTGAACAGGGGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCGATAACTAGAACTc                                                                                    <  2:719022/294‑1 (MQ=255)                                                                                                                                                                                                                         cTGACGGAATAGACAAAACCGCTGACGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCGTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaac                                                                             <  1:24818/298‑1 (MQ=255)                                                                                                                                                                                                                               cAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTTAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTCTAACTATAACTCCCATAACAACTTatca                                                                    >  2:704094/1‑298 (MQ=255)                                                                                                                                                                                                                                         aaacGCTGTAGCGGATAGCAGAAGGCGGGGGAGATTAAGGAGCGGCGATTGCCGGACCACTCACACGTCCAGCAATAACCAAAAGAGACGAGAATACCTGGCGCATAATATTAAGGGAGTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCagag                                                          <  1:975520/298‑1 (MQ=255)                                                                                                                                                                                                                                           ccGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTGGCCATCTGTTCAAAATTCGAACCCGTGATTTCTCTATAAATATAACTAACGGAACAACTTAGAGAGAGGAAGAGCg                                                        >  1:502882/1‑301 (MQ=255)                                                                                                                                                                                                                                                              caggctgtggATATTCAGGCGCGGCGATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGCCGATAATACTGGGCTCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAAGTAAATTTGCCATTACCTTTTGTGGTCAGGGGTCGTAATGAGGGCCGTGTTGTTCGTCTTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACGCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaaa                                        <  2:184044/290‑1 (MQ=255)                                                                                                                                                                                                                                                                       gaaaTTCAGGCGCGGCGAATGCCTTAACACGCACACGTCCAGAAATACCCAAAAGAGACGTTAATACCTGGCGCATAATATTAAGGGATTGTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCAGTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGc                             <  2:452578/297‑1 (MQ=255)                                                                                                                                                                                                                                                                               cgcgcGGCGATTGCCTTACCAATAACAAGTCAAGCAAGCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTTCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAg                    <  2:947079/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                     gcgATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                            >  2:644977/1‑223 (MQ=255)                                                                                                                                                                                                                                                                                     gcgATTGCCTGACCACTCACACGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                            <  1:644988/223‑1 (MQ=255)                                                                                                                                                                                                                                                                                         tttccTGACAACTAAAACGTCAAGAAAACACCAAAAGAGAAGATAATACATGGCGTATAATATTAAGGGAGTTGATGTAAAGAGGCTATCTGAGTAATATTTAATTACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCGCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCAt          <  2:885649/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                         atGCATGAACACTCAGACGTACAGCAATCACCAAAAGAGCAGATAATCCCTGGCGCATAATATTAAGGGATTTGATGTCAAGAGGCTATCTTAATAATATTGAATGACATTTGCCATGACCGTGTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTTAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTc            <  1:949937/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 ccACTCACAAGTCCATCAATCACCAAAAGGAACGATAATACCTGGCCCATAATATTAAGGGATTTTAGTTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCTTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGcc  <  1:1064862/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                    GTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCC  >  NC_000913/1356773‑1357361

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