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%	256.2 / NA	61	T58M (ACG→ATG)	yicC	UPF0701 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/39);  total (22/39)

GGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGG  >  NC_000913/3816704‑3816994                                                                                                                                                 |                                                                                                                                                   ggCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCgg                                                                                                                                              >  1:201946/1‑151 (MQ=255)   cGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTa                                                                                                                                            <  2:55656/151‑1 (MQ=255)      gAAATCAAGGTTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAg                                                                                                                                         <  1:911090/151‑1 (MQ=255)           cAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGaa                                                                                                                                    <  2:654181/151‑1 (MQ=255)           cAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGa                                                                                                                                     <  2:966415/150‑1 (MQ=255)                  gAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTAcc                                                                                                                              <  2:182722/150‑1 (MQ=255)                    aTGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTg                                                                                                                           <  2:787930/151‑1 (MQ=255)                     tGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTg                                                                                                                           >  1:630198/1‑150 (MQ=255)                      gggggAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTgc                                                                                                                          <  2:383834/150‑1 (MQ=255)                      gggggAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTgc                                                                                                                          <  2:308629/150‑1 (MQ=255)                          gAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTAt                                                                                                                     <  2:588302/151‑1 (MQ=255)                          gAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTAt                                                                                                                     <  1:451378/151‑1 (MQ=255)                            gcgcAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGa                                                                                                                   <  2:962554/151‑1 (MQ=255)                                   cTGGGAAATGCGCTCGGTAAATCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGcc                                                                                                                >  1:662927/1‑147 (MQ=255)                                        aaaTGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAg                                                                                                       <  2:602477/151‑1 (MQ=255)                                               cTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCa                                                                                                 >  1:929790/1‑150 (MQ=255)                                                tCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCa                                                                                                 >  2:723451/1‑149 (MQ=255)                                                   gTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTg                                                                                            <  1:1015970/151‑1 (MQ=255)                                                       aCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGc                                                                                         >  1:134521/1‑150 (MQ=255)                                                       aCCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCt                                                                                        <  2:750493/151‑1 (MQ=255)                                                         cAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTg                                                                                       >  1:193834/1‑150 (MQ=255)                                                         cAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGa                                                                                      <  2:193829/151‑1 (MQ=255)                                                             gTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCt                                                                                  <  1:272265/151‑1 (MQ=255)                                                               tATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTc                                                                                 >  2:474128/1‑150 (MQ=255)                                                                  cTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACg                                                                             >  2:239894/1‑151 (MQ=255)                                                                   tGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACg                                                                             >  1:191403/1‑150 (MQ=255)                                                                    ggAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCaa                                                                               <  2:460614/147‑1 (MQ=255)                                                                    ggAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCaa                                                                               <  2:392809/147‑1 (MQ=255)                                                                           tacttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTgg                                                                    >  1:589972/1‑151 (MQ=255)                                                                           tacttTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTgg                                                                    >  1:475847/1‑151 (MQ=255)                                                                                 cGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTaaa                                                               <  1:218210/150‑1 (MQ=255)                                                                                 cGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAAc                                                              <  2:684648/151‑1 (MQ=255)                                                                                 cGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCCGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCt                                                                  <  2:769457/147‑1 (MQ=255)                                                                                  gTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACa                                                             <  1:128413/151‑1 (MQ=255)                                                                                      gCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCt                                                          <  1:346935/150‑1 (MQ=255)                                                                                         ggAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTa                                                      <  2:433565/151‑1 (MQ=255)                                                                                           aGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAAc                                                    >  1:911628/1‑151 (MQ=255)                                                                                             cAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTg                                                  <  2:415541/151‑1 (MQ=255)                                                                                                 tCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcgc                                              >  1:897989/1‑151 (MQ=255)                                                                                                  ccGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCcgcg                                             <  1:529388/151‑1 (MQ=255)                                                                                                           tGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTggg                                      <  1:407705/149‑1 (MQ=255)                                                                                                           tGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGt                                     <  1:467273/150‑1 (MQ=255)                                                                                                            gAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaa                                   <  2:798543/151‑1 (MQ=255)                                                                                                            gAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTa                                    >  2:136785/1‑150 (MQ=255)                                                                                                               ccTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTaaaaa                                <  2:75210/151‑1 (MQ=255)                                                                                                                 tGTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATg                              >  2:1019098/1‑151 (MQ=255)                                                                                                                  gTCGTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGc                             <  2:761350/151‑1 (MQ=255)                                                                                                                     gTTCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGc                             >  1:735589/1‑148 (MQ=255)                                                                                                                       tCGCGAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCagag                         <  2:201941/150‑1 (MQ=255)                                                                                                                        cgcgAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGt                        <  2:451389/150‑1 (MQ=255)                                                                                                                        cgcgAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTg                       <  1:793746/151‑1 (MQ=255)                                                                                                                           gAGCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACg                    >  1:316830/1‑151 (MQ=255)                                                                                                                             gCGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGaa                  <  1:576868/151‑1 (MQ=255)                                                                                                                              cGTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAg                 <  1:246104/151‑1 (MQ=255)                                                                                                                               gTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAgg                >  1:962882/1‑151 (MQ=255)                                                                                                                               gTATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAgg                >  1:822933/1‑151 (MQ=255)                                                                                                                                tATTCGTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAgg                <  1:729976/150‑1 (MQ=255)                                                                                                                                     gTTCTCGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAAt          <  2:990365/151‑1 (MQ=255)                                                                                                                                         tcGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAAc      >  2:914515/1‑151 (MQ=255)                                                                                                                                          cGCCTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCaa       <  1:720768/149‑1 (MQ=255)                                                                                                                                             cTGATGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCgg  >  2:630884/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAACCAGCGTTATCTGGAAACTTACTTTCGTCTGCCGGAGCAGTTCCGTAGCCTTGAACCTGTCGTTCGCGAGCGTATTCGTTCTCGCCTGACGCGCGGTAAAGTGGAATGTACCCTGCGCTATGAGCCAGATGTTAGCGCGCAAGGTGAGCTGATCCTCAACGAAAAACTGGCTAAACAGCTGGTAACTGCCGCGAACTGGGTAAAAATGCAGAGTGACGAAGGGGAAATCAACCCGG  >  NC_000913/3816704‑3816994

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