Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,329,771	G→A	W270* (TGG→TAG)	dacB →	peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,329,771	0	G	A	97.6%	111.5 / NA	41	W270* (TGG→TAG)	dacB	peptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/20);  minor base C (0/1);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCGATATTGGAAACACCAT  >  NC_000913/3329537‑3330017                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       gggCATCTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACActgctg                                                                                                                                                                                                                                        <  2:90051/250‑1 (MQ=255)              ccccGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGGACACTGCTGCGCCCGACTCAgg                                                                                                                                                                                                                           >  2:28798/1‑251 (MQ=255)                ccGTTACGAGGTTCAGCCAGGTACGCGCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGtt                                                                                                                                                                                                                         <  2:25241/251‑1 (MQ=255)                 cgttactaggtgcAGCCAGTTACGTATCCTCCCCCGTGGTTCTGCCGTAGCGCAATGCGGCTAACTGGACGTTTGGCTAGTGGACCGGAACGTCTTTACTCTCACTGGAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTGTTGCCGTGCAGGATGGAGCCCGCTATGCCGGTGCAAGTCTGAAAGATGAGTTAAAACAGGCGAGTAGCACCTAGAGCGGAACACGGGTGCTCCACACTCaactta                                                                                                                                                                                                                        <  2:64621/239‑6 (MQ=255)                                      aCGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCCGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCCCCTAGAGGGGAACCCTGCTGCGCCAGACTCAGGTTAAAGAAACTGGGACGGTAggtg                                                                                                                                                                                                   >  2:393853/1‑248 (MQ=255)                                        gCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGACGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCACCTAGAGCGGAAAACTGCTGCGCCAGAATAAGGTTAAAGAACCTGGAAAGGTAGGTGcc                                                                                                                                                                                                 >  1:289838/1‑251 (MQ=255)                                                           tGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTTCTGCGCCAGACTCAGGTTAACGAAACTGGGCAGGTAATTGCCAGTAAACCGTCGGCCCCGc                                                                                                                                                                              >  1:399140/1‑251 (MQ=255)                                                                   cgcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAAAGGGAGTTGCCCGTAAACAGTCGGCCCCCCTGCACGac                                                                                                                                                                      >  2:327048/1‑250 (MQ=255)                                                                         aCTGCGAACTGGATGCGGTGGCAGGCGACCTGAACCGCTTTCCGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCTGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCtt                                                                                                                                                                <  2:97650/251‑1 (MQ=255)                                                                            gCGAACTGGAGGTGGGGCCAGGCGACCTGAACCTCTTTACGCTGCCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCtt                                                                                                                                                                <  2:11556/248‑1 (MQ=255)                                                                              gAACTGGCTGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCAACTAGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGGAGATGCCATAAAACAGTGGGCCCCGCCGCACCATCTGCTTa                                                                                                                                                               >  1:93578/1‑247 (MQ=255)                                                                                          gtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCAACTAGAGCGGGACACTGCTGCGCCCGAATCAGGTTAACGAACCTGGGACAGGAGGTGCCAGAAAACAGTCGGCCCCGCTGCACGAGCTGCTTAAGaatttggcggaaa                                                                                                                                                >  1:80984/1‑238 (MQ=255)                                                                                            gCCAGGCGATCTGAACCGCTTTACGCTGACGTGATGCCTGCCTCAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAg                                                                                                                                              <  1:181793/250‑1 (MQ=255)                                                                                                  cGACCTGAACCGCTTTTCGCTGCGGGGAGGCCTGCCAGAACGTTCGGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCGGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGaca                                                                                                                                       <  2:93578/251‑1 (MQ=255)                                                                                                     ccTGACCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaaca                                                                                                                                    <  2:232671/251‑1 (MQ=255)                                                                                                                    gCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAAAACTGCTGCGCCCGACTCAGGGTAACGAACCTTGAACGGTCGATGCCAGTAAACAGTCGGCCCCCCTGCACGATCTGCTTAAGATTATGCTGAAAAAATCGGACAAAATGATCGCAGCTACgg                                                                                                                     >  2:95222/1‑251 (MQ=255)                                                                                                                      tGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCCCCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACg                                                                                                                                                                        <  1:122683/198‑1 (MQ=255)                                                                                                                       gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGtt                                                                                                                   <  1:271215/250‑1 (MQ=255)                                                                                                                          gggATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGGAAACAGTCGGCCCCGCTGCACGAACTGCTTAAGATTATGCTGAAAAAGTCGGACAACCTGATCGCCGATACGGTTTTcc                                                                                                               >  1:229694/1‑251 (MQ=255)                                                                                                                                cctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGa                                                                                                         <  1:95222/251‑1 (MQ=255)                                                                                                                                   gccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCACCTAGAGCGGAACACAGCTGCGCCAGACTAAGGTTAAAGAACCTGGAACGGAAGTTGCCAGTAAACAGTAGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCCCCGATACGGTTTTCCGCAAGAt                                                                                                        >  1:401668/1‑249 (MQ=255)                                                                                                                                          cgtgcTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATATACCATgc                                                                                               <  1:113297/247‑1 (MQ=255)                                                                                                                                                             tGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCCGTAAACAGTCGGCCCCCCTGCACGGTCTGCTTAAGATTATGCTGAAAAAGGCGGACAACATGATCGCCGATACGGGTTTCCGCATTATAGGCCCTGAGCGCTTTAATGTGCCTgg                                                                             >  2:82225/1‑250 (MQ=255)                                                                                                                                                                     gCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCg                                                                    <  2:310707/251‑1 (MQ=255)                                                                                                                                                                        gTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCCCCAGACTCAGGGTAACGAACCTGGGACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGGACGAACTGCTTTAGAATATGCTGACAAAGTCGGACAAAATGATCGCCCATACGGGTTTCCGCATGGAAGGCCATGCCAGCTTCCATGGTCCTGGGACATGGCgggg                                                                 >  2:151104/1‑250 (MQ=255)                                                                                                                                                                            gggagggAGCCAGATATGCCGGTGTAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCggg                                                             <  2:93914/246‑1 (MQ=255)                                                                                                                                                                            aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCggg                                                             <  1:28798/251‑1 (MQ=255)                                                                                                                                                                              gATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTc                                                           <  1:299237/251‑1 (MQ=255)                                                                                                                                                                                     ccAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCCAGAGCGGAACACTGCCTCGCCCGACTCAAGGTAACGCACCTGGAACGGGAGGTGCCAGTAAACAGGCGGCCCCGCTGCCCGGTCTGCTTTAGATTATGCTGAAAAAAGCGGACAACATACACTCCCATAAGGGTTTTCGCATGATAGGCCCTGGGCGCTTCAATGTGCCTGGCACCAGGCGGGCAGGGGCggcggcc                                                    >  2:191075/1‑246 (MQ=255)                                                                                                                                                                                       aGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGt                                                  <  2:66963/251‑1 (MQ=255)                                                                                                                                                                                        gCTATGCCGGTGGAAGTCGGAAAGATGAGTGAAAACAGGCGGTTATCACCAAGAGTGGAACACTCCTGCGCCGGACTCAGGTTAACGAACCTGGAAAGGTAGTGTCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGAGTATGCGTAAAAAGGCGGACCGCATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGCTGCAATGTGTCTGGGTCATGGCGGGCCGGGTCGGGTGCCGTg                                                 <  2:44736/251‑1 (MQ=255)                                                                                                                                                                                          tATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAAAACTGCTGCGCCAGAATCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTAGGCCCCGCTGCCCGATCTGCTTAAGATTATGCCGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCCTGAAAGGCCATGCGCGCTTCAATGTGCCTGGAACCTGGCGGGACggggggggacccgcggg                                               >  2:371383/1‑237 (MQ=255)                                                                                                                                                                                           aTGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGAATCAGGTTAACGAACCTGGGACGGTAGTTGCCAGTAAACAGTCGGCCCCGGTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAAAATGATCGCCGATACGGGTTTCCGCATGATAAGACATGCGCGCTTCAAAGGGCCTGGAACATGGCGGGGCGGGGCGGACGCCGGGCgg                                              >  1:64210/1‑250 (MQ=255)                                                                                                                                                                                               cGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGGACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCCATGTGCCTGGAACATGGCGGGCCGGGTCGGGCGCCGGGCGGCAgg                                          >  1:271607/1‑250 (MQ=255)                                                                                                                                                                                                        tCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTAGAGCGGAAAACTGCTGCGCCAGACTCAGGTTAAAGAACCTGGAAAGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGGCGTGCGGCCGATCCTGCGcc                                 >  1:107224/1‑251 (MQ=255)                                                                                                                                                                                                                       aaaaCAGGCGGGTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGTTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCg                  <  2:229694/251‑1 (MQ=255)                                                                                                                                                                                                                            aGGCGGGTATCCCCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGTCGTGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCGATAtt             <  1:373713/251‑1 (MQ=255)                                                                                                                                                                                                                             ggCGGGTATCACCGCGACCGGAACACTGCTGCGCCAGGCTCAGGTTAACGACCCGGGAGCGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAAGTTGTGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATTATAGTCCATGCGCGCTGGAATGTGGCTGGAACATGGCGGGCCGGGTCGGACGCCGGGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCGATATTg            <  2:183850/251‑1 (MQ=255)                                                                                                                                                                                                                                gggTATCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCGATATTGGa          >  1:86696/1‑250 (MQ=255)                                                                                                                                                                                                                                    aTCACCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGACGGGACGGGGCGGACGCCGTGCGCCAGATCCTTCGCCAGCAAACCGGGGTCGAGATTGGaacacc     >  2:392773/1‑247 (MQ=255)                                                                                                                                                                                                                                       aCCTAGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGCTCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGCCGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGGGGTCGATATTGGAAacacaat  >  2:234885/1‑248 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       TGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGCAAGCCGGTGTCGATATTGGAAACACCAT  >  NC_000913/3329537‑3330017

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