Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,236,742	C→T	P161S (CCG→TCG)	yjbF →	lipoprotein YjbF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,236,742	0	C	T	97.8%	160.7 / NA	46	P161S (CCG→TCG)	yjbF	lipoprotein YjbF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/17);  minor base G (0/1);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTGGGAGCGTCATCGGTCTTCGCGGCAGGAACCGTCAAGGTGTTCAGCAAT  >  NC_000913/4236506‑4236983                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  cTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAACGTCGGCTGCGATGAAACATCGGTTCGCGTGCTg                                                                                                                                                                                                                                     >  2:14734/1‑251 (MQ=255)      ctctttACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATATCCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCGCGATGGACTGGACCGAGTACCAGCATGTACGCTACGCCACCGCCCGCTCAGTCTTCAAATGGGATGGCACCGATACCATCAAAATCGGCAGCGATGAAACCTCGGTTCGCGTGCTggtc                                                                                                                                                                                                                                  <  2:210572/245‑3 (MQ=255)      ctcGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAACGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACg                                                                                                                                                                                                                                 >  2:11344/1‑250 (MQ=255)        cGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCGCGGTTCGCGTGCTGGACgaag                                                                                                                                                                                                                              <  2:186930/251‑1 (MQ=255)           tACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACgaagaa                                                                                                                                                                                                                            <  1:239951/250‑1 (MQ=255)            aCCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGaagaag                                                                                                                                                                                                                           <  1:47857/250‑1 (MQ=255)                aaCATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTcc                                                                                                                                                                                                                      <  1:333532/251‑1 (MQ=255)                 agaTGGCCGTCTGGTGAAGACTCAGCTCGGCGGCGACAACCTGATTGCAGAGAATAACCTTGCTGCCGACACGCTCATTGAACCCGAACAAATCGTTGACGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGATACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAACGTCGGCAGCGATGAAACCTCGGTTCGCGTGATGGACGCATAAGTCTcc                                                                                                                                                                                                                      >  2:48089/3‑250 (MQ=255)                     ggCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCg                                                                                                                                                                                                                  >  1:20604/1‑250 (MQ=255)                     ggCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCg                                                                                                                                                                                                                  <  2:20604/250‑1 (MQ=255)                             ggTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGgc                                                                                                                                                                                                           <  2:155803/249‑1 (MQ=255)                                 aaGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTg                                                                                                                                                                                                     <  2:214075/251‑1 (MQ=255)                                    aCTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCa                                                                                                                                                                                                  <  1:325350/251‑1 (MQ=255)                                         gCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAAcc                                                                                                                                                                                             <  2:54051/251‑1 (MQ=255)                                                                          ccTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGcc                                                                                                                                                            >  1:101926/1‑251 (MQ=255)                                                                              gctgccgtaccacTGATTAATCCCTCCCAAATCGTTGATGGCGCAAGCTGGACGCGCACGAGGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCCCGGCACGCTCAGTCGTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAGCCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGTCGAGGCGCGCTGGCATAGCCGCTATTGGATCGACAGGGAAGGGCAAATTCTCCAGt                                                                                                                                                         <  2:407379/238‑1 (MQ=255)                                                                                tgcCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCgg                                                                                                                                                      >  1:10548/1‑251 (MQ=255)                                                                                   cGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAAc                                                                                                                                                   >  2:171024/1‑251 (MQ=255)                                                                                    gACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAAc                                                                                                                                                   >  2:108775/1‑250 (MQ=255)                                                                                              ttAAACACGCACACATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTAAAAGTCGGCAGCGATGAAACCTCGGATCGCGTGCTGGACGAAGAAGTCTCCACCGACAAGGCGCGCTGGCATAACCGCTATTGGATCGCCAGCGAAGGGCAAATTCGCCCGTCGGAACAGTATCTNGgc                                                                                                                                        >  2:378020/1‑251 (MQ=255)                                                                                               tAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTAt                                                                                                                                                                                        >  1:84866/1‑202 (MQ=255)                                                                                               tAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTAt                                                                                                                                                                                        <  2:84866/202‑1 (MQ=255)                                                                                                      gCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGatt                                                                                                                                 >  1:407259/1‑250 (MQ=255)                                                                                                        acaAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGCACAGTATCTCGGCGCGGattatt                                                                                                                              >  2:53231/1‑251 (MQ=255)                                                                                                               gTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCGAAGTCGGCAGCGATGAAACCTCGGTTCGCGTACTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGt                                                                                                                       >  1:67450/1‑251 (MQ=255)                                                                                                                          aaGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGACGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCActc                                                                                                            >  2:192355/1‑251 (MQ=255)                                                                                                                          aaGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCActc                                                                                                            >  1:415073/1‑251 (MQ=255)                                                                                                                            gCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACtctc                                                                                                          >  1:226184/1‑251 (MQ=255)                                                                                                                               ggACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGGCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGACAACCACTCtcatc                                                                                                       <  2:131883/251‑1 (MQ=255)                                                                                                                                gACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCtcatc                                                                                                       >  2:100756/1‑250 (MQ=255)                                                                                                                                 aCGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCaa                                                                                                     <  1:392156/251‑1 (MQ=255)                                                                                                                                     gcACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCACATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAggcg                                                                                                 >  2:330435/1‑251 (MQ=255)                                                                                                                                           tGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCATCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAAAAGTCTCCACCGACCAGGCGCGCTGGCATCACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAATCGGTACAGTATCTCGGCGCGGTTTATTTTCCGTTGAAAACCACTCTCATCAAGTCGGCaaa                                                                                            >  1:131597/1‑250 (MQ=255)                                                                                                                                                       aGTACAAGCAGGTACGCTAAGCCTCCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAACGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGCAGGGCAAATTCGCCAGTCGGCACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTaaac                                                                               >  2:227508/1‑251 (MQ=255)                                                                                                                                                             agcagGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACt                                                                          >  1:362417/1‑250 (MQ=255)                                                                                                                                                                                   aCGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAgtgt                                                   >  1:115119/1‑251 (MQ=255)                                                                                                                                                                                    cGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGtgtg                                                  >  1:85698/1‑251 (MQ=255)                                                                                                                                                                                        cAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTggg                                                <  1:332902/249‑1 (MQ=255)                                                                                                                                                                                        cAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTGGGa                                               <  1:196467/250‑1 (MQ=255)                                                                                                                                                                                              tCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTGGGAGCGTc                                          <  1:335388/249‑1 (MQ=255)                                                                                                                                                                                                    gggATGGCACCGATACCGCCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTGGGAGCGTCATCGGTCt                                  >  2:260440/1‑251 (MQ=255)                                                                                                                                                                                                      gATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTg                                                              <  1:70106/221‑1 (MQ=255)                                                                                                                                                                                                      gATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTg                                                              >  2:70106/1‑221 (MQ=255)                                                                                                                                                                                                                                 aGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTGGGAGCGTCATCGGTCTTCGCGGCAGGAACCGTCAAGGTGTTCAGc     >  1:297697/1‑251 (MQ=255)                                                                                                                                                                                                                                 aGCGATGAAACCTCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCCCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGAGTGGGAGCGTCATCGGTCTTCGCGGAAGGAACCGTCAAGGGGTTCAGc     >  2:94107/1‑251 (MQ=255)                                                                                                                                                                                                                                    gATGAAACCTCGGTTCGCGTGCTCGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTGGGAGCGTCATCGGTCTTCGCGGCAGGAACCGTCAAGGTGTTCAGCAAt  >  1:231114/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  CTACTCTCGTTACCCAACATGGCCGTCTGGTGAAGACTCTGCTCGGCGGCGACAACCTGATTGAAGTGAATAACCTTGCTGCCGACCCGCTGATTAAACCCGCACAAATCGTTGATGGCGCAAGCTGGACGCGCACGATGGGCTGGACCGAGTACCAGCAGGTACGCTACGCCACCGCACGCTCAGTCTTCAAATGGGATGGCACCGATACCGTCAAAGTCGGCAGCGATGAAACCCCGGTTCGCGTGCTGGACGAAGAAGTCTCCACCGACCAGGCGCGCTGGCATAACCGCTATTGGATCGACAGCGAAGGGCAAATTCGCCAGTCGGAACAGTATCTCGGCGCGGATTATTTTCCGGTGAAAACCACTCTCATCAAGGCGGCAAAACAATGATTAAACAAACTATTGTCGCGTTGCTTTTGAGTGTGGGAGCGTCATCGGTCTTCGCGGCAGGAACCGTCAAGGTGTTCAGCAAT  >  NC_000913/4236506‑4236983

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