Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,821,188	A→G	intergenic (‑18/‑240)	ligB ← / → gmk	DNA ligase B/guanylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,821,188	0	A	G	100.0%	173.0 / NA	49	intergenic (‑18/‑240)	ligB/gmk	DNA ligase B/guanylate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/31);  total (18/31)

CAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACCTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGATAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACT  >  NC_000913/3820945‑3821426                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               cAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAg                                                                                                                                                                                                                                         <  1:136901/251‑1 (MQ=255)    agcagcGGTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCa                                                                                                                                                                                                                                      <  2:18798/251‑1 (MQ=255)                           cGGGCGCTTAACTGATCGTAAACACCGTCTTCCGCCTCACTTTTTCCTTCCTGCCAGTAGTCATCGTCCCACTGTTTTAGTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGatagcaa                                                                                                                                                                                                               <  2:58618/251‑1 (MQ=255)                            gggCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGGAATTTTTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCCCCCACACCGATGATTGCCCGCACAAGATACTTATTAATATCGCCATCCAGACTTTCATCCTCACCTCCGCTGTAGGGAGTCAGGCAGATATACAACGTGATGaaagcaa                                                                                                                                                                                                               >  1:38983/1‑250 (MQ=255)                              gcgcTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGaaagcaaagc                                                                                                                                                                                                            <  1:317348/251‑1 (MQ=255)                                              aaaCACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGTCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCATGCAGATATACAACGTCATTAAAGCAAAGCCGTGTGGCAACAACg                                                                                                                                                                                             <  1:171418/250‑1 (MQ=255)                                                             ccTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTTGGGACCAGGCCGGACAGACCGACCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCCTCCATACTTTCATCCTAACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGACAAGCCGCGTGGCAAAAACGGAGTCTGCGAGGAcac                                                                                                                                                                             >  2:249200/1‑249 (MQ=255)                                                               tCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTTTTCCTTTGCTCTGGCTGGCGACCAGGCCGGACAGACCCCCCACACAGATGATTGCCAGCCCAAGAAACTTATTAATATCGCCCTCCCTACTTTCATCATCACCTCCGCTGTTGGTAGTTAGGCAGATTTACAACGTGATGAAAGCAAAACCGAGGGGCAAAAAAGGAGTCTGCGAGGACGCtt                                                                                                                                                                           >  2:215623/1‑251 (MQ=255)                                                               tCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTTTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCtt                                                                                                                                                                           >  1:106579/1‑251 (MQ=255)                                                                               tccAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCt                                                                                                                                                             <  1:299125/249‑1 (MQ=255)                                                                                ccattagtCATCGTCCCACTGTTTTATTGGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCTGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCATATATGCAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCATCCTGAAAATCGTCttt                                                                                                                                                           <  2:282210/246‑1 (MQ=255)                                                                                           cGTCCCACTGGTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCAGCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAAt                                                                                                                                               <  2:225445/251‑1 (MQ=255)                                                                                            gtccctgtGTTTTATTTTCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCTGGACAAACCGCCCACGCGGATGATTGCCAGCACAAGATACTGATTAATATCGCCATCCATACTTTCATCATCACCTCTGCTGTAGGTAGTCAGGCAGATATACAATGTGATGAAAGCAAAGCCGATTGGCAAATACGGAGTCTGCGAGGACGCTTGCTGAAAATCTTCTTTAttacattgattg                                                                                                                                              <  2:199982/244‑11 (MQ=255)                                                                                              cccACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGAc                                                                                                                                            <  2:393082/251‑1 (MQ=255)                                                                                              cccACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAAGCGTCTTTATTGCAGTGAATGAc                                                                                                                                            <  2:303822/251‑1 (MQ=255)                                                                                              cccACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCAAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGa                                                                                                                                             <  1:285007/250‑1 (MQ=255)                                                                                               ccACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACa                                                                                                                                           <  1:378363/251‑1 (MQ=255)                                                                                                 aCTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGGCGGGCAGAACGCCCACACAGATGATTGCCAACACAAGATACTTATTAATATCGCCATCCCTACTTTCATCATCACCTCCGCTGTAGGTAGGCCGGCAGATATACAACGTGATGAAAGCAAAGCCGCGTGGCAAAAACGGAGTCTGCGGGGGCGGCTCCTGACAATCGTCTTTATTGCAGGGAATGAAAgg                                                                                                                                         >  2:388750/1‑251 (MQ=255)                                                                                                    gTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGGCAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCCTACTTTTATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAACCGAg                                                                                                                                                                                                        >  1:318500/1‑185 (MQ=255)                                                                                                    gTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAg                                                                                                                                                                                                        <  2:318500/185‑1 (MQ=255)                                                                                                     ttttATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAAc                                                                                                                                                                                              <  2:103012/194‑1 (MQ=255)                                                                                                     ttttATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAAATGATTGCCCGCACAAGATACTTATTAATATCGCCATCCATACTTTCCTCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAAc                                                                                                                                                                                              >  1:103012/1‑194 (MQ=255)                                                                                                             gCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTGTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAag                                                                                                                             <  1:227431/251‑1 (MQ=255)                                                                                                             gCTGTTGCAGGCGGGAAAATTTTTCCTGTGCTCTGGCTGGCGACCAGGCCGGGCAGACCGCCCACACAGATGATTGCCAGCACAAAATACTTATTAATATCGCCCTCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAAGCAGATATACAACGTGCTGAAAGCAAAGCCCAGTGGCAAAAAAAGAGTCTGCGAGGGCGCGTCCTGAAAATCGTCTTTTTTTCAGTGAATGACAGGCAAATGCGGa                                                                                                                               >  2:139019/1‑249 (MQ=255)                                                                                                                          ggAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCaaa                                                                                                                 >  1:216717/1‑250 (MQ=255)                                                                                                                                    tcCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGcaaca                                                                                                         <  1:121669/248‑1 (MQ=255)                                                                                                                                      cTGTGCTCTGGATGGCGACCAGGACGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCCGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACttt                                                                                                    >  2:318794/1‑251 (MQ=255)                                                                                                                                                       aCCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGc                                                                                   <  2:106579/251‑1 (MQ=255)                                                                                                                                                       aCCAGGCCGGAAAGACCGACCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCCTCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCACAAAGTGTGAGc                                                                                   >  2:66299/1‑251 (MQ=255)                                                                                                                                                                    gACCGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTc                                                                       >  1:207924/1‑250 (MQ=255)                                                                                                                                                                       cGCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACAt                                                                   <  2:207924/251‑1 (MQ=255)                                                                                                                                                                        gCCCACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATGGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGGTACGTCACATg                                                                  <  2:40310/251‑1 (MQ=255)                                                                                                                                                                         cccACACAGATGATTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATgg                                                                 <  1:136708/251‑1 (MQ=255)                                                                                                                                                                                    gatTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCgtgt                                                      >  1:432318/1‑251 (MQ=255)                                                                                                                                                                                    gatTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCgtgt                                                      <  2:299039/251‑1 (MQ=255)                                                                                                                                                                                    gatTGCCAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACGTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCgtg                                                       <  2:368544/250‑1 (MQ=255)                                                                                                                                                                                          cAGCACAAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTataat                                                 >  1:102542/1‑250 (MQ=255)                                                                                                                                                                                               caAGATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGACGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTc                                           >  1:199468/1‑251 (MQ=255)                                                                                                                                                                                                  gATACTTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTa                                        <  2:216717/251‑1 (MQ=255)                                                                                                                                                                                                     acttATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATg                                      <  2:102542/247‑1 (MQ=255)                                                                                                                                                                                                     acttATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGt                                     <  1:147889/248‑1 (MQ=255)                                                                                                                                                                                                               atCGCCATCCATACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTAtt                           >  2:121498/1‑251 (MQ=255)                                                                                                                                                                                                                        catACTTTCATCATCACCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATc                  <  1:340809/251‑1 (MQ=255)                                                                                                                                                                                                                         atACTTTCATCATCACCTCCGCTGTCGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATc                  >  2:155955/1‑250 (MQ=255)                                                                                                                                                                                                                               tcatcatcaCCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATAc            <  1:273322/250‑1 (MQ=255)                                                                                                                                                                                                                                      caCCTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACATGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATa    <  2:416803/251‑1 (MQ=255)                                                                                                                                                                                                                                        ccTCCGCTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACAGACGAAAGATACt  <  1:326489/251‑1 (MQ=255)                                                                                                                                                                                                                                              cTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAACATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAAt                   >  1:187002/1‑228 (MQ=255)                                                                                                                                                                                                                                              cTGTAGGTAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCGTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGGGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAAt                   <  2:187002/228‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               CAAAGCAGCGTTGCCACTGCGTAAGACGGGCGCTTAACTGATCGTAAACACCGTCTTCCACCTCACTTTTTCCTTCCTTCCAGTAGTCATCGTCCCACTGTTTTATTTGCTGTTGCAGGCGGGAAATTTCTTCCTGTGCTCTGGCTGGCGACCAGGCCGGACAGACCGCCCACACAGATGATTGCCAGCACAAGATACCTATTAATATCGCCATCCATACTTTCATCATCACCTCCGCTGTAGATAGTCAGGCAGATATACAACGTGATGAAAGCAAAGCCGAGTGGCAAAAACGGAGTCTGCGAGGACGCTTCCTGAAAATCGTCTTTATTGCAGTGAATGACAGGCAAATGCGGAAGCAGCTACGCAAAACGCAACAACTTTGCGCAAAAAGTGTGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACT  >  NC_000913/3820945‑3821426

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