Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,661,816	T→C	F421L (TTC→CTC)	ynfF →	putative selenate reductase YnfF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,661,816	0	T	C	100.0%	238.7 / NA	67	F421L (TTC→CTC)	ynfF	putative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/36);  total (31/36)

GGGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGT  >  NC_000913/1661679‑1661951                                                                                                                                          |                                                                                                                                        gggagCCACAACGACATTCCACCGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGCCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTc                                                                                                                                       >  2:369552/1‑140 (MQ=255)  ggggCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTcc                                                                                                                                      >  2:528964/1‑140 (MQ=255)  ggggCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTcc                                                                                                                                      >  1:703414/1‑140 (MQ=255)    ggCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCg                                                                                                                                    <  2:1012488/140‑1 (MQ=255)     gCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGa                                                                                                                                   <  1:1021307/140‑1 (MQ=255)       cacaACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATg                                                                                                                                 >  2:336664/1‑140 (MQ=255)       cacaACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATg                                                                                                                                 >  1:898081/1‑140 (MQ=255)         caACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCt                                                                                                                               <  1:62799/140‑1 (MQ=255)                      aCGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGtt                                                                                                                  >  2:368004/1‑140 (MQ=255)                       cGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTa                                                                                                                 >  1:389638/1‑140 (MQ=255)                         gCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTaaa                                                                                                               >  1:259350/1‑140 (MQ=255)                           gAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAAc                                                                                                             >  2:194134/1‑140 (MQ=255)                           gAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAAc                                                                                                             >  1:385057/1‑140 (MQ=255)                            aGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACg                                                                                                            >  1:341487/1‑140 (MQ=255)                            aGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCACGATGCTCGAGAATCCTGTTAAAACg                                                                                                            >  1:849172/1‑140 (MQ=255)                                  cATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGAtt                                                                                                      <  1:1208916/140‑1 (MQ=255)                                   aTCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGAttt                                                                                                     <  2:1206061/140‑1 (MQ=255)                                       cgcgCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGt                                                                                                 >  2:924441/1‑140 (MQ=255)                                        gcgcTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTc                                                                                                <  1:336664/140‑1 (MQ=255)                                        gcgcTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTc                                                                                                <  1:194134/140‑1 (MQ=255)                                        gcgcTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTc                                                                                                <  1:528964/140‑1 (MQ=255)                                        gcgcTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTc                                                                                                <  1:368004/140‑1 (MQ=255)                                         cgcTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCt                                                                                               >  1:733486/1‑140 (MQ=255)                                               tGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACAt                                                                                         >  1:718905/1‑140 (MQ=255)                                                    tGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACa                                                                                    >  2:799123/1‑140 (MQ=255)                                                    tGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACa                                                                                    >  2:1189136/1‑140 (MQ=255)                                                      cTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGa                                                                                  <  1:1189136/140‑1 (MQ=255)                                                      cTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGa                                                                                  <  2:898081/140‑1 (MQ=255)                                                           cGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGcc                                                                              <  2:729867/139‑1 (MQ=255)                                                           cGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGcc                                                                              >  1:729867/1‑139 (MQ=255)                                                            gTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCAt                                                                            <  1:1180075/140‑1 (MQ=255)                                                               ctcACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGa                                                                         <  1:924441/140‑1 (MQ=255)                                                                tcACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGAc                                                                        <  1:799123/140‑1 (MQ=255)                                                                tcACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGAc                                                                        <  2:703414/140‑1 (MQ=255)                                                                  aCCGGCAACGTCGGCATAAACTGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCa                                                                      >  2:862633/1‑140 (MQ=255)                                                                  aCCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCa                                                                      <  2:849172/140‑1 (MQ=255)                                                                  aCCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCa                                                                      <  2:259350/140‑1 (MQ=255)                                                                  aCCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCa                                                                      <  2:946227/140‑1 (MQ=255)                                                                  aCCGGCAACGTC‑GCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCa                                                                      <  2:1035544/139‑1 (MQ=255)                                                                  aCCGGCAACGTC‑GCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCa                                                                      >  1:1035544/1‑139 (MQ=255)                                                                    cGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATg                                                                    <  2:954167/140‑1 (MQ=255)                                                                     ggCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATgg                                                                   >  1:317446/1‑140 (MQ=255)                                                                      gCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGt                                                                  <  2:1207374/140‑1 (MQ=255)                                                                      gCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGt                                                                  <  2:389638/140‑1 (MQ=255)                                                                             cGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGaaa                                                           <  2:907659/140‑1 (MQ=255)                                                                             cGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGaaa                                                           <  1:862633/140‑1 (MQ=255)                                                                                   aaaCGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACcg                                                     <  1:576391/140‑1 (MQ=255)                                                                                     aCGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATGCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACcg                                                     >  2:845618/1‑138 (MQ=255)                                                                                     aCGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGt                                                                  >  2:718115/1‑125 (MQ=255)                                                                                     aCGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGt                                                                  <  1:718115/125‑1 (MQ=255)                                                                                     aCGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACcg                                                     <  1:845618/138‑1 (MQ=255)                                                                                       ggcggcAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGAc                                                 <  1:199114/140‑1 (MQ=255)                                                                                       ggcggcAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGAc                                                 <  1:526797/140‑1 (MQ=255)                                                                                         cggcAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGAccc                                               >  1:437946/1‑140 (MQ=255)                                                                                         cggcAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAA‑GGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCg                                              >  2:1059192/1‑140 (MQ=255)                                                                                          ggcAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCg                                              >  2:231365/1‑140 (MQ=255)                                                                                              aCTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGAt                                          <  2:790762/140‑1 (MQ=255)                                                                                                      gTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCg                                  <  1:369552/140‑1 (MQ=255)                                                                                                         cgcgAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTgg                               >  2:903689/1‑140 (MQ=255)                                                                                                           cgAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGgaa                             <  2:317446/140‑1 (MQ=255)                                                                                                            gAAGGTAGCTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGgaaa                            >  1:909874/1‑140 (MQ=255)                                                                                                                    cTGGGATCTGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAAc                    >  1:837123/1‑140 (MQ=255)                                                                                                                            ttggggTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTc            >  1:559173/3‑140 (MQ=255)                                                                                                                            tGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTc            >  2:133822/1‑140 (MQ=255)                                                                                                                            tGGGGGTAGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTc            <  1:186062/140‑1 (MQ=255)                                                                                                                              ggggTAGAATGGCTCCCGATGCTCGAGAATCCTGCTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTccc          >  2:1078245/1‑140 (MQ=255)                                                                                                                                   aGAATGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCa     <  2:718905/140‑1 (MQ=255)                                                                                                                                      aTGGCTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGt  <  1:231365/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        GGGGGCCACAACGACATTCCAACGGCGAGCAAACATCCCGCGCTATTGCCATGCTTTCCGTTCTCACCGGCAACGTCGGCATAAACGGCGGCAACTCAGGCGTACGCGAAGGTAGCTGGGATCTGGGGGTAGAATGGTTCCCGATGCTCGAGAATCCTGTTAAAACGCAGATTTCCGTCTTTACATGGACAGATGCCATCGACCATGGTACGGAAATGACCGCGACCCGTGATGGTGTTCGTGGAAAAGAAAAACTGGATGTCCCCATCAAGT  >  NC_000913/1661679‑1661951

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