Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,866,783	G→A	Q257* (CAG→TAG)	rpoS ←	RNA polymerase sigma factor RpoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,866,783	0	G	A	97.6%	119.4 / NA	44	Q257* (CAG→TAG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (15/26);  minor base C (1/0);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTA  >  NC_000913/2866539‑2867020                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCat                                                                                                                                                                                                                                         >  1:66731/1‑251 (MQ=255)         aCAGATGCGTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCAt                                                                                                                                                                                                                                 <  1:359034/251‑1 (MQ=255)                     ttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTAt                                                                                                                                                                                                                     >  1:149775/1‑251 (MQ=255)                          cgcgGAACAGCGCTTCGATCTTCAGCCCCTGTGTTTGCACGATTTCGCGCAAACGGGGCAGGCCTTCAACCTGAAGCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTTCATCTGCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGTCGTGCCAGCACTTCACGCCGTTTGGCGTCCAGTGTGGACAGCCATTTGACGAGGCTCTACTTCATATCGTCATCTTGCGGGGAATCTTcc                                                                                                                                                                                                                <  2:176122/251‑1 (MQ=255)                          cgcgGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTc                                                                                                                                                                                                                 <  1:223765/250‑1 (MQ=255)                                     gcTTCGATATTCAGCTCCTGCTGTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTtctc                                                                                                                                                                                                     <  2:359793/251‑1 (MQ=255)                                        tcgagagtCAGCCCCTGCGTTTGCAGGATTTCGCGCGAACGGCGCAGGGCTTCAACCTGAATCTGGCGAACACGTTCACGTGTGAGGCCTATTTCACGTCCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGGAAACCGAATGGACGTTCCATCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTCGACGTTGGGCTATTGCAGGTCGGGATCTTGCGTGGTATCTTCCGGACCTTTCTCtt                                                                                                                                                                                                   <  1:383078/243‑1 (MQ=255)                                         cGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCtttt                                                                                                                                                                                                 <  1:127665/251‑1 (MQ=255)                                            tatTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGAAGGACTTCAACCTGACTCTGGAGACCACGTTCACGGGTGAGGCCAATTACACCACCAACATCTAACACAGTTGCCGATTCGTACCCCAGCAAACCGAATCGACGTGACAGCACTTCACGCTGTTTGACGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTc                                                                                                                                                                                               >  1:250098/1‑250 (MQ=255)                                                cAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGATCACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACTCCAGCAAACCGAATCGACGTGCCAGCACTGCACTCGGTTTCGCGTTCAGCTCGATCAGCCATTTGACGATGCTCTACTTCATATCGTCGTCTGGCGTGGTATCGTCCGGTCCGTTCTCTGTTTCATCgg                                                                                                                                                                                          <  2:89591/251‑1 (MQ=255)                                                cAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCgg                                                                                                                                                                                          >  1:5227/1‑251 (MQ=255)                                                 aGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCACGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGATGCCGCTTCGTACCCCAGCAAACCGACTCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGCACGTTCTCTTTTTCATCGGc                                                                                                                                                                                         >  1:191257/1‑251 (MQ=255)                                                               cAGGATTTCGCGCAAACCCCGAAGGAATTCACCCTGAACATCGCGCACACCTTAACGGGTGCGCCCCATTTCACGACCTCCATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGCCGTGCCAGCCCTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGATCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTcccgc                                                                                                                                                                            >  2:54342/1‑247 (MQ=255)                                                                               cGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCACCCCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCCGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGTCCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGATTTTTCGGAATc                                                                                                                                                           >  2:140451/1‑251 (MQ=255)                                                                                      ggCCTTCAACCTGAATCTGGCGAACACGTTCACGAGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGAAACCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTAGGAATCACCAccc                                                                                                                                                    >  1:345540/1‑251 (MQ=255)                                                                                                  gAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTc                                                                                                                                         <  1:291672/250‑1 (MQ=255)                                                                                                        ggCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGa                                                                                                                                   <  1:161697/250‑1 (MQ=255)                                                                                                           gAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTa                                                                                                                               <  1:243397/251‑1 (MQ=255)                                                                                                                   tCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGGGTTGCCGCTTCGTACCCCAGCAAACCGACTCGGCGTGCCAACACTTCACGCTGTTTGGAGGTAACCTCGAACAGCCATTAGACGATGCTCTCCTTCATATAGTCATCTTGCGTGGTATCTTCCGGAACGTTCTCTTTTTCATCGGCCCGGATGTCCAGCAACCCTTTTTCGGAATCACCACCCAGCGGGGTAGCTACCGAGGGAATgcgatc                                                                                                                       >  2:119725/1‑248 (MQ=255)                                                                                                                     aCGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGGACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGGCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGt                                                                                                                     <  2:156053/251‑1 (MQ=255)                                                                                                                           gAGGCCAAGTTCACGGCCTACATCTTCCAGTGTTGCCTCTTCGTACCCCAGCAAACCGAAACGACGTGCCAGCACTTCCCTCTGTTTGGAGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTGCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGAc                                                                                                               <  2:234566/251‑1 (MQ=255)                                                                                                                                    ttCACGACCTACCTCTTACAGTGTTGCAGATTCATACACCAGCAAACCGAATAGACGCGACAGCAATTAACGCTGTTTGGCGTTCCGCTCGAACAGCCATTTGCCGATGCTCTACTTCATATCGTCATCTTGCGTGGAATCTACCCGACCGTTCTCTTTTTCATCGGCCAGGAAGTCCAACAACGCTTTTTCGGAATCACCACCCAGCGGGGGGTCTAACGAGGTACAGCTCTCGTTAATGCGATGCATAc                                                                                                      >  2:422245/1‑251 (MQ=255)                                                                                                                                    ttCACGACCTACAGCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACTCAGCGGGTTGTCTACCGAGGTAATGCGCTCGTTAAGACGCAGCATa                                                                                                       <  1:268913/250‑1 (MQ=255)                                                                                                                                                    tCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTTCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCTGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGTCTGACGTCATCaa                                                                                       <  2:248742/250‑1 (MQ=255)                                                                                                                                                               cGCTTCGTACCCCAGCAAACCGAAGCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCa                                                                           <  2:5227/251‑1 (MQ=255)                                                                                                                                                                 cTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCAATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGt                                                                         >  1:250597/1‑251 (MQ=255)                                                                                                                                                                          ccAGCAGACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGCTGTCCATCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGCCGAAGCATCCGGCTGACGTCATCAACTGGCTTATCCCGTTGCTCTGCg                                                                >  2:94993/1‑251 (MQ=255)                                                                                                                                                                            aGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAt                                                              <  2:140022/251‑1 (MQ=255)                                                                                                                                                                                  ccGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCTCCACCCAGCGGGGTGTCTACCGAGGTAATGTGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTc                                                        <  2:352562/251‑1 (MQ=255)                                                                                                                                                                                    gAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCg                                                      <  2:56198/251‑1 (MQ=255)                                                                                                                                                                                            gTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGc                                                                                                                                                 >  1:245729/1‑152 (MQ=255)                                                                                                                                                                                            gTGCCAGCACTTCACGCTGTTTGGCGTTCAGCCCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGt                                              <  1:133622/251‑1 (MQ=255)                                                                                                                                                                                                       tCACGCTGTTTGGAGTTCAGCTCGAACAGCCATTTTACGATTCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGTCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCa                                    <  2:358573/250‑1 (MQ=255)                                                                                                                                                                                                       tCACGCTAGTTGGCGTACAGCTCGAACAGCCATTTTACGATGCTCTACTACATATCGTAATCTTGCGTGGTATCTTCCGGACCGTTATCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGCAATCACCACCCAGCGGGCAGTCTACCGCGGTAATGAGCTCGTTAAGACGCAACATACGGCTTCCTTCATCAACTGGCTTATCCCGTTGCTCAGCGTTCTCTTCCGCACGTGGCTCATGCTCCac                                   >  2:355329/1‑250 (MQ=255)                                                                                                                                                                                                           gCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTa                               <  1:270834/251‑1 (MQ=255)                                                                                                                                                                                                              gTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGATCTACTTCATATCGTCATCTTGCGTGGAATCTTCCGGCCCGTTCTCTTTTTAATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATgg                            >  2:447029/1‑251 (MQ=255)                                                                                                                                                                                                                   gCGTTCAGCTCGAACAGCCATTTGACGATGATCTACTTCATACCGTCATCTTGCGAGGTATCTTCCGGAACGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGAAATCACCACCCAGCGGGGCGTCTACCGAGGTAATGCGCTCGTTACGACAAAGCATACGGCTGACGACATCAACTGGCTTATCCAGTTGCTCTGCGCTCCCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACa                        >  1:278803/1‑250 (MQ=255)                                                                                                                                                                                                                       tCAGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTc                    >  2:123751/1‑250 (MQ=255)                                                                                                                                                                                                                         aGCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTAGCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACg                 <  2:66731/251‑1 (MQ=255)                                                                                                                                                                                                                          gCTCGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAATTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGt                <  1:42304/251‑1 (MQ=255)                                                                                                                                                                                                                             cGAACAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCGCTTCCGCACGTGGTTCATGGTCCATCTTATGGGACAACTCACGTGCg             <  2:177869/251‑1 (MQ=255)                                                                                                                                                                                                                               aaCAGCCATTTGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCGACCGAGGTAATGCGCTCGTTAAGACGAAGCATTCGGGTGACGTCATCAACTGGCTTATCCAGTTGCTCGGCGATCTCTTCCGCACTTGGTTCTTGGTCCAGCTTATGGGACAACTCACGTGCGGt           <  1:204666/251‑1 (MQ=255)                                                                                                                                                                                                                                     cATTTGACGATGCTCTACTTCATATCGTCAGCTTGCTTTGTATCTTCCGGACCTTTCTCTTTTTCATCGGCCAGGATGTCCAGCACCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGTTAATGCGCTCGTTAAGACGAAGCATACGTCTGACGTCATCAACTGGCGTATCCGGTTGCTCTGCGATCTCTTCCGCACTTTGTTCATTGTCCAGCTTATGGGACAACTCACTTGCGGTTCGCAg     <  2:166981/251‑1 (MQ=255)                                                                                                                                                                                                                                         tGACGATGCTCTACTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTTCCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTa  <  1:180959/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              CCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTA  >  NC_000913/2866539‑2867020

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