Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,331,507	C→T	G83D (GGT→GAT)	obgE ←	GTPase ObgE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,331,507	0	C	T	96.8%	199.9 / NA	62	G83D (GGT→GAT)	obgE	GTPase ObgE
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (26/34);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTACCGCCATCACCTGCAACGACCAGAATCGATGCTTCATCAACAAAC  >  NC_000913/3331261‑3331749                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aaCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATcgcg                                                                                                                                                                                                                                                <  1:357418/251‑1 (MQ=255)      aCCGTCGTCAGCCAGCAGCATCAGCTCCCGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGCAACGGGTATTTCCCAGACCGTGCCAGCCGCCCTTAGCATCCAGCAGACGCTCACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTTCCCTGGTCTATTACACGCTTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTAc                                                                                                                                                                                                                                           <  1:238493/251‑1 (MQ=255)         gACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCgg                                                                                                                                                                                                                                        >  2:335718/1‑251 (MQ=255)               agccagcAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGt                                                                                                                                                                                                                                  <  1:333253/251‑1 (MQ=255)               agccagcAGCAGCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGt                                                                                                                                                                                                                                  <  2:288886/251‑1 (MQ=255)                gccagcAGCATCAGCTCCAGCGGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTGAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTGCCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTc                                                                                                                                                                                                                                 <  2:47164/251‑1 (MQ=255)                     cagcaTCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATAGCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGc                                                                                                                                                                                                                            <  2:202574/251‑1 (MQ=255)                      agcaTCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCt                                                                                                                                                                                                                           >  1:119830/1‑251 (MQ=255)                           cAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACTTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTgcg                                                                                                                                                                                                                       <  1:189299/250‑1 (MQ=255)                                        agcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACcgcg                                                                                                                                                                                                         >  1:317917/1‑251 (MQ=255)                                             gcgcTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTg                                                                                                                                                                                                    >  1:215258/1‑251 (MQ=255)                                                     cGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGaa                                                                                                                                                                                             <  2:101683/250‑1 (MQ=255)                                                      gCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGaa                                                                                                                                                                                             >  2:280439/1‑249 (MQ=255)                                                       cccGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGa                                                                                                                                                                                          <  2:127049/251‑1 (MQ=255)                                                         cGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGAtt                                                                                                                                                                                        >  2:183085/1‑251 (MQ=255)                                                                 cGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCaaaa                                                                                                                                                                                <  1:280439/251‑1 (MQ=255)                                                                  gTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                                                                                                                                                                               >  1:392986/1‑251 (MQ=255)                                                                      gTTTTCGGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTGTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATa                                                                                                                                                                           <  1:403297/251‑1 (MQ=255)                                                                          tCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCg                                                                                                                                                                       >  2:197573/1‑251 (MQ=255)                                                                            tGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGAt                                                                                                                                                                     >  1:164055/1‑251 (MQ=255)                                                                                   gtgtACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCgtg                                                                                                                                                              <  1:49084/251‑1 (MQ=255)                                                                                     gtACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGtt                                                                                                                                                            <  2:60309/251‑1 (MQ=255)                                                                                      tACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGtt                                                                                                                                                            <  2:215258/250‑1 (MQ=255)                                                                                       aCGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGATTACCGGCACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAACACGATAATCGATAAGCGTGTTc                                                                                                                                                           >  1:409424/1‑250 (MQ=255)                                                                                        cGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCTTGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCGGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAg                                                                                                                                                         <  1:236862/251‑1 (MQ=255)                                                                                              aCGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGTTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCTCGGAAACATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtctc                                                                                                                                                   <  2:44302/251‑1 (MQ=255)                                                                                              aCGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtctc                                                                                                                                                   >  1:292349/1‑251 (MQ=255)                                                                                                                    tGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTc                                                                                                                              >  1:335196/1‑250 (MQ=255)                                                                                                                        cAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTcaccac                                                                                                                         >  2:169287/1‑251 (MQ=255)                                                                                                                         aGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGCCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCaccacc                                                                                                                        <  1:140689/251‑1 (MQ=255)                                                                                                                               tGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCAcc                                                                                                                  >  1:83374/1‑251 (MQ=255)                                                                                                                                    gccgccCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGt                                                                                                             >  1:281547/1‑251 (MQ=255)                                                                                                                                      cgccCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTcg                                                                                                           <  1:281011/251‑1 (MQ=255)                                                                                                                                             aGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCTATTACACGCGTACCTACCGGCACTTTAATCGCCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGTGTGTTCATGTTCCCGTCGGCTTCCATCTATACGTCACCACCATCACCGCCGTCGCCGCCAt                                                                                                    <  1:23854/251‑1 (MQ=255)                                                                                                                                             aGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTTACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCCTACGTCACCACCATCACCGCCGTCGCCGCCac                                                                                                    >  2:115832/1‑250 (MQ=255)                                                                                                                                             aGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATCCGTCACCACCATCACCGCCGTCGCCGCCAt                                                                                                    >  2:121558/1‑251 (MQ=255)                                                                                                                                             aGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCAt                                                                                                    >  1:189423/1‑251 (MQ=255)                                                                                                                                              gCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATc                                                                                                   >  1:50022/1‑251 (MQ=255)                                                                                                                                                    agcagACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGgcc                                                                                             <  1:249407/251‑1 (MQ=255)                                                                                                                                                       agACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTcgccg                                                                                                        >  1:31807/1‑237 (MQ=255)                                                                                                                                                       agACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTcgccg                                                                                                        <  2:31807/237‑1 (MQ=255)                                                                                                                                                        gACGCTGACCGTTTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTCCCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCt                                                                                         <  1:169287/251‑1 (MQ=255)                                                                                                                                                                      ttGGTCATATCGCCCTTTTTTTCACCAGTCCCCTGGTCGATTACACGCGTTCCTACCTGCCCTTTAATCGTCACGTCTGTATCGCGCTTCCCGTTCCAGTCGCGGCTTGCTCCATTCGGACCTGGCTCTGCACGGAAAGATTTTTCAAAACTATAATCGATAAGCGTGTTCGGGTTCTCGTGGGCTTCCATCCATCCGTCACCACAGTCACCGCCGGCGCCGCCATCCGGGCCGCCTTTCGGAATATAcg                                                                            <  2:3221/250‑2 (MQ=255)                                                                                                                                                                      ttGGTCATATCGCCAATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTGTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACAATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACt                                                                            <  1:115832/250‑1 (MQ=255)                                                                                                                                                                        ggTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCCCCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACttt                                                                          <  1:55890/250‑1 (MQ=255)                                                                                                                                                                        ggTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACtttt                                                                         >  2:294215/1‑251 (MQ=255)                                                                                                                                                                           cATATCTCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTcg                                                                       <  1:117961/250‑1 (MQ=255)                                                                                                                                                                              atCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGCTAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCTGCCATCCGGGCCGCCTTTCGGAATATACTTTTCgcggc                                                                   >  1:176947/1‑251 (MQ=255)                                                                                                                                                                              atCGCCCAGGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACTTCTTTCTCTCGCTTACCGGTACATTCGCGGCTTGCGCCATGCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCAATAAGCGTGTTCAGTTTCTCGTCGGCTCCCATCCATACGTCACCTCCATCCCCGCCGTCGCCGCCATCCTGGCCGCCTTTCGGAATATACTTTTCgcggc                                                                   <  1:237099/251‑1 (MQ=255)                                                                                                                                                                                     aTGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCt                                                            <  1:335718/251‑1 (MQ=255)                                                                                                                                                                                          ttCACCAGTACCCTTGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCTCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACg                                                        <  1:32372/250‑1 (MQ=255)                                                                                                                                                                                                 gtccccTGGTCGATTACACGCGTACCTACCGGCACTTTAATCTTCACGTCTTTCTCGCGCTTACCGGTCCAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTGCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCAtt                                                <  2:400512/248‑1 (MQ=255)                                                                                                                                                                                                 gTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCAtt                                                >  1:403486/1‑251 (MQ=255)                                                                                                                                                                                                   ccccTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTAc                                              <  1:121786/250‑1 (MQ=255)                                                                                                                                                                                                     ccTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTACCg                                            <  1:379599/251‑1 (MQ=255)                                                                                                                                                                                                                        ccctcccgcccttttAATCGTCATGCCTTTACCGCCTTTACCGGTAAAGTCGCGGCTTGCGCCTTTCTGACCGCGCTCTGCCCGAAAAGATTTTTCAAAACGAAAATCGATAAGCTTGTACAGGTTCTCGTCGGCTTCCATCCATCCGCCCCCACCACCCCCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTACCGCCATCACCTGCAACGACCa                         <  2:292349/239‑1 (MQ=255)                                                                                                                                                                                                                            accGGCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTACCGCCATCACCTGCAACGACCAGAAt                     >  1:318715/1‑251 (MQ=255)                                                                                                                                                                                                                               ggCACTTTAATCGTCACGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCCCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCCCCCAACCCTTCCCGCCATCACCTGCAACGACCAGAATCg                   >  2:55054/1‑250 (MQ=255)                                                                                                                                                                                                                                     ttAATCGTCACTTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGTCTCCCATCCATACGTCACCCCCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATCCTTTTCGCGGCGGAAGCTCACGCAACCATTACCGCCATCACCTGCAACGACCAGAATCGATGc               <  1:121558/248‑1 (MQ=255)                                                                                                                                                                                                                                         ccGTCACGTCTTTATCGCGCTTACCGGTACGGTCGCGGCTTGCGCCATTCTGACCCCGCTCTGCACGAAAAGATTTTTCAAACCAATAACCGAAAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCACCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTACCGCCATCACCTGCAACGACCAGAATCGATGCTtcatca        <  2:83374/250‑1 (MQ=255)                                                                                                                                                                                                                                               cGTCTTTATCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTACCGCCATCACCTGCAACGACCAGAATCGATGCTTCATCAACAAAc  >  1:204945/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tCTTTATCGCGCTTACGGGTACGTTCGCGGCTTGCCCCATCCGGCCCGCGCGCGGCAGGGAAAGATTTTTCAAAACGATAATGGATAAGCGTGTTCAGGTTCTCGGCGGCTTCCATCCATACGTCCCCACCATCTCCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCATCCATTACCGCCATCACCTGCAACGACCAGAATCGATGCTTCATCAACAAAc  <  1:167385/249‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   AACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCCCTTAGCAACCAGCAGACGCTGACCGTGTTTGGTCATATCGCCCATGGTTTCACCAGTACCCTGGTCGATTACACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGTCGGCTTCCATCCATACGTCACCACCATCACCGCCGTCGCCGCCATCCGGGCCGCCTTTCGGAATATACTTTTCGCGGCGGAAGCTCACGCAACCATTACCGCCATCACCTGCAACGACCAGAATCGATGCTTCATCAACAAAC  >  NC_000913/3331261‑3331749

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