Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,114,152	T→C	P67P (CCA→CCG)	pppA ←	prepilin peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,114,152	0	T	C	83.6%	140.1 / 18.8	67	P67P (CCA→CCG)	pppA	prepilin peptidase
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base C (26/30);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

CCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGGTTGTTTGTTGTAACTTGCTGGATGCGGCGTAAAA  >  NC_000913/3113909‑3114395                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    cctttGGGGATCGAGGCCGATACCCCTCGCGGCAATCACCCAGGCGGAAAATAACACCACCGCCAGCCCCCATCCACTTTCGGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGt                                                                                                                                                                                                                                              <  2:485548/248‑1 (MQ=255)      ggTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGa                                                                                                                                                                                                                                          >  1:430838/1‑250 (MQ=255)       gTGATCGAGGTCGATGACACCCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAAc                                                                                                                                                                                                                                        <  1:25779/251‑1 (MQ=255)         gATCGAGGTCGATGACACTCGCGGCATTCAGCCAGGCGGATAAGATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCATACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGGCACGTATCCGGATGGTCTGCTGACAAGGTGGGCTGTGCGAAcgc                                                                                                                                                                                                                                      <  2:586113/251‑1 (MQ=255)         gATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATAGCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAAcg                                                                                                                                                                                                                                       <  1:193778/250‑1 (MQ=255)         gATCGAGGCCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTGTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAAcg                                                                                                                                                                                                                                       <  1:612143/250‑1 (MQ=255)             gAGGTCGATGACCCGGGCGGCAACCAGCCAGGCGTATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACGCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCGCCAGCGTATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTa                                                                                                                                                                                                                                  <  1:204916/251‑1 (MQ=255)              aGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGt                                                                                                                                                                                                                                   >  2:230298/1‑249 (MQ=255)               ggTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGAAgc                                                                                                                                                                                                                                >  1:60806/1‑251 (MQ=255)                gTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAgcg                                                                                                                                                                                                                               <  1:463927/251‑1 (MQ=255)                    atgccacTCGCGGCAACCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCGGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACGCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAg                                                                                                                                                                                                                           <  2:405110/247‑1 (MQ=255)                      gACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGc                                                                                                                                                                                                                         <  1:203523/251‑1 (MQ=255)                      gACACTCGCGGCAATCAGCCAGGCGGATAATATCAACACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGc                                                                                                                                                                                                                         <  2:117950/251‑1 (MQ=255)                       acacTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCt                                                                                                                                                                                                                        >  1:255411/1‑251 (MQ=255)                                           ggCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGAAGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCg                                                                                                                                                                                                    <  1:466201/251‑1 (MQ=255)                                                 taataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACAttt                                                                                                                                                                                              >  1:500252/1‑251 (MQ=255)                                                    tatCATAACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAAATCCACCAGCGGATAACGATTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCGGGAGCACCGCGGAATATTGTACCGAATCCGGATGGTCTGATGACACTGCGGGCAGTGCGAACGCGGGAGAGAAAGGCTTATTTTTGACTGCGCACTCGACATTTCAc                                                                                                                                                                                           >  2:577141/1‑251 (MQ=255)                                                         tcaCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTg                                                                                                                                                                                       >  1:186853/1‑250 (MQ=255)                                                          caCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGa                                                                                                                                                                                      >  1:24178/1‑250 (MQ=255)                                                              gCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTc                                                                                                                                                                                 >  1:191024/1‑251 (MQ=255)                                                               ccAGCCCCCATCCACTGTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTcc                                                                                                                                                                                <  1:426268/251‑1 (MQ=255)                                                                      ccATCCCCTTCCTGGCCAGACCAGGCTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGATTACAATGCGGGCAGTGCGAACGCGATAGCGCAAGGCTTATTTTTGTCTGCGCACTCGACATTTCACCGTGAAACTCCGCCAttt                                                                                                                                                                         <  2:158026/251‑1 (MQ=255)                                                                                   ggCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCgcat                                                                                                                                                               <  2:18028/248‑2 (MQ=255)                                                                                           cagacTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCACCGCATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGCGCACAGCGGCATCTTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCCTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGt                                                                                                                                                    >  1:31012/1‑251 (MQ=255)                                                                                               cTCGCCAGCAACAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCATGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAAcgc                                                                                                                                                                                                                                      >  2:319281/1‑165 (MQ=255)                                                                                               cTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAAcgc                                                                                                                                                                                                                                      <  1:319281/165‑1 (MQ=255)                                                                                                 cGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGcc                                                                                                                                              >  1:75079/1‑251 (MQ=255)                                                                                                 cGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGAGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGcc                                                                                                                                              >  2:247702/1‑251 (MQ=255)                                                                                                   caaGCAAAAAAGAAAGGGCTGTCAATAACTCCACCAGCGGATAACGATGGCTGAGTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTGGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTGGTTGGCGCAGCAGGATTGGGTAACGCCaa                                                                                                                                            <  2:304132/249‑1 (MQ=255)                                                                                                        aaaaaaGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATcac                                                                                                                                       >  1:228882/1‑251 (MQ=255)                                                                                                         aaaaaGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGGGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATcacc                                                                                                                                      >  2:158798/1‑251 (MQ=255)                                                                                                         aaaaaGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATcacc                                                                                                                                      >  1:152849/1‑251 (MQ=255)                                                                                                                aaGTGCTGTCAATAACTCCGCCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAGTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCTCATTCa                                                                                                                               <  2:192365/251‑1 (MQ=255)                                                                                                                   tGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCaaa                                                                                                                             >  1:314043/1‑250 (MQ=255)                                                                                                                          aaTAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGcc                                                                                                                     >  1:345503/1‑251 (MQ=255)                                                                                                                                tCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCgatgat                                                                                                               >  1:544728/1‑251 (MQ=255)                                                                                                                                 ccaccaGCGGAGAACTCTGGCGGATGTTCGCCTGACAGGCGCGGCGGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATAGTGTGACGGATCCGGATGGTCTGCGGACAATGCGGGCAGTGCGAACGCGTTAGCGCAAGGCTTAGTGTTGACTGCGCACGCGACATGTCACCGTGAGACTCCGCCATTGGTTGGCGCAGCATGATTTTTTAAGGCGAAAGCACCACAGCCAAAAAACTGCCGGTGATc                                                                                                              <  2:31012/251‑1 (MQ=255)                                                                                                                                    ccaGCGGATAACGCTTTCTGATCGTCGCCTGACCGTCGCGGTAGCACCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATCGTCACGTATCCGGATGGTCTTCCGACAATGCTGGCAGTGCGAACGCGGTAGCGCAAGTCTTATTTTTGTCGGCGCACTCGAAATTACACCGTGAAACTCCTCCATTTGTTGGCGCAGCTTGATTGGGTAACGCCGTATCACCCCATTCAAAAACCTTCCTATGTTCAAt                                                                                                           <  2:367469/251‑1 (MQ=255)                                                                                                                                     caGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAAtc                                                                                                          >  1:135572/1‑251 (MQ=255)                                                                                                                                        cGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGATCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACCTTCAATAAACTGCCGATGCTCActcctc                                                                                                       >  2:514221/1‑251 (MQ=255)                                                                                                                                             aaCGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGCGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAACACTGCCGATGATCAATCCTCCGACg                                                                                                  >  2:425008/1‑251 (MQ=255)                                                                                                                                             aaCGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACg                                                                                                  >  2:179941/1‑251 (MQ=255)                                                                                                                                                   tGCTGATTTTCGCCTGACAGTAGCGGAAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGAGCGCACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGCTGATCAATCCTCCGACGGTTgac                                                                                            >  2:113884/1‑249 (MQ=255)                                                                                                                                                        aTTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGAc                                                                                       <  2:257496/251‑1 (MQ=255)                                                                                                                                                          tttCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGCGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACgg                                                                                     >  1:51049/1‑251 (MQ=255)                                                                                                                                                                    cAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATGGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGt                                                                           <  2:314043/251‑1 (MQ=255)                                                                                                                                                                           ggCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTAt                                                                    >  1:466207/1‑251 (MQ=255)                                                                                                                                                                            gCAGCGCCCTTTTAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGTTAGCGCAAGTCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTAtt                                                                   <  1:577141/251‑1 (MQ=255)                                                                                                                                                                                      ttGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCAGGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGaaaaa                                                          <  1:99080/250‑1 (MQ=255)                                                                                                                                                                                      ttGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTACCCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGaaaaa                                                          <  1:300643/250‑1 (MQ=255)                                                                                                                                                                                          gCATCATCCAGGAGAACAGCGTAATATTGGCACGTATCCGGATGGTCTGCTTTCAATGCGGGCAGTGCGAACTCGGTAGCGCAAGGCTTATTTTTTACTGCGCACTCGACATTTCACCGTGAAACACCGCCATTTGTAGAAGCATCATGATTGGGTAACGCCAAATCGCCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCa                                                     <  2:466207/251‑1 (MQ=255)                                                                                                                                                                                          gCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCa                                                     <  2:135572/251‑1 (MQ=255)                                                                                                                                                                                            aTCAGCCAGGAGAACAGCGGAATATTTTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACGGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTTTTGGCGCAGCATGATTGGGTAACGCCAAATCACCTCATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTAGTGTTGAAAAACATCaaa                                                   <  2:414232/251‑1 (MQ=255)                                                                                                                                                                                            aTCAGCCAGGAGAAAAGCGGAATATTGTCACGTTGCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGAAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAGCTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACGTCaaa                                                   <  2:619556/251‑1 (MQ=255)                                                                                                                                                                                                  cAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCAt                                             >  1:473177/1‑251 (MQ=255)                                                                                                                                                                                                     gggaACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGt                                          <  1:158798/249‑1 (MQ=255)                                                                                                                                                                                                       gaACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGtt                                         >  1:82126/1‑250 (MQ=255)                                                                                                                                                                                                         aCAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACTCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTTAAACTCCGCCAGTTGTTGGCGCAGCAGGATTGGGTAACGCCAAAGCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGATGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTaaa                                      <  1:565415/251‑1 (MQ=255)                                                                                                                                                                                                              ggAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGATTATTTTTGACTGCGCACTCGCCATTTCACCGTGACACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTCCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAACCATCAAAAAGCATGGTTAAAGgttg                                 >  1:295611/1‑251 (MQ=255)                                                                                                                                                                                                                aaTAGTGTCACGTATCCGGATGTGCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCTGAGGGCTTATTGTTGTCTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCCCCACAGTCAAAAAGCTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGTGTATGCCGGAAAAACATCAAAAAGCATGGTTAAAGgttgtt                               <  2:31738/251‑1 (MQ=255)                                                                                                                                                                                                                      gTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATGTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGGTTGTTtgttgt                         <  1:247702/251‑1 (MQ=255)                                                                                                                                                                                                                          cGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGGTTGTTTGTTGTAACt                     <  1:429829/251‑1 (MQ=255)                                                                                                                                                                                                                           gTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGTTTGTTTGTTGTAACtt                    <  1:137180/251‑1 (MQ=255)                                                                                                                                                                                                                                   aTGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATTGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGGTTGTTTGTTGTAACTTGCTGGATg            >  2:150927/1‑251 (MQ=255)                                                                                                                                                                                                                                   aTGGTCTGCTGACAATGCGGGCAGTGCGAACGCGGGAGCGCAAGGCTTATTTTTGTCTGCGCACTCGACATTTCACCGGGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGGTTGTTTGTTGTAACTTGCTGGATg            <  2:228882/251‑1 (MQ=255)                                                                                                                                                                                                                                          gctgACAATGCGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGGTTGTTTGTTGTAACTTGCTGGATGCGGCGt      <  2:24178/250‑1 (MQ=255)                                                                                                                                                                                                                                             gACAATGCGGGCATTGGGAACGCGTTAGCTCAAGGTTTATTGTTGACGTCGCACTCGACATTTCACCGTTAAACTCCGCCAGTTGTTGGCGCAGCGTGATTTGGTAACGCCAAATTACCACATTCAAAAAACGGCCGATGATCAATCCTCCGACTGTGGCCAGGACGGGCATCGCCGTGGTGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGTTTGGTTGTTGTATCTTGCTGGATgcgttttaaaa  <  2:570379/251‑9 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    CCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGGTAACGCCAAATCACCACATTCAAAAAACTGCCGATGATCAATCCTCCGACGGTTGCCAGGACGGGCATCGCCGTGGGGTATTGCTGAAAAACATCAAAAAGCATGGTTAAAGGTTGTTTGTTGTAACTTGCTGGATGCGGCGTAAAA  >  NC_000913/3113909‑3114395

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