Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	91,916	C→T	Y168Y (TAC→TAT)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	91,916	0	C	T	92.7%	106.1 / 2.0	41	Y168Y (TAC→TAT)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (16/22);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

GGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACA  >  NC_000913/91673‑92160                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ggTAAAAGCGATTGGGGCTGCCCCTAAAGAAGGGCAGGGCGTTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTACCGCGCGCAATATTCCGCTGGAGCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAATAACTGAAACTGCCGGGGTGTCATCTGCGTGAAGAATCTCGCCGTTATTAtccat                                                                                                                                                                                                                                               <  2:121601/251‑3 (MQ=255)                          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)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGG‑CTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAACTGAACA  >  NC_000913/91673‑92160

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: ‑