Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,135,541	T→A	G52G (GGA→GGT)	unknown ←	Type IV secretion system protein VirD4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,135,541	0	T	A	100.0%	127.2 / NA	48	G52G (GGA→GGT)	unknown	Type IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (26/22);  total (26/22)

TTCCCGCGCTTCCGAGGCAGTCGCAAATTCGGCGTTACCATAAAGATTCTTTTTAGGCTGATAGAAGAGTTTCAGCATAATCAGTAGGCAGACGACAAAACTAATCAATAAG‑GATTGCCAGATTTTCGGATAATTTCCCAGGATGGCATCCAGCGTAAGAGTAGGGTAGTAACGCTTGAACGCTTCAATCAGACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGAAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCGTGTTCATGACTGTTGGTA  >  1/1135406‑1135780                                                                                                                                         |                                                                                                                                                                                                                                                tttttgtacttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAg                                                                                                                              <  1:878219/229‑1 (MQ=255)    ttgtacttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGa                                                                                                                           >  1:339015/20‑251 (MQ=255)      gtacttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACa                                                                                                                         >  1:377257/18‑251 (MQ=255)         cttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACaaa                                                                                                                       >  1:175940/15‑250 (MQ=255)         cttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAAc                                                                                                                      <  1:102095/237‑1 (MQ=255)         cttctgattttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCCCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAAc                                                                                                                      <  1:438044/237‑1 (MQ=255)                 tttgtggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGaataa                                                                                                              <  1:820065/245‑1 (MQ=255)                     tggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGc                                                                                                           <  1:492742/248‑1 (MQ=255)                     tggCAAATTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGc                                                                                                           <  1:724521/248‑1 (MQ=255)                         aaaTTCAGCGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAatat                                                                                                      <  1:106068/251‑1 (MQ=255)                                 cGTTACCATAAAGGGATTTTTTAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTACGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCtttt                                                                                              >  1:275186/1‑251 (MQ=255)                                     accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGaataa                                                                                                              <  1:322697/218‑1 (MQ=255)                                     accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTtca                                                                                          >  1:665924/14‑251 (MQ=255)                                     accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTtca                                                                                          >  1:342106/14‑251 (MQ=255)                                     accataaagggattttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTtca                                                                                          >  1:504395/14‑251 (MQ=255)                                                 ttttttATGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGaataa                                                                                                              >  1:839431/2‑219 (MQ=255)                                                 ttttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTcatcat                                                                                      >  1:785354/2‑243 (MQ=255)                                                  tttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGaaaa                                                                                                                                              <  1:160864/186‑1 (MQ=255)                                                  tttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACg                                                                                                                                          >  1:802996/1‑190 (MQ=255)                                                   ttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCGCTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAacctacct                                                                            <  1:679974/251‑1 (MQ=255)                                                   ttttAGGTTGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAacctacct                                                                            <  1:248688/251‑1 (MQ=255)                                                        ggttGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGc                                                                                                           >  1:147112/4‑215 (MQ=255)                                                           tGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACtt                                                                                                                                                    <  1:294301/171‑1 (MQ=255)                                                           tGATAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGaa                                                                    >  1:476304/1‑251 (MQ=255)                                                             aTAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAATTGGCATCTAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTaa                                                                                                                                                              >  1:388841/1‑159 (MQ=255)                                                             aTAGAATAGTTTCAGCATGATTAGTAAGCTCACTACAAAACTAAT‑AAGGAGAGACTGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAg                                                                                                                              >  1:797784/1‑191 (MQ=255)                                                                                              ctacaaaactaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCg                                  >  1:392351/24‑251 (MQ=255)                                                                                              ctacaaaactaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCg                                  >  1:486176/24‑251 (MQ=255)                                                                                                   aaactaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCAtg                             <  1:786281/233‑1 (MQ=255)                                                                                                   aaactaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATCAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCAtg                             <  1:463098/233‑1 (MQ=255)                                                                                                     actaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAg                                         <  1:314391/221‑1 (MQ=255)                                                                                                     actaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCCTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATACAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCAtgtg                           >  1:443342/17‑251 (MQ=255)                                                                                                      ctaatatggatagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATgtgt                          >  1:869315/16‑251 (MQ=255)                                                                                                      ctaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAatat                                                                                                      <  1:229109/160‑1 (MQ=255)                                                                                                      ctaataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTcatcat                                                                                      >  1:326821/16‑191 (MQ=255)                                                                                                        aataaggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAac                                                                                  >  1:317060/14‑193 (MQ=255)                                                                                                            aggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCg                    >  1:543871/10‑251 (MQ=255)                                                                                                             ggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCt                                                                                                 >  1:120654/9‑173 (MQ=255)                                                                                                             ggagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGAtt                                      >  1:276380/9‑232 (MQ=255)                                                                                                              gagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCg                    <  1:344569/242‑1 (MQ=255)                                                                                                              gagagactGCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCg                    <  1:318898/242‑1 (MQ=255)                                                                                                                      gCCATATTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAacctacct                                                                            >  1:93268/1‑185 (MQ=255)                                                                                                                       ccataTTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCGTGTTCATGAc          <  1:642513/246‑1 (MQ=255)                                                                                                                         ataTTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCGTGTTCATGACTGt       <  1:361354/249‑1 (MQ=255)                                                                                                                           aTTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTg                                                                                                    >  1:220848/1‑156 (MQ=255)                                                                                                                           aTTTTCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCGTGTTCATGAc          <  1:782245/246‑1 (MQ=255)                                                                                                                            ttttCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCGTGTTCATGACTGTTgg    <  1:320157/251‑1 (MQ=255)                                                                                                                              ttCGGATAATTACCCCAAATGGCATCCAGCGTTAGAGTAGGGTAGTAGCGCTTAAATGCTTCGATTAAACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGGAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCGTGTTCATGACTGTTGGTa  >  1:581010/1‑251 (MQ=255)                                                                                                                                         |                                                                                                                                                                                                                                                TTCCCGCGCTTCCGAGGCAGTCGCAAATTCGGCGTTACCATAAAGATTCTTTTTAGGCTGATAGAAGAGTTTCAGCATAATCAGTAGGCAGACGACAAAACTAATCAATAAG‑GATTGCCAGATTTTCGGATAATTTCCCAGGATGGCATCCAGCGTAAGAGTAGGGTAGTAACGCTTGAACGCTTCAATCAGACGAATTTTATTTGCCAGTAGAAGTAAAAGCCCACTTAGAAAATACGTCGGAAGCAAAGTGACAAACAGAAATAATGCAATATTGGCTTTTTTCATCATTAACCTACCTTGAGCGAATTCAGTTTAGCTTGCGGATCAAAATAGATTTCCGTCATGTGTCTGCCGTGTTCATGACTGTTGGTA  >  1/1135406‑1135780

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