Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,400,811	G→T	H93N (CAC→AAC)	mreB ←	cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,400,811	0	G	T	100.0%	381.2 / NA	98	H93N (CAC→AAC)	mreB	cell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (49/49);  total (49/49)

TCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGAT  >  NC_000913/3400663‑3400957                                                                                                                                                     |                                                                                                                                                   tcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCt                                                                                                                                                  >  2:1043613/1‑151 (MQ=255) tcAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCt                                                                                                                                                  >  1:350681/1‑151 (MQ=255)   aaTCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTgg                                                                                                                                                >  2:252171/1‑151 (MQ=255)     tCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTgg                                                                                                                                                >  2:568153/1‑149 (MQ=255)      cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAg                                                                                                                                              <  1:200516/150‑1 (MQ=255)      cAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGc                                                                                                                                             >  1:847657/1‑151 (MQ=255)        ggAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCAt                                                                                                                                           <  2:1121680/151‑1 (MQ=255)        ggAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCAt                                                                                                                                           <  1:202433/151‑1 (MQ=255)          aaGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCAttt                                                                                                                                         <  2:668668/151‑1 (MQ=255)            gACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCAttttt                                                                                                                                       >  2:126135/1‑151 (MQ=255)             aCTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCAtttttt                                                                                                                                      <  2:238901/151‑1 (MQ=255)                 cACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAg                                                                                                                                   <  1:137813/150‑1 (MQ=255)                   cGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTc                                                                                                                                 >  1:586449/1‑150 (MQ=255)                     ggCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACg                                                                                                                              <  2:187721/151‑1 (MQ=255)                         ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACgaaga                                                                                                                          <  2:143353/151‑1 (MQ=255)                         ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACgaaga                                                                                                                          <  1:891918/151‑1 (MQ=255)                         ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACga                                                                                                                             >  1:663646/1‑148 (MQ=255)                           aGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATACAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACgaaga                                                                                                                          >  2:603235/1‑149 (MQ=255)                            gcgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGt                                                                                                                       <  2:701816/151‑1 (MQ=255)                             cgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTc                                                                                                                      <  2:161022/151‑1 (MQ=255)                             cgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTc                                                                                                                      <  1:389651/151‑1 (MQ=255)                               cccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCgg                                                                                                                    >  1:1089302/1‑151 (MQ=255)                                ccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGc                                                                                                                   >  1:89621/1‑151 (MQ=255)                                   gcgcGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGa                                                                                                                 <  1:746307/150‑1 (MQ=255)                                      cgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAAc                                                                                                             <  1:487465/151‑1 (MQ=255)                                       ggATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACg                                                                                                            >  2:674665/1‑151 (MQ=255)                                        gATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGc                                                                                                           >  1:54914/1‑151 (MQ=255)                                         aTTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGcc                                                                                                          >  1:614646/1‑151 (MQ=255)                                           tCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGt                                                                                                        >  1:562686/1‑151 (MQ=255)                                             aCGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCt                                                                                                      <  2:571156/151‑1 (MQ=255)                                               gAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCttt                                                                                                    >  2:865061/1‑151 (MQ=255)                                                  ttGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTTAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCAt                                                                                                 >  1:606555/1‑151 (MQ=255)                                                    gcgcGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTg                                                                                               >  1:1125175/1‑151 (MQ=255)                                                     cgcgGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTgg                                                                                              >  1:155212/1‑151 (MQ=255)                                                        ggcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCg                                                                                           <  1:1099469/151‑1 (MQ=255)                                                         gcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGACGTCGGCGATAACGCCGTCTTTCATTGGGCGa                                                                                          <  1:475820/151‑1 (MQ=255)                                                          cgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGa                                                                                          <  2:770123/150‑1 (MQ=255)                                                           gTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAAt                                                                                        >  1:186301/1‑151 (MQ=255)                                                           gTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAAt                                                                                        <  1:308146/151‑1 (MQ=255)                                                           gTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAAt                                                                                        <  2:1081207/151‑1 (MQ=255)                                                             tCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATgg                                                                                      >  1:636706/1‑151 (MQ=255)                                                               aaCCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgca                                                                                    <  1:1255891/151‑1 (MQ=255)                                                               aaCCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgc                                                                                     >  1:707331/1‑150 (MQ=255)                                                                 ccTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcagc                                                                                  <  2:923026/151‑1 (MQ=255)                                                                    gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGgcagca                                                                                 <  2:1132947/149‑1 (MQ=255)                                                                    gggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAat                                                                               >  2:661905/1‑151 (MQ=255)                                                                     ggTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAata                                                                              <  2:642634/151‑1 (MQ=255)                                                                        cgcgCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATAtt                                                                            <  1:363509/150‑1 (MQ=255)                                                                          cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGcc                                                                         >  2:124177/1‑151 (MQ=255)                                                                           gcCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGccc                                                                        >  1:758562/1‑151 (MQ=255)                                                                            ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCg                                                                       <  1:57035/151‑1 (MQ=255)                                                                              aaCCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGc                                                                     <  2:72460/151‑1 (MQ=255)                                                                              aaCCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGc                                                                     >  1:1071149/1‑151 (MQ=255)                                                                                ccGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGt                                                                   >  2:423096/1‑151 (MQ=255)                                                                                 cGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTa                                                                  <  2:1017066/151‑1 (MQ=255)                                                                                    cacacaAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACgg                                                               <  1:163047/151‑1 (MQ=255)                                                                                    cacacaAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACgg                                                               <  1:946018/151‑1 (MQ=255)                                                                                        caAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc                                                            >  1:1297006/1‑150 (MQ=255)                                                                                         aaaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGccc                                                            >  1:897464/1‑149 (MQ=255)                                                                                           aCCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCa                                                        <  1:1242887/151‑1 (MQ=255)                                                                                            ccAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCAt                                                       <  1:1062585/151‑1 (MQ=255)                                                                                              aGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCt                                                     <  2:1218884/151‑1 (MQ=255)                                                                                              aGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCt                                                     <  2:1196397/151‑1 (MQ=255)                                                                                              aGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCt                                                     <  1:541475/151‑1 (MQ=255)                                                                                                  cgcgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCtt                                                 >  2:582009/1‑151 (MQ=255)                                                                                                  cgcgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCtt                                                 <  2:481418/151‑1 (MQ=255)                                                                                                   gcgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTc                                                <  2:1112914/151‑1 (MQ=255)                                                                                                        ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTc                                           >  2:1081020/1‑151 (MQ=255)                                                                                                        ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTc                                           >  1:158682/1‑151 (MQ=255)                                                                                                         gCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCa                                          <  2:456801/151‑1 (MQ=255)                                                                                                           ttGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATg                                        >  2:323732/1‑151 (MQ=255)                                                                                                               aCGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCt                                    <  1:1181066/151‑1 (MQ=255)                                                                                                                 gCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTAc                                  >  1:88139/1‑151 (MQ=255)                                                                                                                  cATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                                 >  2:91374/1‑151 (MQ=255)                                                                                                                  cATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                                 >  1:1088797/1‑151 (MQ=255)                                                                                                                  cATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                                 <  2:778768/151‑1 (MQ=255)                                                                                                                     aaaGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGc                               >  2:491251/1‑150 (MQ=255)                                                                                                                     aaaGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGc                               >  1:534862/1‑150 (MQ=255)                                                                                                                       aGCTGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGc                            <  1:244754/151‑1 (MQ=255)                                                                                                                          tGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTAc                         >  2:590206/1‑151 (MQ=255)                                                                                                                          tGTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTAc                         >  2:458294/1‑151 (MQ=255)                                                                                                                           gTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTa                          <  2:1001145/149‑1 (MQ=255)                                                                                                                           gTTGCTGTGCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACg                        <  1:907057/151‑1 (MQ=255)                                                                                                                                tgtgCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCtttt                   >  1:612073/1‑151 (MQ=255)                                                                                                                                 gtgCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCttt                    <  2:217573/149‑1 (MQ=255)                                                                                                                                 gtgCACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTc                  <  2:594389/151‑1 (MQ=255)                                                                                                                                    cACTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGt               >  2:577345/1‑151 (MQ=255)                                                                                                                                     aCTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGt               <  2:664781/150‑1 (MQ=255)                                                                                                                                     aCTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTg              <  1:127024/151‑1 (MQ=255)                                                                                                                                     aCTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTg              >  2:785070/1‑151 (MQ=255)                                                                                                                                     aCTTGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTg              >  2:175165/1‑151 (MQ=255)                                                                                                                                        tGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGaa            <  1:1229541/150‑1 (MQ=255)                                                                                                                                        tGTTTGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGa             <  2:922712/149‑1 (MQ=255)                                                                                                                                           ttGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCgg        >  1:452566/1‑151 (MQ=255)                                                                                                                                           ttGATGAAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCgg        >  1:96727/1‑151 (MQ=255)                                                                                                                                            tgatgaAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGc       >  2:156898/1‑151 (MQ=255)                                                                                                                                             gatgaAGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGc       >  2:53387/1‑150 (MQ=255)                                                                                                                                                 aaGTTCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGAt  >  1:600024/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TCAATCAGGAAGACTTCACGGGCACCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGAT  >  NC_000913/3400663‑3400957

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