Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	92,096	A→T	A228A (GCA→GCT)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	92,096	0	A	T	100.0%	137.2 / NA	38	A228A (GCA→GCT)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (20/18);  total (20/18)

ACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCGGGCTCAACGGTTAA  >  NC_000913/91849‑92341                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      tcATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGATATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGCGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTgcg                                                                                                                                                                                                                                                    >  2:196344/2‑251 (MQ=255)                                    gggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCa                                                                                                                                                                                                                 <  1:245006/251‑1 (MQ=255)                                             cTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGt                                                                                                                                                                                                        >  2:253824/1‑251 (MQ=255)                                             cTGCGCGAAGACTGTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTggg                                                                                                                                                                                                        >  2:371472/1‑250 (MQ=255)                                               gCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGttt                                                                                                                                                                                                      <  1:169325/251‑1 (MQ=255)                                                  tGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGGGACTCCTCAACTCATTAGCTTTACTAACGTCGATAGTCAAGAGATTGAGGGCGTTGAGACGCGTTTCGCTAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGtttctc                                                                                                                                                                                                   >  1:371891/1‑248 (MQ=255)                                                      gaatctcGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATcgcg                                                                                                                                                                                                >  1:264946/4‑250 (MQ=255)                                                                              ggCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGCTAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACACCCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACATTGCGGTGGCCTTTa                                                                                                                                                                        >  2:236038/1‑250 (MQ=255)                                                                                cGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACa                                                                                                                                                                     <  2:320991/251‑1 (MQ=255)                                                                                          gCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGaa                                                                                                                                                             <  1:332205/249‑1 (MQ=255)                                                                                             cACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCt                                                                                                                                                          <  1:64042/249‑1 (MQ=255)                                                                                              aCCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGt                                                                                                                                                       >  1:308964/1‑251 (MQ=255)                                                                                              aCCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGt                                                                                                                                                       >  2:71669/1‑251 (MQ=255)                                                                                                          ttACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATCAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCtg                                                                                                                                           >  2:359621/1‑251 (MQ=255)                                                                                                                aCGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGAt                                                                                                                                     <  1:21371/251‑1 (MQ=255)                                                                                                                  gTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGt                                                                                                                                   >  1:128915/1‑251 (MQ=255)                                                                                                                         gTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACAcc                                                                                                                            <  1:107116/251‑1 (MQ=255)                                                                                                                                 aTTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGt                                                                                                                    <  2:128915/251‑1 (MQ=255)                                                                                                                                                      ttCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGcc                                                                                                <  1:198860/250‑1 (MQ=255)                                                                                                                                                            aaaTGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCAt                                                                                          >  2:34589/1‑250 (MQ=255)                                                                                                                                                            aaaTGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCAt                                                                                          >  1:80361/1‑250 (MQ=255)                                                                                                                                                                 gCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAAcc                                                                                    >  1:139983/1‑251 (MQ=255)                                                                                                                                                                    tACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAAccc                                                                                   <  1:345065/249‑1 (MQ=255)                                                                                                                                                                     aCCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTa                                                                                 <  2:308964/250‑1 (MQ=255)                                                                                                                                                                      ccGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAAc                                                                               <  1:268234/251‑1 (MQ=255)                                                                                                                                                                        gggCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACa                                                                              >  2:82846/1‑250 (MQ=255)                                                                                                                                                                                             gtgCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGaa                                                           >  1:43229/1‑248 (MQ=255)                                                                                                                                                                                                       aCCGCTATGGTCGCGTAATTGAAGATATTTCTTCTGCTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGt                                              <  1:47560/251‑1 (MQ=255)                                                                                                                                                                                                       aCCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGt                                              >  1:182824/1‑251 (MQ=255)                                                                                                                                                                                                             aTGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGTCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGGTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTTCTGGCGATGGCTAACAGCCCGTCATACAACCTTAACAATCTTAGCGGCACGCCGAAAGAGGCGATGCGTAACCGt                                        <  1:369354/251‑1 (MQ=255)                                                                                                                                                                                                              tGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACa                                                                              >  2:333463/1‑212 (MQ=255)                                                                                                                                                                                                              tGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACa                                                                              <  1:333463/212‑1 (MQ=255)                                                                                                                                                                                                                gTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTAcc                                     >  2:341976/1‑251 (MQ=255)                                                                                                                                                                                                                        aTTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGa                             <  1:137411/251‑1 (MQ=255)                                                                                                                                                                                                                                             gACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCTGGCTCa          <  2:71972/249‑1 (MQ=255)                                                                                                                                                                                                                                                    aGGCTGCGCACATCCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCGGGTAGCGCCTTGCTGGTAGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGc                                                <  1:174922/204‑1 (MQ=255)                                                                                                                                                                                                                                                    aGGCTGCGCACAACCTGGCGCTGAGAATTGCTGCACGCCTGCAGGCGCTGGTTTATCGCGCACTGCACAATGCGGTGGCCTTTAACAAGGCTGACTCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGc                                                >  2:174922/1‑204 (MQ=255)                                                                                                                                                                                                                                                    aGGCTGCGCAAAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCTGGCTCAACGGTTaa  >  2:142719/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGTGGCCTTTAACAAGGCTGAATCTGGTAGCGCCGTGCTGGTGGATGTCAACACCGGTGAAGTGCTGGCGATGGCTAACAGCCCGTCATACAACCCTAACAATCTGAGCGGCACGCCGAAAGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCGGGCTCAACGGTTAA  >  NC_000913/91849‑92341

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