Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,744,369	C→T	A14A (GCC→GCT)	yiaM →	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,744,369	0	C	T	100.0%	142.7 / NA	42	A14A (GCC→GCT)	yiaM	2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/20);  total (22/20)

CACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCTTATTTAT  >  NC_000913/3744130‑3744599                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                       cACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCtg                                                                                                                                                                                                                             >  1:137437/1‑251 (MQ=255)        cAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCtgtattgt                                                                                                                                                                                                                       >  1:48261/1‑250 (MQ=255)          ggTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTgtattgtatt                                                                                                                                                                                                                    >  2:544806/1‑251 (MQ=255)          ggTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTgtattgtatt                                                                                                                                                                                                                    >  1:278787/1‑251 (MQ=255)              gTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAAGCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTtat                                                                                                                                                                                                                 <  2:10437/250‑1 (MQ=255)                attatGCAAACATCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATaa                                                                                                                                                                                                               <  2:100999/250‑1 (MQ=255)                                                          gTCGCATTAACTGCGATGAATTAATCGAATAATCGTCTACGACAGAACGACTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGATGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGt                                                                                                                                                                    <  1:241605/251‑1 (MQ=255)                                                            cGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTtg                                                                                                                                                                  >  1:509297/1‑251 (MQ=255)                                                             gCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTtga                                                                                                                                                                 <  1:281710/251‑1 (MQ=255)                                                                                 aTCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTAGTGCTCTAATTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCAGTTTATCTGTTGATGAATTGTCACGttatttat                                                                                                                                             <  2:549510/251‑1 (MQ=255)                                                                                   cGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTtatttatt                                                                                                                                            <  1:544806/250‑1 (MQ=255)                                                                                           cGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGt                                                                                                                                    <  1:426322/250‑1 (MQ=255)                                                                                                                               tAACTATCTTAATAATTCTGGATGCTACACATGAGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACCTAATTTTAAGCTATGGTTTTCAGACAAGCATTTTATCTGTTGATGACTTGTCACGATATTTCTTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATa                                                                                                >  2:22597/1‑250 (MQ=255)                                                                                                                                aaCTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAAc                                                                                              >  1:415905/1‑251 (MQ=255)                                                                                                                                              ttCAGGTTGCTGCACATGCGACTTGAATTATCACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTAAGACAAGCATTTTATCTGTTGATGAATTGTAACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTCGCTTTTATGGTTAACGCacactctacggt                                                                                >  2:437831/1‑243 (MQ=255)                                                                                                                                                cAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGCACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTc                                                                               >  1:467639/1‑250 (MQ=255)                                                                                                                                                 aGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCAc                                                                             >  1:350720/1‑251 (MQ=255)                                                                                                                                                   gTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCt                                                                           >  1:17661/1‑251 (MQ=255)                                                                                                                                                     ttctgcAAATGCGACTTGAATTATGACAGACATAGCCTCAAGTAATAGCTATGAAAAAAATTCTCGAAGCAATTCTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTGATCTGTTGATGAATTGTCACGTTAGTTATTTGTCTGGTTAACGTTTATTGGCGAGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCtt                                                                          <  2:326974/248‑1 (MQ=255)                                                                                                                                                              tGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACattatt                                                                                                                                                                                                       >  1:162081/1‑116 (MQ=255)                                                                                                                                                              tGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACattatt                                                                                                                                                                                                       <  2:162081/116‑1 (MQ=255)                                                                                                                                                               gCGACTTGAATTATGACAGACATAGCCTCAATGAATAGCTATGAAAAAAATTCTCGAAGCAATACATGCGATTAATCTCGCTGTACTTGCCTTTATTTTTTCTATAAACATTATTGTaa                                                                                                                                                                                                   <  2:515455/119‑1 (MQ=255)                                                                                                                                                               gCGACGTGAATTATGACAGACATAGCCACAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCAGTATTGTATTTATAAACATTATTTTaa                                                                                                                                                                                                   >  1:515455/1‑119 (MQ=255)                                                                                                                                                                     tGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAActct                                                         >  2:14138/1‑251 (MQ=255)                                                                                                                                                                           aTGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCg                                                   >  1:536503/1‑251 (MQ=255)                                                                                                                                                                           aTGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCg                                                   <  1:336496/251‑1 (MQ=255)                                                                                                                                                                              acagacaTAGACTCAAGGAATAGCTATGAAAAACATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCACAAGTTCAGGTCACCTTTCTGGGGGAGAAActcatcaccgcat                                                >  2:371223/1‑241 (MQ=255)                                                                                                                                                                                        ccTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGcgacg                                      >  1:377100/1‑251 (MQ=255)                                                                                                                                                                                            aaGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTc                                  <  1:14138/251‑1 (MQ=255)                                                                                                                                                                                             aGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCg                                 >  1:321189/1‑251 (MQ=255)                                                                                                                                                                                                  tAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGATTCTCTCTg                            >  2:440666/1‑251 (MQ=255)                                                                                                                                                                                                       aTGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGtt                         <  2:137437/249‑1 (MQ=255)                                                                                                                                                                                                       aTGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTAc                       <  2:498190/251‑1 (MQ=255)                                                                                                                                                                                                         gAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTAccc                     <  1:36462/251‑1 (MQ=255)                                                                                                                                                                                                          aaaaaaaTACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTAccc                     >  1:455953/1‑250 (MQ=255)                                                                                                                                                                                                                aTACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCtt              <  1:191708/251‑1 (MQ=255)                                                                                                                                                                                                                aTACTCGAAGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTAACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTTTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGGGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCtt              >  2:142499/1‑251 (MQ=255)                                                                                                                                                                                                                        aGCAATACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCttat      <  2:468262/251‑1 (MQ=255)                                                                                                                                                                                                                           aaTACTGGCGATTAATCTCTCTTGACTTACCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACTATCATTTTATTTGTTGATGAATTGTTACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTATCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCTTTCCTTTATCttattta   <  2:424835/251‑1 (MQ=255)                                                                                                                                                                                                                            aTACTGGCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCttatttat  <  1:254170/251‑1 (MQ=255)                                                                                                                                                                                                                                  gCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCttt                                                                         >  1:53112/1‑174 (MQ=255)                                                                                                                                                                                                                                  gCGATTAATCTCGCTGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCttt                                                                         <  2:53112/174‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                       CACCGACCAGGTTGTATTATGCAAACAGCCTCTGAGATACGAAAAATTGTGGTTAAAGTCGCATTAACGGCGCTGAATTAATCGAATAATCGTCTACGCCAGAACGCCTGGGTAATGTTATTGCTCTAACTATATTAATAATTCAGGTTGCTGCACATGCGACTTGAATTATGACAGACATAGCCTCAAGGAATAGCTATGAAAAAAATACTCGAAGCAATACTGGCGATTAATCTCGCCGTACTTTCCTGTATTGTATTTATAAACATTATTTTAAGATATGGTTTTCAGACAAGCATTTTATCTGTTGATGAATTGTCACGTTATTTATTTGTCTGGTTAACGTTTATTGGCGCGATTGTAGCTTTTATGGATAACGCCCACGTTCAGGTCACCTTTCTGGTGGAGAAACTCTCCCCCGCATGGCAGCGACGAGTCGCTCTGGTTACCCATTCCTTAATCTTATTTAT  >  NC_000913/3744130‑3744599

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: ‑