Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,528,272	C→T	P602S (CCG→TCG)	igaA →	inner membrane protein IgaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,528,272	0	C	T	100.0%	169.0 / NA	51	P602S (CCG→TCG)	igaA	inner membrane protein IgaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/23);  total (28/23)

GGGACGCGCTGGTAAAACGCGCCAACGCCGGGAAGCTCGATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGC  >  NC_000913/3528027‑3528512                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 gggACGCGCTGGTAAAACGCGCCAACGCCGGGAAGCTCGATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACtcgccg                                                                                                                                                                                                                                             >  1:321554/1‑251 (MQ=255)        gcTGGTAAAACGCGCCAACGCCGGGCAGCTCGATGGCGTGAATGTGTTAGTACCCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCATCCTCCACCGCGTCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCGCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCATCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAAc                                                                                                                                                                                                                                      <  2:328079/251‑1 (MQ=255)            gTAAAACGCGCCAACGCCGGGAAGCTCGATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGAATCACTTTACGACTACTCGCCGCAAGAACAGTg                                                                                                                                                                                                                                  >  2:275876/1‑251 (MQ=255)                             cGGGAAGCTCGCTGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACATGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCCGAAAc                                                                                                                                                                                                                 >  2:43592/1‑251 (MQ=255)                             cGGGAAGCTCGATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAAc                                                                                                                                                                                                                 >  2:306280/1‑251 (MQ=255)                                    cTCGATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGcac                                                                                                                                                                                                           >  1:124796/1‑250 (MQ=255)                                    cTCGATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGcac                                                                                                                                                                                                           >  2:381805/1‑250 (MQ=255)                                       gATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGCTTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAAt                                                                                                                                                                                                       >  2:221004/1‑251 (MQ=255)                                               gAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGACCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGCTTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCACGAACAGTGGAACGCTTTCCAGAAACTGGCACACATGCTgcacc                                                                                                                                                                                               >  2:340109/1‑247 (MQ=255)                                                   gtgtTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATa                                                                                                                                                                                            <  1:398885/250‑1 (MQ=255)                                                       tattaCGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTGCTCGCCGCATGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGcc                                                                                                                                                                                        <  1:166607/250‑1 (MQ=255)                                                        attaCGGCCGGTGATTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCCCCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATCCGCCGt                                                                                                                                                                                      <  2:135164/251‑1 (MQ=255)                                                            cGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTa                                                                                                                                                                                   >  1:273571/1‑250 (MQ=255)                                                             gCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAAc                                                                                                                                                                                 <  1:384967/251‑1 (MQ=255)                                                             gCCCGGTGAGTGCGGAATCGCGGGATAACCTGCTGGCAACCTCCACCGCGCCGTTTATCCCGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGACCAGCCCTGGCCTTCGGCATCCCTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAAc                                                                                                                                                                                 <  2:294852/251‑1 (MQ=255)                                                              cccGGTGAGACCGGAATGGCTGGATAACCTGGTGGCAACCTCCACCGCGACGTTTATAACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCTGGCACCCGGCGGATTCCTGATTGGCAGCGACGAAGGCAGCGATTTTGTTGTGCAGCCCCGGCCTTCGGCATCACTTTACGAGTACTCGCAGCAAGAACAGAGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCAAACACCGGTTAAcc                                                                                                                                                                                >  2:229459/1‑250 (MQ=255)                                                                cGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCAACGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGcc                                                                                                                                                                              >  2:248831/1‑251 (MQ=255)                                                                        gCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTa                                                                                                                                                                       <  2:124796/250‑1 (MQ=255)                                                                            aaTCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTc                                                                                                                                                                  <  1:106419/251‑1 (MQ=255)                                                                                             gtggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGa                                                                                                                                                 <  1:22051/251‑1 (MQ=255)                                                                                             gtggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGa                                                                                                                                                 <  2:142128/251‑1 (MQ=255)                                                                                                ccAACCTCCACCTCGCCGTGTATCACGCATGCCACGGCCCGCGCGGCACAATCACTAAACACTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCTAGGGTACCGTCACAAAAATCTTCCCTGACGc                                                                                                                                              <  2:340538/250‑1 (MQ=255)                                                                                                           cgcgCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACgc                                                                                                                                   <  2:106506/251‑1 (MQ=255)                                                                                                                cGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCat                                                                                                                              <  2:256418/251‑1 (MQ=255)                                                                                                                         cGCATGAAACGGCGCGCGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCtt                                                                                                                     >  1:333155/1‑251 (MQ=255)                                                                                                                          gCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAAAGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTc                                                                                                                    >  1:256217/1‑251 (MQ=255)                                                                                                                           cATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCa                                                                                                                   <  1:351359/251‑1 (MQ=255)                                                                                                                                       cgAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGCTTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCgg                                                                                                        >  2:383884/1‑250 (MQ=255)                                                                                                                                        gagcgacaaaatgcctaaaccggcaggCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGACAGCGATATTGTTGATCAACCCTGCCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGCGGAACGCTTTCCAGCAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTAAAAAAACTCTTCACTGAAGCCAATGGAACGAAGCATATTGGCCTTCATCCGATCCCGGAt                                                                                                      >  2:396461/26‑251 (MQ=255)                                                                                                                                           cGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGt                                                                                                   >  2:99344/1‑251 (MQ=255)                                                                                                                                                 aaTCACTAAACAGTCCGGCACCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGc                                                                                             >  1:342525/1‑251 (MQ=255)                                                                                                                                                     aCTAAACAGTCCGGCCCCCGGCGGATTCATGCTTGTCAGAGAAGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCAATTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGGATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATAGTGGCCTTCATCCGATCCCGGATCGTTCCGGCCtgt                                                                                         >  1:214150/1‑251 (MQ=255)                                                                                                                                                     aCTAAACAGTCCGGCCCCCGGCGGAGTCCTGATTGTCAGCGACGAAGGCAGCGATGTTGTTGATCAGCCCTGGCCTTCGGGATCACTTTACGACTACTCGCCGCAAGAACAGTGGAAAGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTGAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCATCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCtgt                                                                                         <  1:340109/251‑1 (MQ=255)                                                                                                                                                             gTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTAt                                                                                 >  1:207777/1‑251 (MQ=255)                                                                                                                                                                         gcggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACcaca                                                                     <  1:17986/251‑1 (MQ=255)                                                                                                                                                                                   gATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTAc                                                           <  1:221004/251‑1 (MQ=255)                                                                                                                                                                                    aTTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACt                                                          >  2:374189/1‑251 (MQ=255)                                                                                                                                                                                      tGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTg                                                         <  2:11158/250‑1 (MQ=255)                                                                                                                                                                                                     cAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAgcg                                         >  2:18787/1‑251 (MQ=255)                                                                                                                                                                                                               gttgATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACa                                <  1:381805/250‑1 (MQ=255)                                                                                                                                                                                                                 tgATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAAt                              <  2:200498/250‑1 (MQ=255)                                                                                                                                                                                                                  gATCAGCCCTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATgg                            >  1:396298/1‑251 (MQ=255)                                                                                                                                                                                                                          caggCCTTCGGCATCAATTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAAAGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGACTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGAAGATGCTGGGTAGCGCCATTTACCATGGCGTCCAgg                    >  2:112860/3‑251 (MQ=255)                                                                                                                                                                                                                          cTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAgg                    >  1:252831/1‑251 (MQ=255)                                                                                                                                                                                                                          cTGGCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAgg                    >  2:291766/1‑251 (MQ=255)                                                                                                                                                                                                                            ggCCTTCGGAATCCCTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTATCGCCATTTACAATGGCGTACAGGcc                  <  2:13973/251‑1 (MQ=255)                                                                                                                                                                                                                             gCCTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCt                 >  1:63491/1‑251 (MQ=255)                                                                                                                                                                                                                               cTTCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTgg               >  1:279704/1‑251 (MQ=255)                                                                                                                                                                                                                                 tCGGCATCACTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTgg               <  2:203238/249‑1 (MQ=255)                                                                                                                                                                                                                                   ggCATCCCTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGACCCCGGGTCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGcgtc           <  2:203807/251‑1 (MQ=255)                                                                                                                                                                                                                                          cTTTACGACTACTCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTAcc     >  2:18673/1‑250 (MQ=255)                                                                                                                                                                                                                                            ttACGACTACTCGCAGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGc  >  1:65185/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 GGGACGCGCTGGTAAAACGCGCCAACGCCGGGAAGCTCGATGGCGTGAATGTGTTATTACGCCCGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCATATTGGCCTTCATCCGATCCCGGATCGTTCCGGCCTGTGGCGCTATCTCAGCACCACATTGCTGCTACTGACGATGCTGGGTAGCGCCATTTACAATGGCGTACAGGCCTGGCGTCGTTACCAGC  >  NC_000913/3528027‑3528512

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