Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,479,642	T→C	G151G (GGT→GGC)	unknown →	Flp pilus assembly protein RcpC/CpaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,479,642	0	T	C	100.0%	188.2 / NA	66	G151G (GGT→GGC)	unknown	Flp pilus assembly protein RcpC/CpaB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/31);  total (35/31)

TGATTATATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTAGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAACGATCCTAGATTCTTTACATCCAGCCTAGATCCGAATCAAGAAGTTGCATATCGCATCTATGTTAAATCTGGTGATGATTACATTCTAAATACTATTTTAAGTGGAGACTATGTTGCTATCTATAGCCAACAAATTGCTTCAGATTCGCGCAACACTTATGATCGTAAAGATTTAGTTAAAATTTTAGGTAAAGTATTGGTTTTACAAGCAAGAAACTTTACTCAATCAACAGAACAAGGTGCTCCTGAAGCTGA  >  1/1479402‑1479824                                                                                                                                                                                                                                                 |                                                                                                                                                                                       tgatTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACa                                                                                                                                                                              >  1:736366/1‑251 (MQ=255)   atTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACAtt                                                                                                                                                                            >  1:348346/1‑251 (MQ=255)         ttttAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAAt                                                                                                                                                                      >  1:568385/1‑251 (MQ=255)                    aTTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAg                                                                                                                                                           >  1:855304/1‑251 (MQ=255)                     tttCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAgc                                                                                                                                                          <  1:306672/251‑2 (MQ=255)                     tttCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAgc                                                                                                                                                          <  1:114208/251‑2 (MQ=255)                      ttCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAgcg                                                                                                                                                         >  1:682349/1‑249 (MQ=255)                            cccTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGgaga                                                                                                                                                     >  1:699484/1‑249 (MQ=255)                             ccTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATAc                                                                                                                                                                    >  1:343120/1‑233 (MQ=255)                             ccTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTa                                                                                                                                                  >  1:13135/1‑251 (MQ=255)                               tGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATg                                                                                                                                                <  1:193416/251‑1 (MQ=255)                                gAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGt                                                                                                                                               >  1:375888/1‑251 (MQ=255)                                       gTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTAttt                                                                                                                                                               <  1:83111/228‑1 (MQ=255)                                        tGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatctgtt                                                                                                                                        <  1:813078/250‑8 (MQ=255)                                        tGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatc                                                                                                                                            <  1:429730/246‑4 (MQ=255)                                        tGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatc                                                                                                                                            <  1:119626/246‑4 (MQ=255)                                        tGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatc                                                                                                                                            <  1:306997/246‑4 (MQ=255)                                              ggTCGACAATGACCTAGGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatctgtttatagc                                                                                                                                  >  1:111620/2‑237 (MQ=255)                                              ggTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatctgtttatagc                                                                                                                                  >  1:544388/2‑237 (MQ=255)                                               gTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatctgtttatagcc                                                                                                                                 <  1:447034/250‑15 (MQ=255)                                                        gACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatctgtttatagc                                                                                                                                  <  1:840304/240‑14 (MQ=255)                                                         aCCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACAt                                                                                                                                                                             >  1:764722/1‑196 (MQ=255)                                                         aCCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCaaaccg                                                                                                                        >  1:144065/1‑246 (MQ=255)                                                          ccTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTAt                                                                                                                                                                 >  1:530847/1‑207 (MQ=255)                                                             aaGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAg                                                                                                                  >  1:125742/1‑251 (MQ=255)                                                                          aCTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAgc                                                                                                                                                          <  1:873273/198‑2 (MQ=255)                                                                          aCTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCaatac                                                                                                     <  1:769804/251‑4 (MQ=255)                                                                           cTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATAc                                                                                                                                                                    <  1:186589/187‑1 (MQ=255)                                                                               aaaTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGAc                                                                                                                                                    <  1:326446/199‑1 (MQ=255)                                                                               aaaTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatc                                                                                                                                            >  1:297848/1‑204 (MQ=255)                                                                               aaaTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATg                                                                                                <  1:98335/251‑1 (MQ=255)                                                                                 aTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCaatac                                                                                                     >  1:187440/1‑241 (MQ=255)                                                                                 aTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGAt                                                                                              <  1:626647/251‑1 (MQ=255)                                                                                 aTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGAt                                                                                              <  1:10668/251‑1 (MQ=255)                                                                                    gTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGt                                                                                           >  1:218883/1‑251 (MQ=255)                                                                                           gggATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATg                                                                                                <  1:209340/239‑1 (MQ=255)                                                                                            ggATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGa                                                                                      <  1:667185/248‑1 (MQ=255)                                                                                              aTTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAg                                                                                 >  1:553043/1‑251 (MQ=255)                                                                                               ttACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGta                                                                                                                                              >  1:80914/1‑188 (MQ=255)                                                                                                    aGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGtt                                                                               <  1:8770/247‑1 (MQ=255)                                                                                                     ggCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatc                                                                                                                                            <  1:550103/185‑4 (MQ=255)                                                                                                      gCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTaaagtt                                                                         >  1:283197/1‑248 (MQ=255)                                                                                                       cTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATAc                                                                                                                                                                    <  1:541637/159‑1 (MQ=255)                                                                                                       cTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGttt                                                                        >  1:263498/1‑251 (MQ=255)                                                                                                         aTTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTa                                                                      >  1:105407/1‑251 (MQ=255)                                                                                                          tttAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAg                                                                     <  1:35820/251‑1 (MQ=255)                                                                                                                  cAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggc                                                              <  1:366099/250‑4 (MQ=255)                                                                                                                        aTGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggcagtggtt                                                       >  1:605186/1‑241 (MQ=255)                                                                                                                                   ggTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTAt                                                                                                                                                                 <  1:838395/134‑1 (MQ=255)                                                                                                                                   ggTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacg                                            >  1:412097/1‑249 (MQ=255)                                                                                                                                   ggTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaag                                                >  1:141345/1‑247 (MQ=255)                                                                                                                                    gTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacgg                                           >  1:682175/1‑248 (MQ=255)                                                                                                                                       cTTTTCTATCCAACAAATCTATTATCTTTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGAt                                                                                              <  1:432130/197‑1 (MQ=255)                                                                                                                                         tttCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattt                                      >  1:564018/1‑243 (MQ=255)                                                                                                                                                 caacaaATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCTTATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggc                                                              <  1:841199/219‑4 (MQ=255)                                                                                                                                                      aaTCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggcagtggtt                                                       >  1:716179/1‑211 (MQ=255)                                                                                                                                                       aTCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCCAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattttgctcaagtagcag                        >  1:87582/1‑229 (MQ=255)                                                                                                                                                        tCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAgc                                                                                                                                                          <  1:139195/120‑2 (MQ=255)                                                                                                                                                        tCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGtatc                                                                                                                                            >  1:85459/1‑131 (MQ=255)                                                                                                                                                        tCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattttgctcaagtagcag                        >  1:154701/1‑228 (MQ=255)                                                                                                                                                          tattatCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattttgctcaagtagcagaac                     <  1:47873/251‑26 (MQ=255)                                                                                                                                                               tctctcCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggc                <  1:740715/251‑2 (MQ=255)                                                                                                                                                                  ctcCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggca               >  1:91141/1‑247 (MQ=255)                                                                                                                                                                  ctcCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCCATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggca               >  1:54289/1‑247 (MQ=255)                                                                                                                                                                           aTCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAgc                                                                                                                                                          <  1:344984/101‑2 (MQ=255)                                                                                                                                                                             ccTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggcactgaaactgagaa  >  1:167761/1‑236 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                       TGATTATATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTAGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAACGATCCTAGATTCTTTACATCCAGCCTAGATCCGAATCAAGAAGTTGCATATCGCATCTATGTTAAATCTGGTGATGATTACATTCTAAATACTATTTTAAGTGGAGACTATGTTGCTATCTATAGCCAACAAATTGCTTCAGATTCGCGCAACACTTATGATCGTAAAGATTTAGTTAAAATTTTAGGTAAAGTATTGGTTTTACAAGCAAGAAACTTTACTCAATCAACAGAACAAGGTGCTCCTGAAGCTGA  >  1/1479402‑1479824

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑