Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,848	C→T	T58M (ACG→ATG)	yicC →	UPF0701 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,848	0	C	T	100.0%	241.6 / NA	60	T58M (ACG→ATG)	yicC	UPF0701 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/25);  total (35/25)

CCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCC  >  NC_000913/3816701‑3816992                                                                                                                                                    |                                                                                                                                                 cccGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGAtgcg                                                                                                                                               >  1:582075/1‑151 (MQ=255)  ccGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGAtgcgc                                                                                                                                              <  1:231129/151‑1 (MQ=255)      cGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGt                                                                                                                                           >  1:737411/1‑150 (MQ=255)      cGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTa                                                                                                                                          <  1:854889/151‑1 (MQ=255)       gTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTaa                                                                                                                                         <  2:936220/151‑1 (MQ=255)         gAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAg                                                                                                                                       <  1:388676/151‑1 (MQ=255)              cAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGaa                                                                                                                                  <  2:68446/151‑1 (MQ=255)               aaGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAAt                                                                                                                                 <  2:765941/151‑1 (MQ=255)                aGGGTGAATGGGGGAGCGCAACCTGGGAAATGCACTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAAtt                                                                                                                                <  1:129871/151‑2 (MQ=255)                   gTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTAc                                                                                                                             >  2:932738/1‑151 (MQ=255)                    tGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTAcc                                                                                                                            >  1:346011/1‑151 (MQ=255)                      aaTGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCt                                                                                                                          >  1:1006757/1‑151 (MQ=255)                       aTGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTg                                                                                                                         <  1:559053/151‑1 (MQ=255)                           gggAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCt                                                                                                                     <  1:588447/151‑1 (MQ=255)                            ggAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTa                                                                                                                    <  1:945905/151‑1 (MQ=255)                                 gcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGc                                                                                                               >  1:239197/1‑151 (MQ=255)                                 gcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGc                                                                                                               >  2:264160/1‑151 (MQ=255)                                        gggAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGt                                                                                                        >  2:33908/1‑151 (MQ=255)                                               gcgcTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAgcgc                                                                                                  <  1:1121350/150‑1 (MQ=255)                                                cgcTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAgcgcgc                                                                                                >  2:129659/1‑151 (MQ=255)                                                   tCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCa                                                                                               >  1:513015/1‑149 (MQ=255)                                                      gTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTg                                                                                          <  2:1169440/151‑1 (MQ=255)                                                       tAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTg                                                                                          >  1:1110136/1‑150 (MQ=255)                                                       tAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTg                                                                                          >  1:1011608/1‑150 (MQ=255)                                                       tAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTACCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGTTGCGCGGTAAAGAGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAg                                                                                             >  1:396901/1‑147 (MQ=255)                                                        aaaCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTg                                                                                          >  1:44/1‑149 (MQ=255)                                                                gTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCt                                                                                >  2:990263/1‑151 (MQ=255)                                                                gTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCt                                                                                <  1:488273/151‑1 (MQ=255)                                                                     cTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACg                                                                           >  2:271173/1‑151 (MQ=255)                                                                       ggAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACg                                                                           <  2:439770/149‑1 (MQ=255)                                                                       ggAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGaa                                                                         <  1:833545/151‑1 (MQ=255)                                                                       ggAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGaa                                                                         <  1:397648/151‑1 (MQ=255)                                                                            cttacttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACt                                                                    >  2:391685/1‑151 (MQ=255)                                                                                 tttCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGGGAGCTGATCCTCAACGAAAAACTGGCt                                                                >  2:367221/1‑150 (MQ=255)                                                                                  ttCGTCCGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTa                                                               >  1:1119957/1‑150 (MQ=255)                                                                                         gCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCt                                                        <  1:1145630/150‑1 (MQ=255)                                                                                         gCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTg                                                       >  2:90645/1‑151 (MQ=255)                                                                                           cGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGt                                                     <  1:173331/151‑1 (MQ=255)                                                                                                cAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTg                                                >  2:1051682/1‑151 (MQ=255)                                                                                                 aGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGc                                               <  1:993921/151‑1 (MQ=255)                                                                                                  gTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGcc                                              >  1:177062/1‑151 (MQ=255)                                                                                                   ttCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAACGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcg                                             >  2:346988/1‑151 (MQ=255)                                                                                                         aGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACt                                       <  1:949740/151‑1 (MQ=255)                                                                                                                aaCCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaa                                 >  2:580354/1‑150 (MQ=255)                                                                                                                aaCCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaa                                 >  1:156215/1‑150 (MQ=255)                                                                                                                 aCCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaaaa                               <  2:798844/151‑1 (MQ=255)                                                                                                                  ccTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaaaaa                              <  1:1032014/151‑1 (MQ=255)                                                                                                                  ccTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaaaa                               >  1:161350/1‑150 (MQ=255)                                                                                                                   cTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAAt                             >  1:335351/1‑151 (MQ=255)                                                                                                                    tGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATg                            >  2:695288/1‑151 (MQ=255)                                                                                                                     gTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATg                            <  2:33742/150‑1 (MQ=255)                                                                                                                       cGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCag                         <  1:848488/151‑1 (MQ=255)                                                                                                                         ttCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCag                         >  1:373437/1‑149 (MQ=255)                                                                                                                           cgcgAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTg                     <  1:851100/151‑1 (MQ=255)                                                                                                                             cgAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGAc                   >  1:1064199/1‑151 (MQ=255)                                                                                                                                      tCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGa           >  2:687010/1‑150 (MQ=255)                                                                                                                                        gTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAAt        >  2:63183/1‑151 (MQ=255)                                                                                                                                          tctcGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATc       <  1:44573/150‑1 (MQ=255)                                                                                                                                          tctcGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATc       >  1:619970/1‑150 (MQ=255)                                                                                                                                           ctggCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCa      >  1:828099/4‑150 (MQ=255)                                                                                                                                             cGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAAcc   >  2:81427/1‑151 (MQ=255)                                                                                                                                              gCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAAccc  >  2:604859/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCC  >  NC_000913/3816701‑3816992

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