Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,456,707	G→A	K145Q (AAG→CAA)	unknown →	ParA‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,456,707	0	G	A	100.0%	166.9 / NA	50	K145Q (AAG→CAA)	unknown	ParA‑like protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/26);  total (24/26)

TGCGCTTGGTCGCCAAGCGGATTTTTTAATCATTGATGGTGCGCCACGTATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAACGAAAACTTTGAAATACCGTTATTAAAGGGTAGTACGGCGGTACGCGTCAGCTATGCAGAAGTCTTAAGTTCCGGTAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGCACTAGGAGTTCAATATGAGTAAATTAACAAAAGGTACACTAA  >  1/1456485‑1456934                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    ggtgcTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttgag                                                                                                                                                                                                          <  1:661241/247‑4 (MQ=255)  gtgcTTGGTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttgagtt                                                                                                                                                                                                        >  1:508998/3‑246 (MQ=255)         gTCGCCAGGCGGATTTTTTAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTAttaa                                                                                                                                                                                                  <  1:461247/250‑4 (MQ=255)                      ttttttAATTATTGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttgag                                                                                                                                                                                                          >  1:592423/1‑226 (MQ=255)                                 ttGATGGAGCTCCCCGCATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttgag                                                                                                                                                                                                          <  1:780512/218‑4 (MQ=255)                                     tggagctccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaag                                                                                                                                                                     <  1:454616/238‑5 (MQ=255)                                       gagctccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACtt                                                                                                                                                                                         >  1:487132/12‑229 (MQ=255)                                        agctccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggc                                                                                                                                                                   <  1:756345/240‑7 (MQ=255)                                           tccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGaaa                                                                                                                                                                                    <  1:79786/223‑1 (MQ=255)                                            ccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggcagtac                                                                                                                                                              <  1:806558/245‑12 (MQ=255)                                            ccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggcagtac                                                                                                                                                              <  1:423631/245‑12 (MQ=255)                                                gcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttg                                                                                                                                                                                                            >  1:1276/3‑200 (MQ=255)                                                           gCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAg                                                                                                                                               >  1:227459/1‑251 (MQ=255)                                                              aTGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTa                                                                                                                                            <  1:542422/251‑1 (MQ=255)                                                              aTGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTa                                                                                                                                            <  1:480745/251‑1 (MQ=255)                                                               tGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggc                                                                                                                                                                   >  1:675372/1‑221 (MQ=255)                                                                         gattAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCt                                                                                                                                             >  1:397924/1‑239 (MQ=255)                                                                              aaCTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggcagtacaacgg                                                                                                                                                         <  1:735347/222‑17 (MQ=255)                                                                                cTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAgctc                                                                                                                          <  1:85244/251‑4 (MQ=255)                                                                                       aTTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAt                                                                                                                                                                                                              >  1:127949/1‑160 (MQ=255)                                                                                            aTTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTAttaa                                                                                                                                                                                                  <  1:573837/167‑4 (MQ=255)                                                                                               cTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGtttt                                                                                                           >  1:177639/1‑251 (MQ=255)                                                                                                            aaCCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggc                                                                                                                                                                   >  1:291841/1‑176 (MQ=255)                                                                                                                 tCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                         <  1:683301/251‑1 (MQ=255)                                                                                                                 tCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                         <  1:590249/251‑1 (MQ=255)                                                                                                                      aTTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaagc                                                                                                                                                                   >  1:437352/1‑166 (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)                                                                                                                                                          ggTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttg                                                                                                                                                                                                            <  1:311085/95‑2 (MQ=255)                                                                                                                                                          ggTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCg                                                                                     <  1:662222/214‑1 (MQ=255)                                                                                                                                                           gTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTAttg                                                                                                                                                                                                            <  1:561121/94‑2 (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)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    TGCGCTTGGTCGCCAAGCGGATTTTTTAATCATTGATGGTGCGCCACGTATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAACGAAAACTTTGAAATACCGTTATTAAAGGGTAGTACGGCGGTACGCGTCAGCTATGCAGAAGTCTTAAGTTCCGGTAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGCACTAGGAGTTCAATATGAGTAAATTAACAAAAGGTACACTAA  >  1/1456485‑1456934

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