Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	2	27,305	T→C	L20L (TTG→CTG)	unknown →	Helicase PriA essential for oriC/DnaA‑independent DNA replication

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	2	27,305	0	T	C	100.0%	208.2 / NA	75	L20L (TTG→CTG)	unknown	Helicase PriA essential for oriC/DnaA‑independent DNA replication
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/41);  total (34/41)

TTCATTCATCCTGTACGCATTTCTTGCTTCTTAAATTTGACAGTATTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCACTATTTTGATTATATCTTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCGCAGGCGTTGC  >  2/27083‑27545                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  ttcattcaTCCTGTACGCATTTCTTGCTTCTTAAATTTGACAGTATTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑ggg                                                                                                                                                                                                                      <  1:22450/251‑1 (MQ=255)          ccTGTACGCATTTCTTGCTTCTTAAATTTGACAGTATTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCTCTGAATGCTGTCGT‑GGGGGCGAGGgt                                                                                                                                                                                                             >  1:804732/1‑251 (MQ=255)             gTACGCATTTCTTGCTTCTTAAATTTGACAGTATTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGtgtg                                                                                                                                                                                                          >  1:416800/1‑251 (MQ=255)                             tcttAAATTTGACAGTATTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGgcg                                                                                                                                                                                          >  1:382399/1‑251 (MQ=255)                                             aTTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGc                                                                                                                                                                                                                                <  1:396367/198‑1 (MQ=255)                                             aTTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCa                                                                                                                                                                          >  1:489436/1‑251 (MQ=255)                                                 tATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCAtggt                                                                                                                                                                      <  1:169291/251‑1 (MQ=255)                                                    cccTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAggg                                                                                                                                                                                                              <  1:4908/208‑1 (MQ=255)                                                    cccTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCAtggtgg                                                                                                                                                                    >  1:841895/1‑250 (MQ=255)                                                    cccTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATggtggt                                                                                                                                                                   >  1:106555/1‑251 (MQ=255)                                                        tATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAg                                                                                                                                                               <  1:786387/251‑1 (MQ=255)                                                         aTCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGCCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAgg                                                                                                                                                              >  1:44698/1‑251 (MQ=255)                                                                    aTTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAggg                                                                                                                                                                                                              <  1:698429/192‑1 (MQ=255)                                                                         aCCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAg                                                                                                                                                                                    <  1:23511/213‑1 (MQ=255)                                                                           cccAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGAtg                                                                                                                                            <  1:723524/251‑1 (MQ=255)                                                                                   acaGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGC‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGGGGCGGAGCAGTCGGa                                                                                                                                              >  1:381330/1‑241 (MQ=255)                                                                                     ggagCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGt                                                                                                                                                             <  1:303116/223‑1 (MQ=255)                                                                                      gagCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAgg                                                                                                                                                              <  1:329096/222‑1 (MQ=255)                                                                                      gagCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAg                                                                                                                                 >  1:614531/1‑251 (MQ=255)                                                                                      gagCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAg                                                                                                                                 >  1:510697/1‑251 (MQ=255)                                                                                           aTTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGAtgt                                                                                                                                           >  1:889957/1‑236 (MQ=255)                                                                                            ttCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTg                                                                                                                           <  1:667678/251‑1 (MQ=255)                                                                                               tcCCCCAAATAGACAACAAGCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAg                                                                                                                        <  1:494790/251‑1 (MQ=255)                                                                                                cccccAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑gg                                                                                                                                                                                                                       <  1:330281/155‑1 (MQ=255)                                                                                                 ccccAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGAc                                                                                                                      <  1:879609/251‑1 (MQ=255)                                                                                                 ccccAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGAc                                                                                                                      <  1:393547/251‑1 (MQ=255)                                                                                                 ccccAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGAc                                                                                                                      <  1:131291/251‑1 (MQ=255)                                                                                                 ccccAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGAc                                                                                                                      <  1:222231/251‑1 (MQ=255)                                                                                                  cccAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGt                                                                                                                                                                                                                             <  1:392846/148‑1 (MQ=255)                                                                                                  cccAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTACCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACa                                                                                                                     >  1:437130/1‑251 (MQ=255)                                                                                                    cAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGt                                                                                                                                                                                                  <  1:375051/172‑1 (MQ=255)                                                                                                    cAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTgggg                                                                                                                                                                                            <  1:303591/178‑1 (MQ=255)                                                                                                      aaTAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCa                                                                                                                                                                          >  1:530946/1‑194 (MQ=255)                                                                                                         aGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGAtgtg                                                                                                                                          >  1:806557/1‑223 (MQ=255)                                                                                                           acaacaATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCt                                                                                                            <  1:195368/251‑1 (MQ=255)                                                                                                           acaacaATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCt                                                                                                            <  1:24239/251‑1 (MQ=255)                                                                                                           acaacaATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCt                                                                                                            <  1:473478/251‑1 (MQ=255)                                                                                                                           aTGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGa                                                                                                                                              <  1:865053/201‑1 (MQ=255)                                                                                                                             gCTATATTCCGAGCCCATTTTATTTCACTGAATTTCTGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGCGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAg                                                                                                                                                               >  1:537142/1‑182 (MQ=255)                                                                                                                             gCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGc                                                                                                                                                                           <  1:405886/170‑1 (MQ=255)                                                                                                                             gCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGt                                                                                                                                                             >  1:503640/1‑184 (MQ=255)                                                                                                                             gCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGTGGTTTc                                                                                           >  1:196446/1‑250 (MQ=255)                                                                                                                             gCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTc                                                                                           >  1:384687/1‑250 (MQ=255)                                                                                                                                  atTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCt                                                                                                            <  1:151769/228‑1 (MQ=255)                                                                                                                                          cccATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAgg                                                                                               >  1:283435/1‑233 (MQ=255)                                                                                                                                           ccATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCa                                                                            >  1:250455/1‑251 (MQ=255)                                                                                                                                                ttATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGTGGGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCtt                                                                                        <  1:285116/235‑1 (MQ=255)                                                                                                                                                        cTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGAc                                                                                                                      <  1:726495/196‑1 (MQ=255)                                                                                                                                                                   ttCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCt                                                                                                            <  1:565964/195‑1 (MQ=255)                                                                                                                                                                   ttCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTggg                                                     >  1:265747/1‑250 (MQ=255)                                                                                                                                                                                cccGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGtattat                                       <  1:311922/251‑1 (MQ=255)                                                                                                                                                                                  cGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCa                                     <  1:892435/251‑1 (MQ=255)                                                                                                                                                                                  cGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCa                                     <  1:780210/251‑1 (MQ=255)                                                                                                                                                                                  cGAATTGCGCTTCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCa                                     <  1:245259/251‑1 (MQ=255)                                                                                                                                                                                            ttCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTg                           >  1:471265/1‑251 (MQ=255)                                                                                                                                                                                             tCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATtgat                                                                                                   >  1:361616/1‑178 (MQ=255)                                                                                                                                                                                             tCCTGTGCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTgg                          >  1:599089/1‑251 (MQ=255)                                                                                                                                                                                                tgtgCCGTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTggggg                       <  1:867413/251‑1 (MQ=255)                                                                                                                                                                                                      gTTACATCGCTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAgtgt                 >  1:531304/1‑251 (MQ=255)                                                                                                                                                                                                           atcgcTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCg            <  1:59698/247‑1 (MQ=255)                                                                                                                                                                                                            tcgcTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGtattat                                       <  1:392998/220‑1 (MQ=255)                                                                                                                                                                                                            tcgcTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGta                                           <  1:609835/216‑1 (MQ=255)                                                                                                                                                                                                            tcgcTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCg            >  1:452921/4‑250 (MQ=255)                                                                                                                                                                                                            tcgcTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCg            >  1:427635/4‑250 (MQ=255)                                                                                                                                                                                                              gcTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTa                                                                                       >  1:304763/2‑173 (MQ=255)                                                                                                                                                                                                               cTATTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCGCAgg        >  1:882583/1‑251 (MQ=255)                                                                                                                                                                                                                 aTTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGc                                                           <  1:516400/198‑1 (MQ=255)                                                                                                                                                                                                                 aTTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGt                                            <  1:587434/213‑1 (MQ=255)                                                                                                                                                                                                                 aTTTTGATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCg            >  1:402541/1‑245 (MQ=255)                                                                                                                                                                                                                      gATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGttt                                                                                                                                                                                                <  1:842848/60‑1 (MQ=255)                                                                                                                                                                                                                      gATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGa                      >  1:152683/1‑230 (MQ=255)                                                                                                                                                                                                                      gATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCGCAGGCGTTGc  >  1:73066/1‑250 (MQ=255)                                                                                                                                                                                                                      gATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCGCAGGCGTTGc  >  1:146798/1‑250 (MQ=255)                                                                                                                                                                                                                      gATTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGGGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGAtgtg                                                                                                                                          >  1:735429/1‑114 (MQ=255)                                                                                                                                                                                                                       aTTATATCCTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTggggg                       <  1:489376/228‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  TTCATTCATCCTGTACGCATTTCTTGCTTCTTAAATTTGACAGTATTTTATCCCTTATCCAAGTGAAATTAAACCCCAAACAACAGAGCTATTCTCCCCCAAATAGACAACAATCCCCCAAAATGCTATATTCCGAGCCCATTTTATTTCACTGAATTTCCGTTCATGTCATTTGCCCGAATTGCGCTTCCTGTGCCGTTACATCACTATTTTGATTATATCTTGCCGGCGTCGCTGAATGCTGTCGT‑GGGGGCGAGGGTGTGGGTGCCGTTTGGGGCGCAGAAGCGTGTTGGCATGGTGGTTGAGGTGGCGGAGCAGTCGGATGTGCCTGCGGAGCAGTTGAAGACAGTGCTGGCTTGCATTGATGAGGTTTCCTTATTTTCTCCGCAAAGCTGGCAGTTTTTGCAGTGGGCGGCGGGGTATTATCACGCGCCCTTGGGGGAAGTGTTGTCGCAGGCGTTGC  >  2/27083‑27545

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