Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	155,463	A→T	*247K (TAA→AAA)	yadV ←	putative fimbrial chaperone YadV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	155,463	0	A	T	100.0%	137.4 / NA	40	*247K (TAA→AAA)	yadV	putative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (22/18);  total (22/18)

GACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAAC  >  NC_000913/155229‑155699                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             gACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCTTGATTTCAGCCTGGCATTACCt                                                                                                                                                                                                                              >  2:482378/1‑251 (MQ=255) gACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCt                                                                                                                                                                                                                              <  1:167575/251‑1 (MQ=255)       gtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                                                                                                                                                                                                        >  1:117454/1‑251 (MQ=255)        taactaccagttaaaGCGGGGTTAAATTTACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCCTATTCAACAAGTTTAGCACTGAAAGCAGTATTGCACTACAGCAGAGCGCAAAACGTGGAGATGCGGCTGAGATGATGATAATTTTTAGTATATTCTATAG‑AACGGCTATGCTTCCTGTTGAATTTGATCTTTNCCTGATTTCAGCCTGGACTTACCTTCatttc                                                                                                                                                                                                                       >  2:248562/15‑247 (MQ=255)           acacCCGGGATAGCGTGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAGAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATGCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcgc                                                                                                                                                                                                                    <  1:265842/251‑1 (MQ=255)             acCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGcc                                                                                                                                                                                                                   >  2:429386/1‑250 (MQ=255)                           ggTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                                                                                                                                                                                                    <  1:10886/251‑1 (MQ=255)                              tACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                                                                                                                                                                                                  <  1:218984/250‑1 (MQ=255)                                 cTTCACTGTAACGGCTGAGATCAATATTAGATGCATTCTGCCCCATCCGGAAGGAATGGTCATATTCAACGAGTTCAGAACTGAAAGCAGTATTGCAATACATCAGCGAGGAAAACGTGGCGATAAGGGTCTGATGGTGATAATTTTTAGTATATTCTATAGTGACG‑GTGTGCTTCCTGCGGTATTATATCCTTCCCTGATTTCAGCTTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGcct                                                                                                                                                                                               >  1:306669/1‑248 (MQ=255)                                          aCCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCAc                                                                                                                                                                                     >  1:363158/1‑251 (MQ=255)                                                                     tcttccccATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTTAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCTCCAAAGTCATTAATGGCGTAAAAATGCTCTTTTGCAGAATTCTCTTTGCCATTaa                                                                                                                                                           <  2:293292/246‑1 (MQ=255)                                                                     tCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTa                                                                                                                                                            >  1:261706/1‑249 (MQ=255)                                                                              tCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTaa                                                                                                                                                           <  2:336895/241‑1 (MQ=255)                                                                              tCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTaa                                                                                                                                                           >  1:336895/1‑241 (MQ=255)                                                                                aGGAAGGTATGTGCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGAGACGGGTCAGATGATTATAATTTTTAGTTTATTCTATAGTCACG‑TTATGCTTCCTGCTGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTGGCAGAATTCGCTTTGCCATTAAGGCCATTGACtt                                                                                                                                               <  2:344961/251‑1 (MQ=255)                                                                                             cATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCtcatc                                                                                                                                  >  1:449484/1‑251 (MQ=255)                                                                                                atTCAACGAGTGCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCAc                                                                                                                                <  2:197661/250‑1 (MQ=255)                                                                                                   cAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTaaa                                                                                                                            <  1:50864/251‑1 (MQ=255)                                                                                                            cAGCACTGCAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTAAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCATTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAGATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAAt                                                                                                                   >  2:275786/1‑251 (MQ=255)                                                                                                                aCTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCAtt                                                                                                               <  1:429386/251‑1 (MQ=255)                                                                                                                     aaGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTcac                                                                                                          <  2:289096/251‑1 (MQ=255)                                                                                                                                     cagcagCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGt                                                                                          >  1:459120/1‑251 (MQ=255)                                                                                                                                            gcgcAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATNATTTTCACATCAATCGGATAGCGTTTACCGc                                                                                   >  2:228187/1‑251 (MQ=255)                                                                                                                                              gcAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGGGCAATCATTTTCTCATCAATCGGATAGCGTTTGCCGCTa                                                                                 <  1:420020/251‑1 (MQ=255)                                                                                                                                                        gCGATACGGGTCAGATTATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTTATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAAt                                                                       <  2:256261/251‑1 (MQ=255)                                                                                                                                                             aCGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATcacc                                                                   <  2:108787/250‑1 (MQ=255)                                                                                                                                                                ggTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTg                                                               >  1:169364/1‑251 (MQ=255)                                                                                                                                                                 gTCAGTTGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAACTCACCACTGc                                                              >  2:340707/1‑251 (MQ=255)                                                                                                                                                                 gTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGAAATCATTTTCACATCAATCGGATAGCGTTTACAGCTAGCTTCTAAATCACCACTGc                                                              >  1:382838/1‑251 (MQ=255)                                                                                                                                                                   cAGATTATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTTGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTa                                                            <  2:471251/251‑1 (MQ=255)                                                                                                                                                                      atgatgatAATTTTCAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCATTGCCTGATTTCAGCCTGGCATTACCTTCCATTGCGCCACTAAAGTCATTAATGGCGTAAAAATGAACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTaca                                                          >  2:428714/1‑248 (MQ=255)                                                                                                                                                                         atgatAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCCTTAATGGCGTAAAAATGc                                                                                                                                                                                       >  2:114244/1‑122 (MQ=255)                                                                                                                                                                         atgatAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGc                                                                                                                                                                                       <  1:114244/122‑1 (MQ=255)                                                                                                                                                                         atgatAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAag                                                        >  1:59712/1‑249 (MQ=255)                                                                                                                                                                                   ttAGTATATTCTATAGTCACG‑CCATGCTTCATGCGGACTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAACAATTCAATTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAATTGTTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTGGCTTCTCAATCACCACTGCTATAAGAGACggcgta                                              >  2:349500/1‑244 (MQ=255)                                                                                                                                                                                               ataGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTgg                                >  1:272943/1‑251 (MQ=255)                                                                                                                                                                                                 aGTCACG‑CTATGCTTCCTGCGTAATTATATNCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGt                               <  2:245349/250‑1 (MQ=255)                                                                                                                                                                                                           aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTaa                      <  1:468912/250‑1 (MQ=255)                                                                                                                                                                                                            ttcTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAACGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAACGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAAt                     >  2:148888/3‑250 (MQ=255)                                                                                                                                                                                                                              atCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGAAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAAc  >  1:320766/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             GACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACG‑CTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAAC  >  NC_000913/155229‑155699

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