Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	91,922	G→A	P170P (CCG→CCA)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	91,922	0	G	A	92.3%	102.4 / 0.6	40	P170P (CCG→CCA)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (16/20);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

AAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGT  >  NC_000913/91698‑92167                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      aaataagTGCATGAAGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTACCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAACGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCAt                                                                                                                                                                                                                             <  2:368017/247‑1 (MQ=255)    gaagTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCAttgg                                                                                                                                                                                                                          >  1:82771/1‑248 (MQ=255)        ggcactccGCTGGCGGTATCAGCGTCGGTGACAGCTGGAAGGCGCTGG‑CTCACTCTCTCAATATCCCGCTGGATCAGCTTTCCGCCCGCATTAACGCCAACCCGAAAGGTCGTTGTATTTATCTGGCGCGTCAGTTGACCCCTGACATTGCTGACTCCATCAAAAAAATGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCAGCCGGCGAAGTGACTGCTCACCTTATTGGCttt                                                                                                                                                                                                                      <  2:190604/244‑1 (MQ=255)                  gtgtgtttgcagagccGGTGACCGCTGGAAGGCG‑TGGTCTCACGCGCTCAATATTCCGCTGGTTCAGCTTTCATACCGTATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGGCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCCTCCGGCGAAGTGACTGCGCTCCTCATTGGCTTTACTAACGTc                                                                                                                                                                                                             <  2:221609/235‑1 (MQ=255)                                         gCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCAGTAACGCCAACCCGAAAGGGCGCTTTATGTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACGGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGtt                                                                                                                                                                                     <  2:461516/251‑1 (MQ=255)                                           gggAAGGCGCTTG‑CTAACGCGCTCAATATTCCGCTTGATCAGCTTTCAGCCCTCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTCACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCTTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGTTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTga                                                                                                                                                                                   <  1:420124/250‑1 (MQ=255)                                              aaGGCGCTGG‑CTAACGAGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGACGTTCTTATCCATCCGGCGAAGTGACAGCTCTCCTCATTGGCTTTACTAACGTCGTTAGTCCAGGGATTGAGGGCGTTgaga                                                                                                                                                                                 >  1:163090/1‑250 (MQ=255)                                                           aaCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGAt                                                                                                                                                                      <  1:337149/249‑1 (MQ=255)                                                                 cTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCt                                                                                                                                                              >  1:362644/1‑251 (MQ=255)                                                                  tCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCAGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATc                                                                                                                                                                                                                            >  2:66756/1‑188 (MQ=255)                                                                  tCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCAGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATc                                                                                                                                                                                                                            <  1:66756/188‑1 (MQ=255)                                                                   cAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGCCTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATACATGGCtt                                                                                                                                                             >  2:42944/1‑250 (MQ=255)                                                                    actatTCCGCTGGAACAGCTTTCAGCCCGCATTAAGGCCAAACCGTCAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCAGTTATTATCCATCCGGCGAAATGACTGCTCACCTCATTGGCTTACCTACCGCCGATAGTCAAGGGGTTGAGGGCGTTCagaagtatatcgataacaggctaa                                                                                                                                                            >  2:161372/3‑231 (MQ=255)                                                                          ccGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCa                                                                                                                                                     <  1:28649/251‑1 (MQ=255)                                                                                          cAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCAtt                                                                                                                                     >  1:199488/1‑251 (MQ=255)                                                                                          cAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTAATTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTAAATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTAGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGACTTACCGGGCAGCCGGGTGCGGGCAtt                                                                                                                                     >  2:314244/1‑251 (MQ=255)                                                                                           aGCCCGCATTAACACCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGTGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAATAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATtg                                                                                                                                    <  1:204362/251‑1 (MQ=255)                                                                                             cccGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGGCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATtgtg                                                                                                                                  <  1:215398/251‑1 (MQ=255)                                                                                                 cATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTa                                                                                                                              <  1:176689/251‑1 (MQ=255)                                                                                                     aaCGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTaa                                                                                                                             >  1:59638/1‑248 (MQ=255)                                                                                                                gAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTc                                                                                                               >  1:397548/1‑251 (MQ=255)                                                                                                                           ttatttatCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGa                                                                                                     <  1:42944/250‑1 (MQ=255)                                                                                                                                            cAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGCCTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTg                                                                                    >  2:272427/1‑250 (MQ=255)                                                                                                                                                tGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGc                                                                               >  1:426016/1‑251 (MQ=255)                                                                                                                                                acgctcttccgatctcGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGAc                                                                                                                         >  1:361199/16‑209 (MQ=255)                                                                                                                                                acgctcttccgatctcGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGAc                                                                                                                         <  2:361199/194‑1 (MQ=255)                                                                                                                                                   aCCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAg                                                                            <  1:195099/251‑1 (MQ=255)                                                                                                                                                        gACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGctg                                                                        <  1:51037/250‑3 (MQ=255)                                                                                                                                                              gCGGAATACATCAAACAACTGAAACCGCCGGGGTTTCATCTGCGTGAAGAATCTCTCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTTCTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACaa                                                                 <  2:121433/251‑1 (MQ=255)                                                                                                                                                                  aCTACATCAAAAAACGGAAACTGCCGGGGATTCATCTGCGGGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAAGTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTg                                                             <  2:397548/251‑1 (MQ=255)                                                                                                                                                                                   aaactgCCGGGGATTCATCTGCGTGAAGAATCGCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCCCAACCTGGCGCTGAGTATtgatga                                            <  2:355140/251‑1 (MQ=255)                                                                                                                                                                                     actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGCGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATtgatgacc                                          >  2:22779/1‑249 (MQ=255)                                                                                                                                                                                        gCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATACATGGCTTACCGGGCAACCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTAGCGTAATTGAAGATATTTCTTCAACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGcc                                       >  2:251837/1‑251 (MQ=255)                                                                                                                                                                                        gCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACg                                         <  2:603/249‑1 (MQ=255)                                                                                                                                                                                                 tCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGgcgc                              >  1:73602/1‑251 (MQ=255)                                                                                                                                                                                                       gCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGATTATTGATGAACGCCTGCAGGCGCTGGttt                        <  1:381720/251‑1 (MQ=255)                                                                                                                                                                                                         gTGAAGAATCTGGCCGTTATTATCCATCCGGCGAAGTGACTGCTAGCCGCAATGGCTTTACTAACGGCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGGCAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTAt                      <  2:425796/251‑1 (MQ=255)                                                                                                                                                                                                           gaagaaTCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATcg                    >  1:1269/1‑251 (MQ=255)                                                                                                                                                                                                                    cGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGaa           <  1:241318/251‑1 (MQ=255)                                                                                                                                                                                                                      ccGTTATTATCCATCCGGCGACGTGACTGCTCTCCTCAGTGGCTTTTCTAACGTCGATAGTCAAGTGAGTGAGGGCGTTGAGAAGAGGTTCGATAAATGGCTTACCGTGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTTACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCTCTGGTTGTTCGCTAACTGaaca         <  2:290343/251‑1 (MQ=255)                                                                                                                                                                                                                              atccatccGGCGAAGTGACTGCTCCCCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGAAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAATGCGGt  <  2:332031/250‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      AAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACAACGCGGT  >  NC_000913/91698‑92167

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