Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,771,979	G→T	intergenic (‑235/‑302)	yibI ← / → mtlA	DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,771,979	0	G	T	100.0%	149.1 / NA	42	intergenic (‑235/‑302)	yibI/mtlA	DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/18);  total (24/18)

TTAGAAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCAC  >  NC_000913/3771735‑3772219                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 ttAGAAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTca                                                                                                                                                                                                                                             <  1:241702/250‑1 (MQ=255)       aCATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGTCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAAcc                                                                                                                                                                                                                                      <  1:184640/251‑1 (MQ=255)       aCATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAAGATCTTGTTATTCAGATCACAAATATTCAACAAAcc                                                                                                                                                                                                                                      <  2:130378/251‑1 (MQ=255)             cTGCCTCGGCTTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATcaa                                                                                                                                                                                                                                 >  2:274229/1‑250 (MQ=255)              gtccgtgggtgcACGATTATCATTGTTCTTTTTGGGCTGGCTGCCATAAGGGATTCAGGTTAATATTTTTACAATAATCCGGGTTCGCTATTCTTGCCATAACACCAAAGAATAATTTTTAGAGGTGATGTGTTGCTTAGTTACATAACGATTGTATGACGTAGGCATAACATTCTGTAGATCACATCACTTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACCAATATTCAACAACCCAGcaaaa                                                                                                                                                                                                                               <  2:723/240‑1 (MQ=255)                            gTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGacac                                                                                                                                                                                                                 <  2:233168/251‑1 (MQ=255)                                   aGTGTTGTTTTTAGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTGGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTcac                                                                                                                                                                                                          <  2:95587/251‑1 (MQ=255)                                     tgttgtTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGGGTTGCTTAGATACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGCTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTcacgt                                                                                                                                                                                                        >  1:273551/1‑249 (MQ=255)                                              tGGGCTGGCAGCCAGCAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACCCCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGcc                                                                                                                                                                                               >  2:117422/1‑251 (MQ=255)                                              tGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGcc                                                                                                                                                                                               >  1:273031/1‑251 (MQ=255)                                                             aaGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAAt                                                                                                                                                                                >  1:206687/1‑251 (MQ=255)                                                               gggAGGCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTGTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGAGTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGGGACGCTACTCACATTTAAATGCCATTTTTAGCTAAAATcg                                                                                                                                                                              <  1:19159/251‑1 (MQ=255)                                                                           gATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTTTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCTCATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttgtt                                                                                                                                                                    <  2:231361/249‑1 (MQ=255)                                                                                  tGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCtttttt                                                                                                                                                            >  1:27491/1‑250 (MQ=255)                                                                                                       aTTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGCATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGACCGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAACTGTGACACCACTCACATTTAACTGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGCTGAac                                                                                                                                      >  1:271187/1‑251 (MQ=255)                                                                                                            cGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATTTTCAACAAACCATCAAAACAAAAATGTGACGCTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtc                                                                                                                                   <  1:244392/249‑1 (MQ=255)                                                                                                            cGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtca                                                                                                                                  <  1:141960/250‑1 (MQ=255)                                                                                                                  aaCACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATcaa                                                                                                                                                                                                                                 >  2:219181/1‑149 (MQ=255)                                                                                                                  aaCACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGAGTCAGATCACAAATATTCAACAAACCATcaa                                                                                                                                                                                                                                 <  1:219181/149‑1 (MQ=255)                                                                                                                      ccAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGATGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAATCAAAAATGTGCCACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGATCGTCACGTCAATTAcc                                                                                                                       >  2:263010/1‑251 (MQ=255)                                                                                                                         aaGAATAATTTGTAGAGGATATGGGTTGCTTAGTTACAAAAAGATGGAATGACGAAGGCATAACATGATGTAGATAAAATCAGGTGAACGCCGTAAGAAAATCTCTTGTGATTCAGATCACATATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTCAATGCCATTTTAAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGGTGAACGCAACGCCAATTACctct                                                                                                                    >  2:249907/1‑251 (MQ=255)                                                                                                                          aGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctctc                                                                                                                   >  2:147176/1‑251 (MQ=255)                                                                                                                                  ttgtAGAGGTGATGAGTTGCTGAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTTCTTTTTGACACAAGCGTTTTGTGATGATCGTCACGTCAATTACCTCTCGAccccc                                                                                                            <  2:271788/247‑1 (MQ=255)                                                                                                                                   tttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCt                                                                                                           <  1:28072/250‑1 (MQ=255)                                                                                                                                     tAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCtat                                                                                                         <  1:117422/250‑1 (MQ=255)                                                                                                                                        agGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATAttt                                                                                                     >  2:214892/1‑251 (MQ=255)                                                                                                                                            gatgaGTTGCTTCTTTACATAACGATTGTATGACGACGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTAtgt                                                                                                 >  2:56242/1‑251 (MQ=255)                                                                                                                                            gatgaGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTAtg                                                                                                  <  1:274229/250‑1 (MQ=255)                                                                                                                                             atgaGTTGCTTCCTTACATACCGATGGTATGACGAAGGAATCACATGCTGTAGCTCACATCAGGTGAACGCCGTAAGCAAATATCTTGAGATTCAGATCAAAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACCAACATTAAAATGCAATTTTTAGCGAAAATCGCCGCATTGATGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTATCTACACCCTATATTTAtgtg                                                                                                >  2:154084/1‑251 (MQ=255)                                                                                                                                                 gTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTa                                                                                                                                                                                                     <  2:79322/146‑1 (MQ=255)                                                                                                                                                 gTTGCTCAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACCCTACTCACATTTa                                                                                                                                                                                                     >  1:79322/1‑146 (MQ=255)                                                                                                                                                                                                       aTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTc                                      >  1:38104/1‑251 (MQ=255)                                                                                                                                                                                                              gAACGCCGTAAGAAAATATCTTGTGATGCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTCGCGAAAATCGCCGCCTTGTTACTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCACTTACCTCTCTACCCCCTCTATTTATGTGATTGCTATCACCCAAACGGCCGTCGACTGGACAGCTAACCGATTCAGTCCCAGATTTCGCAgcat                               >  2:110031/1‑248 (MQ=255)                                                                                                                                                                                                                    cttAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACa                          >  1:69165/3‑250 (MQ=255)                                                                                                                                                                                                                     gTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACa                          >  2:84051/1‑249 (MQ=255)                                                                                                                                                                                                                          aaaaTATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTCACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCg                   >  1:233750/1‑251 (MQ=255)                                                                                                                                                                                                                            aaTATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGc                 >  1:2264/1‑251 (MQ=255)                                                                                                                                                                                                                                tCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTAcc             >  1:178601/1‑251 (MQ=255)                                                                                                                                                                                                                                   tgtgATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACct            <  2:27491/249‑1 (MQ=255)                                                                                                                                                                                                                                       aTTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTgccg      >  1:138702/1‑251 (MQ=255)                                                                                                                                                                                                                                            gATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTCCCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCcac  <  2:60533/250‑1 (MQ=255)                                                                                                                                                                                                                                            gATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCcac  >  1:60533/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 TTAGAAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCAC  >  NC_000913/3771735‑3772219

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