Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,456,773	T→C	G167G (GGT→GGC)	unknown →	ParA‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,456,773	0	T	C	100.0%	155.1 / NA	45	G167G (GGT→GGC)	unknown	ParA‑like protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/17);  total (28/17)

GCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAACGAAAACTTTGAAATACCGTTATTAAAGGGTAGTACGGCGGTACGCGTCAGCTATGCAGAAGTCTTAAGTTCCGGTAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGCACTAGGAGTTCAATATGAGTAAATTAACAAAAGGTACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACGGCTATGCCACAGCTCAAAGA  >  1/1456543‑1456992                                                                                                                                                                                                                                       |                                                                                                                                                                                                                            gCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAg                                                                                                                                                                                                         >  1:227459/1‑251 (MQ=255)    aTGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTa                                                                                                                                                                                                      <  1:480745/251‑1 (MQ=255)    aTGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTa                                                                                                                                                                                                      <  1:542422/251‑1 (MQ=255)               gattAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCt                                                                                                                                                                                                       >  1:397924/1‑239 (MQ=255)                      cTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAgctc                                                                                                                                                                                    <  1:85244/251‑4 (MQ=255)                                     cTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGtttt                                                                                                                                                                     >  1:177639/1‑251 (MQ=255)                                                       tCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   <  1:683301/251‑1 (MQ=255)                                                       tCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   <  1:590249/251‑1 (MQ=255)                                                            aTTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAAGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGa                                                                                                                                              >  1:425982/1‑251 (MQ=255)                                                                 atatTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGc                                                                                                                                          <  1:546383/250‑1 (MQ=255)                                                                 atatTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCa                                                                                                                                         >  1:404339/1‑251 (MQ=255)                                                                                 aGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:49132/1‑225 (MQ=255)                                                                                  gattcgatagcgCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTaatcaaatca                                                                                                                         >  1:740660/12‑250 (MQ=255)                                                                                  gATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAgc                                                                                                                                                                                      >  1:33978/1‑188 (MQ=255)                                                                                      cgatagAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   <  1:242778/215‑1 (MQ=255)                                                                                           gAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAAt                                                                                                               >  1:237924/1‑251 (MQ=255)                                                                                                ggTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCg                                                                                                                                               <  1:662222/214‑1 (MQ=255)                                                                                                  tGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAAc                                                                                                         <  1:8772/250‑1 (MQ=255)                                                                                                       gCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAgc                                                                                                                                                                                      <  1:620220/168‑2 (MQ=255)                                                                                                        cACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTa                                                                                                  <  1:525285/251‑1 (MQ=255)                                                                                                          cGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAgg                                                                                                >  1:200031/1‑251 (MQ=255)                                                                                                            tATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGa                                                                                               <  1:96834/250‑1 (MQ=255)                                                                                                                 aaaTTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCaa                                                                                         >  1:47557/1‑251 (MQ=255)                                                                                                                    ttAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAAGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAatat                                                                                      <  1:299934/251‑1 (MQ=255)                                                                                                                    ttAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTAGTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAatat                                                                                      <  1:430768/251‑1 (MQ=255)                                                                                                                              gCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGa                                                                                               >  1:530502/1‑232 (MQ=255)                                                                                                                                        ggtttAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:351420/4‑170 (MQ=255)                                                                                                                                        ggtttAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAggca                                                                  >  1:50307/4‑249 (MQ=255)                                                                                                                                        ggtttAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAggca                                                                  >  1:703895/4‑249 (MQ=255)                                                                                                                                              gTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTaa                                                            >  1:674597/1‑251 (MQ=255)                                                                                                                                                      tcaaacaaataaacaaaccaaTTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATcc                                                    <  1:335871/232‑1 (MQ=255)                                                                                                                                                           caaataaacaaaccaaTTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAgg                                                                                                >  1:471212/15‑202 (MQ=255)                                                                                                                                                           caaataaacaaaccaaTTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCggg                                                 >  1:334441/15‑249 (MQ=255)                                                                                                                                                                                    agatgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:381207/13‑126 (MQ=255)                                                                                                                                                                                      atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAg                                                                                              >  1:587947/11‑177 (MQ=255)                                                                                                                                                                                      atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAcag                     >  1:508724/11‑248 (MQ=255)                                                                                                                                                                                      atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAcag                     >  1:666918/11‑248 (MQ=255)                                                                                                                                                                                      atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGGATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAcagc                    <  1:438021/241‑4 (MQ=255)                                                                                                                                                                                        gttattgagTTTATTAAAGATAACTTTGAAATCCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACAGCTa                  >  1:551538/9‑251 (MQ=255)                                                                                                                                                                                                    attaaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTaaa                                                                               >  1:517246/10‑178 (MQ=255)                                                                                                                                                                                                    attaaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACAGCTATGCCACAGCTCa      >  1:585959/10‑251 (MQ=255)                                                                                                                                                                                                       aaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:58835/7‑107 (MQ=255)                                                                                                                                                                                                       aaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATggg                                      >  1:653384/7‑216 (MQ=255)                                                                                                                                                                                                        aagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACAGCTATGCCACAGCTCAAaga  <  1:334682/246‑1 (MQ=255)                                                                                                                                                                                                               cTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAc                       >  1:835765/1‑223 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                            GCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAACGAAAACTTTGAAATACCGTTATTAAAGGGTAGTACGGCGGTACGCGTCAGCTATGCAGAAGTCTTAAGTTCCGGTAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGCACTAGGAGTTCAATATGAGTAAATTAACAAAAGGTACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACGGCTATGCCACAGCTCAAAGA  >  1/1456543‑1456992

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