Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	881,531	A→T	I269F (ATC→TTC)	dacC →	D‑alanyl‑D‑alanine carboxypeptidase DacC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	881,531	0	A	T	89.5%	99.5 / 2.8	38	I269F (ATC→TTC)	dacC	D‑alanyl‑D‑alanine carboxypeptidase DacC
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (18/16);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

GCATTGCTGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTATCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTACCG  >  NC_000913/881288‑881753                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                gCATTGCTGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGAccgtttccgttt                                                                                                                                                                                                                          >  1:462931/1‑251 (MQ=255)        tGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTTCTTCGGCTACCCAGTGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGATAACCTGCCGTTTCCGTTTTAATg                                                                                                                                                                                                                     >  1:181264/1‑249 (MQ=255)                gtaTTGATCTACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGGTCACCTTTCACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGCATGTTGATGTCCTGAAGTGAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGACTCGGCTTCCCATGGTGATATGCTTTTAATCTCCGTAGTGATGGGGTCTAAAACCGACCGTTTCCTGTTTAATGAGTCTgagta                                                                                                                                                                                                           <  2:279447/249‑3 (MQ=255)                    tGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑att                                                                                                                                                                                                       >  1:195298/1‑250 (MQ=255)                     gATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAAAAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGATTTAATGAGTCTGAGAA‑attat                                                                                                                                                                                                     >  1:557523/1‑251 (MQ=255)                        ccACGATGTGCCGGAAGAGTACGTCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTATCCGTAACCGGCTGCTGGGGAGCAGCAATCTGAATGTTGATGGCATGAAGTCAGGAACCACGGCAGGCTCGGTATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGTCTAAAACCGCCCGTTTCCGTTTTAATTAGTCTGAGAA‑ATTATTGa                                                                                                                                                                                                  <  1:188999/251‑1 (MQ=255)                               gtgCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTggg                                                                                                                                                                                            <  1:478142/250‑1 (MQ=255)                                            aCGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGACAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCCTGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGAAATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTATCCGTTTTAATGCGTCTG‑GAATATTATTGACCTTGGTGTTCCCCTTCttt                                                                                                                                                                              >  2:185964/1‑251 (MQ=255)                                                    cATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTACCCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGt                                                                                                                                                                      <  2:445667/251‑1 (MQ=255)                                                          gagcaagagTTCACCTGCAACAAACTTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCACTCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATATGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCAGTAGTGTTGGGGGCGAAACCCGACCGTTTCCGTTTAAATGTGGCAGAGGA‑ATTATTGGCCTGGGGTTTCAGCTTCTTTGAAACCGTGActac                                                                                                                                                                >  2:24963/5‑248 (MQ=255)                                                             aaagagTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAAt                                                                                                                                                             >  1:315764/1‑251 (MQ=255)                                                              aagagTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAAtt                                                                                                                                                            <  1:230869/251‑1 (MQ=255)                                                                 agTTCCCCTTCAACAAAATGCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTaaa                                                                                                                                                         <  1:177917/251‑1 (MQ=255)                                                                   ttCACCTTCAACAAAATGCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGCATGTTGATGGCATGAAGACAGGAACAACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTCCCCCGGGCGATATGCGTTTAATCACCGAAGTGTTGGGGGCGACAACCGACCGTTTCCGTTTTAATGAGGCTGAGAA‑ATTATTGACCTGGGGTTTCAGCTTTTTTGAAACCGTGACGcacattcaaac                                                                                                                                                       >  2:372896/1‑241 (MQ=255)                                                                        cTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATg                                                                                                                                                  >  1:413608/1‑251 (MQ=255)                                                                                          gCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTATCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAgc                                                                                                                                <  1:185964/251‑1 (MQ=255)                                                                                           ccTATCCGTATCCGTCTGCTGTGGAGCCGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGTGATATAATCTGGTTGCGTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTGGGGGGCGAAAACCGTCCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTTACCTGGGGTTTCCTCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAgcg                                                                                                                               <  2:333597/251‑1 (MQ=255)                                                                                           ccTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGGTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGCAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAgcg                                                                                                                               >  2:309113/1‑251 (MQ=255)                                                                                                    aaccgtCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTTGGGGCGAAAACCGACCGTATCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAAcc                                                                                                                                                                        >  2:232834/1‑201 (MQ=255)                                                                                                    aaccgtCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGtt                                                                                                                      <  1:191252/251‑1 (MQ=255)                                                                                                    aaccgtCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTATCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAAcc                                                                                                                                                                        <  1:232834/201‑1 (MQ=255)                                                                                                            gctgTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGAAGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATa                                                                                                              <  1:149927/251‑1 (MQ=255)                                                                                                                 ggAGCAGCAATCTGAATGTTGTTGGCATTAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGTGGTTTCCGCTTCTTTGAAACCGTGACGTCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCGGGTTTGGTGATAAGAGc                                                                                                         <  2:315764/251‑1 (MQ=255)                                                                                                                             tGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAAtctc                                                                                             >  1:379178/1‑251 (MQ=255)                                                                                                                              gaaggTTGATGGCATGAAGACCGGAACCACGGCAGGAGCTGGATATAATCTGGTTGCGTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAATTGAAtct                                                                                              <  2:467384/245‑1 (MQ=255)                                                                                                                                      aTGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGCCGCCAATTAAACCTGCTGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGGGATAAGAGCGTAGGGAATCTCGGGGcgggg                                                                                    >  2:279856/1‑247 (MQ=255)                                                                                                                                        ggCATGAGGCCAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAACCCGACCGTTTCCGTTTTAATGCTTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAACCCTGATGCCCCCTTGGTGACTCAGCGCGTTTGGTTTGGTGATATGAGAGAAGGGAATATCGGGGCAGGTg                                                                                   >  2:384214/1‑250 (MQ=255)                                                                                                                                            tGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCg                                                                              >  1:258493/1‑251 (MQ=255)                                                                                                                                                     gAACCACGGCAGGCGCGGGATATAAGCTGGTTGCTTCGGCTACCCAGGGCTATATGCGTTTAATCTCCGTAGTGTGGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAATTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTg                                                                     <  1:177700/251‑1 (MQ=255)                                                                                                                                                                        atatAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGc                                                  >  1:4323/1‑251 (MQ=255)                                                                                                                                                                                     gCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGaa                                     <  1:301362/251‑1 (MQ=255)                                                                                                                                                                                                   gggCGATATGCGGTTAATCTCCGTCGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAATGACCATACCGCGTGGGCCGCTGAAAAACCTGAACGCGAGTTATACGt                       >  2:135058/1‑251 (MQ=255)                                                                                                                                                                                                       gATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAAc                   >  1:305860/1‑251 (MQ=255)                                                                                                                                                                                                                 tAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGc           >  1:439891/1‑249 (MQ=255)                                                                                                                                                                                                                      tcCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGATTCTGAGAA‑ATTAATGACCGGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACGCAGCGCGGCTGGTGTGGTGATAAGAGCGAAGTGaa                                                                                                 <  2:230020/158‑1 (MQ=255)                                                                                                                                                                                                                      tcCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGCTGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGaa                                                                                                 >  1:230020/1‑158 (MQ=255)                                                                                                                                                                                                                        cGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTaccg  <  2:413608/251‑1 (MQ=255)                                                                                                                                                                                                                                              gaccgTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTg                                       >  1:533951/1‑192 (MQ=255)                                                                                                                                                                                                                                              gaccgTTTCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTg                                       <  2:533951/192‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                GCATTGCTGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTATCCGTTTTAATGAGTCTGAGAA‑ATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTACCG  >  NC_000913/881288‑881753

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