Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,574,257	G→A	intergenic (+326/+2655)	yjiT → / ← mcrC	putative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,574,257	0	G	A	100.0%	140.2 / NA	41	intergenic (+326/+2655)	yjiT/mcrC	putative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/18);  total (23/18)

GCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGATATTGAACAGGTTTGTAGTG  >  NC_000913/4574012‑4574495                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               gCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGaaa                                                                                                                                                                                                                                            >  1:16658/1‑250 (MQ=255)      tCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGATATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATaaa                                                                                                                                                                                                                                       <  1:340182/250‑1 (MQ=255)              tCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGAt                                                                                                                                                                                                                              >  1:18688/1‑251 (MQ=255)                      cATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATtctc                                                                                                                                                                                                                      <  1:592604/251‑1 (MQ=255)                        tatCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGCTTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTCTGACTGGGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGACAATAAGGGTGGCGATGAATTCTCtt                                                                                                                                                                                                                    >  2:126815/1‑251 (MQ=255)                                            aGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGAACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGtt                                                                                                                                                                                                >  2:221143/1‑251 (MQ=255)                                                aaTATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAAc                                                                                                                                                                                            >  1:184510/1‑251 (MQ=255)                                                  tatGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGAGATTGTGGAAGCACTGTATAAGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTACGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACAt                                                                                                                                                                                          >  1:453895/1‑251 (MQ=255)                                                   atGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACAt                                                                                                                                                                                          <  2:248986/250‑1 (MQ=255)                                                          tgtgTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTAc                                                                                                                                                                                  >  1:184057/1‑251 (MQ=255)                                                                       aTTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTc                                                                                                                                                                     <  1:49635/251‑1 (MQ=255)                                                                                               cGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTGCAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCt                                                                                                                                             <  1:542277/251‑1 (MQ=255)                                                                                                       aaaGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGat                                                                                                                                     >  2:203894/1‑251 (MQ=255)                                                                                                        aaGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAATCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGatt                                                                                                                                    >  2:346257/1‑251 (MQ=255)                                                                                                               gCCCTTTGGCTCTGTGTATGTGAACCAGCGACAAAATTGAGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGACTCATATTCTGGCATGGGCCAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAATTCGAAAAAGACAACAGTTTTCAGCATCCGTTGATTATTNCCg                                                                                                                             >  2:591183/1‑251 (MQ=255)                                                                                                                cccTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGc                                                                                                                            >  1:657532/1‑251 (MQ=255)                                                                                                                    ttGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAgg                                                                                                                         <  1:276839/250‑1 (MQ=255)                                                                                                                      ggCTATGTGTATGTGAACCAGCGAAAAAATTGTGACTGTGGAAGCTCTGTATACGAAGTTGCATTCTGGAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGGGAGATAATAAAGGTGGCGAAGAATTCTCTTTGCAGGATGAAGTACCTGTTGTACATGACGCTACAGAAGAAAAAGTCGAAAACGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAgg                                                                                                                         >  2:17322/1‑248 (MQ=255)                                                                                                                        cTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGgag                                                                                                                       >  1:336060/1‑248 (MQ=255)                                                                                                                        cTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGgag                                                                                                                       >  2:61192/1‑248 (MQ=255)                                                                                                                                atgtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGa                                                                                                             <  2:128932/250‑1 (MQ=255)                                                                                                                                    gAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCa                                                                                                         <  1:164989/250‑1 (MQ=255)                                                                                                                                      aCCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGt                                                                                                      <  2:437527/251‑1 (MQ=255)                                                                                                                                           cGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATAc                                                                                                 <  2:483604/251‑1 (MQ=255)                                                                                                                                             aCAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATAcc                                                                                                <  1:649026/250‑1 (MQ=255)                                                                                                                                             aCAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATAcc                                                                                                >  2:649026/1‑250 (MQ=255)                                                                                                                                                    tgtgACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAAc                                                                                         <  1:375097/250‑1 (MQ=255)                                                                                                                                                                               gTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTgg                                                             >  2:62485/1‑251 (MQ=255)                                                                                                                                                                                  gCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTTGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTCCAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCg                                                           <  2:191392/250‑1 (MQ=255)                                                                                                                                                                                     ttCTGAAATGATTGAAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGAAGCGCCGTATTGGCGttat                                                       >  2:113935/1‑251 (MQ=255)                                                                                                                                                                                               aTTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGGAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGAc                                             >  1:513437/1‑251 (MQ=255)                                                                                                                                                                                                       gAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTcc                                     >  2:281268/1‑251 (MQ=255)                                                                                                                                                                                                                     ggCACGGGACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGAttatta                          <  2:441344/248‑1 (MQ=255)                                                                                                                                                                                                                          gggACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAGAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGATAtt                  >  1:265049/1‑251 (MQ=255)                                                                                                                                                                                                                           ggACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGATATTg                 <  1:568166/251‑1 (MQ=255)                                                                                                                                                                                                                           ggACAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGAGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGATATTg                 >  1:42627/1‑251 (MQ=255)                                                                                                                                                                                                                              cAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCt                                                                                                                                             >  1:338003/1‑124 (MQ=255)                                                                                                                                                                                                                              cAAAAAGGGCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCt                                                                                                                                             <  2:338003/124‑1 (MQ=255)                                                                                                                                                                                                                                     ggCAAACTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGATATTGAACAGGTTTg       <  1:155032/251‑1 (MQ=255)                                                                                                                                                                                                                                       cAAACTAGTCAAGTGAGAAAATAAAGGGGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGCAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGGTATTGAACAGGTTTGt      >  2:89606/1‑250 (MQ=255)                                                                                                                                                                                                                                          aCTAGTCCAGTGAGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCGCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGATATTGAACAGGTTTGTAGTg  <  2:83864/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               GCATTTCTTAAAGTCAGGGCGCATATCTTCCAGCGTAATACCTAGGGAATATGGGCTTGTGTTGATAAAGATTGCAGATTAAAGCACGGTACACCGCTCGACAAAGGCTGGCCCTTTGGCTATGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCAGGTTATACCCTACAACGCCTCGACCGTCAGACGCGCCGTATTGGCGTTATTAACAATGACAGCATTCCGCTGATTATTAATGATATTGAACAGGTTTGTAGTG  >  NC_000913/4574012‑4574495

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