Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,947,489	A→C	L160W (TTG→TGG)	hdfR ←	flhDC operon transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,947,489	0	A	C	100.0%	250.1 / NA	65	L160W (TTG→TGG)	hdfR	flhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (30/35);  total (30/35)

TCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATAC  >  NC_000913/3947346‑3947635                                                                                                                                                |                                                                                                                                                   tCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGact                                                                                                                                             >  2:1021478/1‑151 (MQ=255)  cAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGacta                                                                                                                                            <  1:1073706/151‑1 (MQ=255)   aGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTAAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCTCTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGact                                                                                                                                             <  2:860111/149‑1 (MQ=255)           cACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCg                                                                                                                                   <  1:101660/151‑1 (MQ=255)           cACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCg                                                                                                                                   >  2:268438/1‑151 (MQ=255)              ttCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGGCGCAGATAATTAAGATCGCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTc                                                                                                                                 >  2:535181/1‑150 (MQ=255)               tCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTcc                                                                                                                                <  2:829191/150‑1 (MQ=255)               tCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTc                                                                                                                                 <  1:256069/149‑1 (MQ=255)                cGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTc                                                                                                                                 <  2:1019944/148‑1 (MQ=255)                cGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTc                                                                                                                                 <  2:546353/148‑1 (MQ=255)                   cAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCAtttt                                                                                                                           <  1:608086/151‑1 (MQ=255)                   cAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCAtttt                                                                                                                           >  2:370832/1‑151 (MQ=255)                   cAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCAttt                                                                                                                            <  2:470246/150‑1 (MQ=255)                     gCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTg                                                                                                                          >  1:454032/1‑150 (MQ=255)                      cACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTggg                                                                                                                        <  2:1052663/151‑1 (MQ=255)                       aCCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGgc                                                                                                                       >  1:132426/1‑151 (MQ=255)                        ccGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGgcg                                                                                                                      <  1:370918/151‑1 (MQ=255)                             cAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTc                                                                                                                  >  1:403689/1‑150 (MQ=255)                             cAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTc                                                                                                                  <  1:196540/150‑1 (MQ=255)                                    gCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCGTTTTGGGCGCTTCAGTGGta                                                                                                           <  1:486488/150‑1 (MQ=255)                                     ccTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAATCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGtaat                                                                                                         <  2:906690/151‑1 (MQ=255)                                     ccTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGtaat                                                                                                         >  1:582765/1‑151 (MQ=255)                                        cATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGtaata                                                                                                        <  1:129620/149‑1 (MQ=255)                                        cATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGtaata                                                                                                        >  1:78367/1‑149 (MQ=255)                                        cATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATaag                                                                                                      <  1:132158/151‑1 (MQ=255)                                             tGTTGAAAATCTGGCCCCCACTCAAGTCGCAGCTAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATaaga                                                                                                     >  1:275629/1‑147 (MQ=255)                                                 gAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCaa                                                                                              <  2:677338/150‑1 (MQ=255)                                                 gAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCaa                                                                                              <  2:32584/150‑1 (MQ=255)                                                 gAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAg                                                                                             <  1:858453/151‑1 (MQ=255)                                                  aaaaTCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGc                                                                                            <  1:607562/151‑1 (MQ=255)                                                        tGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCGTCAGTGGTAATAAGAAGATCAAGCTGGCGt                                                                                      >  1:661027/1‑151 (MQ=255)                                                         ggCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGtt                                                                                     >  2:557948/1‑151 (MQ=255)                                                         ggCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGtt                                                                                     >  2:912029/1‑151 (MQ=255)                                                         ggCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGtt                                                                                     >  1:365424/1‑151 (MQ=255)                                                          gCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTc                                                                                    <  1:359454/151‑1 (MQ=255)                                                               cACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTAGCGTTCATGc                                                                                >  1:809321/1‑150 (MQ=255)                                                                   cAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTg                                                                           >  2:958857/1‑151 (MQ=255)                                                                     aGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGtt                                                                         <  2:868471/151‑1 (MQ=255)                                                                      gTCGCAGTTAATTAAGATCTCCCTTTAGTTTTGAAGCGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGttt                                                                        >  2:835236/1‑151 (MQ=255)                                                                      gTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGttt                                                                        <  2:880199/151‑1 (MQ=255)                                                                       tCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGttt                                                                        <  1:665268/150‑1 (MQ=255)                                                                         gCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGCGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTAc                                                                     >  1:180794/1‑151 (MQ=255)                                                                          cAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTAcc                                                                    <  1:762366/151‑1 (MQ=255)                                                                           aGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCa                                                                   <  1:405273/151‑1 (MQ=255)                                                                             aTAATTAAGATCTCCCTTGAGTTTTGAAGGGGCACTGGTATAAAGCGCGAAAGTGAAATATCCCAACCACTGACTACTAAATTCGTCCATTTGGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCag                                                                  <  2:1082492/150‑1 (MQ=255)                                                                              tAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGTCGTTCATGCAGCTGTTTTAcc                                                                    >  2:247451/1‑147 (MQ=255)                                                                               aaTTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCagag                                                                <  1:1038259/150‑1 (MQ=255)                                                                                     gATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCg                                                         >  2:858787/1‑151 (MQ=255)                                                                                          cccTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGt                                                        <  1:945759/147‑1 (MQ=255)                                                                                            cTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGc                                                   >  1:1046439/1‑150 (MQ=255)                                                                                             tttAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCa                                                  >  2:263923/1‑150 (MQ=255)                                                                                               tAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGAGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGCTTTACCAGAGACTGCCGTTGGGCAAtt                                               >  1:221796/1‑151 (MQ=255)                                                                                                     tGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCgcgc                                          <  1:114791/150‑1 (MQ=255)                                                                                                              aCTGGTATAAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGt                                >  1:1081200/1‑151 (MQ=255)                                                                                                                   tataAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGcc                           >  2:266533/1‑151 (MQ=255)                                                                                                                     taAAGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCt                          >  1:1033851/1‑150 (MQ=255)                                                                                                                      aaaGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGt                        >  1:902883/1‑151 (MQ=255)                                                                                                                      aaaGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGt                        <  1:1148086/151‑1 (MQ=255)                                                                                                                        aGCGCTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTAt                      <  1:566667/151‑1 (MQ=255)                                                                                                                         gcgcTAAAGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATg                     <  2:409978/151‑1 (MQ=255)                                                                                                                              aaaGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCa                 <  2:80108/150‑1 (MQ=255)                                                                                                                              aaaGTGAAATATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCAt                >  2:71279/1‑151 (MQ=255)                                                                                                                                     aaTATCCCAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGAtt         <  2:922495/151‑1 (MQ=255)                                                                                                                                           ccAGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGAt    >  1:256521/1‑150 (MQ=255)                                                                                                                                             aGCCACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATAc  >  2:57412/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                   TCAGAATGGGCACTTCGTCAGCACCGATCAAACCTGCCTCATGCTGTTGAAAATCTGGCCCCCACTCAAGTCGCAGATAATTAAGATCTCCCTTTAGTTTTGAAGGGGCACTGGTATAAAGCGCTAAAGTGAAATATCCCAGCAACTGACTACTAAATTCGTCCATTTTGGGCGCTTCAGTGGTAATAAGAAGATCAAGCTGGCGTTCATGCAGCTGTTTTACCAGAGACTGCCGTTGGGCAATTCGCGCTTCGAACTGTAAGCCTGTATGGGCATCCTGATTTTGATAC  >  NC_000913/3947346‑3947635

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