Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,587,531	G→A	L195F (CTC→TTC)	ugpQ ←	glycerophosphodiester phosphodiesterase UgpQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,531	0	G	A	98.6%	217.2 / NA	72	L195F (CTC→TTC)	ugpQ	glycerophosphodiester phosphodiesterase UgpQ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (43/27);  minor base C (1/0);  total (45/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGAT  >  NC_000913/3587288‑3587777                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cagcagCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGCCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACtt                                                                                                                                                                                                                                                  >  2:245165/1‑250 (MQ=255)    cagcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATg                                                                                                                                                                                                                                               >  1:303858/1‑250 (MQ=255)        atatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTg                                                                                                                                                                                                                                           >  1:136401/1‑250 (MQ=255)          atGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTgag                                                                                                                                                                                                                                         <  1:109322/250‑1 (MQ=255)               gcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAAt                                                                                                                                                                                                                                   >  2:452363/1‑251 (MQ=255)               gcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGGTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAAt                                                                                                                                                                                                                                   <  2:521153/251‑1 (MQ=255)               gacGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAAGCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGa                                                                                                                                                                                                                                     <  1:308719/247‑1 (MQ=255)                      aggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGCCATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑AATCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACg                                                                                                                                                                                                                            >  1:253931/1‑251 (MQ=255)                      aggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACg                                                                                                                                                                                                                            >  1:350487/1‑251 (MQ=255)                         agaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGc                                                                                                                                                                                                                           >  1:468156/1‑249 (MQ=255)                          gaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGcagc                                                                                                                                                                                                                        >  2:167860/1‑251 (MQ=255)                                gCGACATGCTGAAGCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGcagcccagcc                                                                                                                                                                                                                  <  2:425055/251‑1 (MQ=255)                                  gaaaTGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCCTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCAAGTCAATCGCATCGGAACAGATGCAATCAACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGAAGTCCGGCGTCTTTCAACAGCATC‑ACTAGCGCTTTATTAAATAACTTATGATTGAGTTGAATAGAGACGaagcccagcagc                                                                                                                                                                                                                >  1:473713/4‑248 (MQ=255)                                  gACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCcgc                                                                                                                                                                                                                >  1:391670/1‑251 (MQ=255)                                    cATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTAAATCGCATCGGTCCAGCTGAAATCAACACCCCAGCGCAGCAACACTGCTGCGCGCTGGGGCTTCTTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCctcgc                                                                                                                                                                                                              >  1:570526/1‑251 (MQ=255)                                       gCTGAAACAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCACTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGCGTAATAGAGCCGCAGCCCAGCCGCGCggg                                                                                                                                                                                                           >  2:259095/1‑250 (MQ=255)                                                     gCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAgt                                                                                                                                                                                             <  2:177001/251‑1 (MQ=255)                                                          ccgcAGCCTGACATTCCGTTGCCAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGAGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGtcgtcg                                                                                                                                                                                        >  2:538116/1‑251 (MQ=255)                                                           cgcAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTcgc                                                                                                                                                                                       >  2:205818/1‑251 (MQ=255)                                                               gCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCa                                                                                                                                                                                   <  1:372809/251‑1 (MQ=255)                                                                ccTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCAc                                                                                                                                                                                  >  2:453275/1‑251 (MQ=255)                                                                        tCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGAtt                                                                                                                                                                                                                                            >  1:389279/1‑185 (MQ=255)                                                                        tCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGAtt                                                                                                                                                                                                                                            <  2:389279/185‑1 (MQ=255)                                                                                     aTTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATACACACCACAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGACGCGCCCCTCCTCCAGAACCAAACcgcac                                                                                                                                                             >  2:476211/1‑249 (MQ=255)                                                                                          gCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTACGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCACAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACAGCGCGGc                                                                                                                                                          >  1:377824/1‑249 (MQ=255)                                                                                             gTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTc                                                                                                                                                     >  2:387166/1‑251 (MQ=255)                                                                                             gTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAACGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTc                                                                                                                                                     >  1:325255/1‑251 (MQ=255)                                                                                               aaaGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACGCCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTcc                                                                                                                                                    <  1:13079/250‑1 (MQ=255)                                                                                                     cGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTgccg                                                                                                                                             >  2:146910/1‑251 (MQ=255)                                                                                                     cGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTgcc                                                                                                                                              <  1:270228/250‑1 (MQ=255)                                                                                                       gACCAATCACGTCAATCGCAGCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTGATGATTGAGATGAATAGAGACGCAGCCCGGCCGCGCGGTCAGTGCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGATgccgcc                                                                                                                                           <  2:239163/251‑1 (MQ=255)                                                                                                         ccAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCtg                                                                                                                                         <  1:205510/251‑1 (MQ=255)                                                                                                          cAATCACGTCAATCGCATCGGAACAGATGCCATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCAGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCCACTGCATC‑ACTCGCGCTTTATTGACTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCtgt                                                                                                                                        >  2:391770/1‑251 (MQ=255)                                                                                                            aTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGCCGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGAAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCAt                                                                                                                                                                              >  1:82711/1‑211 (MQ=255)                                                                                                            aTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGCCGTCGCGCCACTCAt                                                                                                                                                                              <  2:82711/211‑1 (MQ=255)                                                                                                                cGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGc                                                                                                                                   >  1:549290/1‑250 (MQ=255)                                                                                                                 gTCAAGCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACGCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAg                                                                                                                                 <  2:465488/251‑1 (MQ=255)                                                                                                                   cAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGG‑TTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATa                                                                                                                                                                                                                                  >  2:504346/1‑151 (MQ=255)                                                                                                                   cAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGG‑TTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATa                                                                                                                                                                                                                                  <  1:504346/151‑1 (MQ=255)                                                                                                                     ctcgcatcgatcccgataccaTCCAAACCCCAGCCCAGCAAATCTGCAGCGCGATCGGGTTTATTGACGGTCTAAACCCGAATACGCAGTCCGGCGTCTTTAAACTGCATC‑ACTCGCCCTTTATTGAATAACTTGTGATTGAGATGATTAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGCCGCGCCGCTCATCCAGCAACAACCCGCGCGGCCGTTCCGCTGCCGCCTGCTGTTCCGCTTc                                                                                                                             >  2:87893/21‑251 (MQ=255)                                                                                                                      tCGCATCGGTACAGATGCACTCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCt                                                                                                                            >  1:498731/1‑251 (MQ=255)                                                                                                                          aTCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTaa                                                                                                                          >  1:70187/1‑249 (MQ=255)                                                                                                                             ggTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGATTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGAAGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCAt                                                                                                                     >  2:307400/1‑251 (MQ=255)                                                                                                                                    aTGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAAtctc                                                                                                              <  1:185194/251‑1 (MQ=255)                                                                                                                                              cacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGAc                                                                                                    <  1:300101/251‑1 (MQ=255)                                                                                                                                                     gcgcAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAgcg                                                                                             <  2:416324/251‑1 (MQ=255)                                                                                                                                                            aaCTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGGGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAgcggcggcgc                                                                                      >  2:234125/1‑250 (MQ=255)                                                                                                                                                                ctgctgCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑CCTCGCGCTTTATTGCATAACTTCTGATTGAGCTGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCAGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGCTGACAGCAGCGGCGGCGTCAt                                                                                   >  1:64553/1‑250 (MQ=255)                                                                                                                                                                            ggggTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGGCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCcaca                                                                        <  1:198973/249‑1 (MQ=255)                                                                                                                                                                              ggTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTc                                                                    >  1:508463/1‑251 (MQ=255)                                                                                                                                                                                  tATTGACGGTATAAAACAGCATCCGCAGTCCGGCGTCTTTCAACTGCATA‑ACTCGCGCTTTATTGAATAACTTATGCTTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACCGCAGCGGCGGCGTCAAACCTGCCCACAGCTCGCGt                                                                >  1:195901/1‑251 (MQ=255)                                                                                                                                                                                     tGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACACACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACCGCTCGCGTgcc                                                             >  2:581509/1‑251 (MQ=255)                                                                                                                                                                                      gACGGTATAAACCATCATCCGCAGTCCGGCGTCTTTCAACTGCATC‑CCTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTgcc                                                             <  2:324006/250‑1 (MQ=255)                                                                                                                                                                                       aCGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCCGTTCCGGTGCCGCCTGTTGTGCCGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTgccgc                                                           >  1:315480/1‑251 (MQ=255)                                                                                                                                                                                         ggTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTGCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGt                                                                                                                                                 <  2:380291/163‑1 (MQ=255)                                                                                                                                                                                          gTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGAGTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCtt                                                                                                                              <  1:353250/181‑1 (MQ=255)                                                                                                                                                                                          gTATAAACCAGAATCAGCAGTACGGAGTCTTTCAACTGCCTC‑ACTCGAGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCCGTTCCGGTGCCGCCTGTTGTGCAGCtt                                                                                                                              >  2:353250/1‑181 (MQ=255)                                                                                                                                                                                                aCCAGAATCCGCAGTCCGGCGTCTTTCAACTGAATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTGCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGccac                                                  <  2:353356/251‑1 (MQ=255)                                                                                                                                                                                                     aaTCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCAttt                                             <  1:234125/251‑1 (MQ=255)                                                                                                                                                                                                            aGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGCCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGGTACGGTGCCGACTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTATTACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTAGccact                                      >  1:32236/1‑248 (MQ=255)                                                                                                                                                                                                                      cTTTCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccg                            >  1:528834/1‑251 (MQ=255)                                                                                                                                                                                                                        ttCAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTAGCGCCAGTAGTCGCGCCACTCATCAAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggt                          >  2:49666/1‑251 (MQ=255)                                                                                                                                                                                                                          cAACTGCCTC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCAAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgc                        <  2:126120/251‑1 (MQ=255)                                                                                                                                                                                                                          cAACTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtg                         <  2:391670/250‑1 (MQ=255)                                                                                                                                                                                                                           aaCTGCATCGAC‑CGCGCTTTAATGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCATCAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAAGCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgcc                       <  2:178947/251‑1 (MQ=255)                                                                                                                                                                                                                            aCTGCATC‑ACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGATCGCGCAAGTCGTCGCGCCACTCATCCCGCAACAAACCGCGCGGCAGTTACGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgacg                      >  1:338956/1‑248 (MQ=255)                                                                                                                                                                                                                            aCTGCATC‑ACTAGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCTCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCCCCATTTTGCCCGTTAATGGCccggtgccg                      >  2:479512/1‑251 (MQ=255)                                                                                                                                                                                                                                      cTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGgggtgggt             >  2:161150/1‑251 (MQ=255)                                                                                                                                                                                                                                      cTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCAACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATAGCCCGGCGCCGGTGGTgggg             >  2:430868/1‑250 (MQ=255)                                                                                                                                                                                                                                        cgcgCTTTATTGAATAACTTATGATTGAGATGAATAGACACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATAACGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGGTAATGGCCCGGTGCCGGGGGGGGGttt           >  2:318554/1‑251 (MQ=255)                                                                                                                                                                                                                                           gcTTTATGGAATAACTTATGATTGAGATGAATAGAGCCGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGAt        <  2:350487/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ttGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGat  >  1:459239/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATC‑ACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGAT  >  NC_000913/3587288‑3587777

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