Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,789,146	G→C	V312V (GTC→GTG)	waaH ←	UDP‑glucuronate:LPS(HepIII) glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,789,146	0	G	C	97.7%	141.6 / NA	43	V312V (GTC→GTG)	waaH	UDP‑glucuronate:LPS(HepIII) glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (13/29);  minor base A (1/0);  total (14/29)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCCTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCATAATTTTGTCGGCATAATTTCGATTTAATTTCTCCAGCAGGCGGGTAATCTTAATATAGTGACGTTGATAATTAAGGTTTT  >  NC_000913/3788900‑3789390                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAATGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCt                                                                                                                                                                                                                                                  >  2:290308/1‑251 (MQ=255)  cAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATTGTGCGACCGCGTTTTGTCTCGCGATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCtt                                                                                                                                                                                                                                                 <  1:197092/251‑1 (MQ=255)                 gATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTGCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACGTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACAtt                                                                                                                                                                                                                                  <  2:334818/251‑1 (MQ=255)                   tAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTgg                                                                                                                                                                                                                                >  1:202015/1‑251 (MQ=255)                   tAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGACAGAGCACGCTGATCGTTTGCTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGGCACTTAACGCTAGCGCAGATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGATACCCACCTTCACGCGGCGCACATTgg                                                                                                                                                                                                                                >  1:250811/1‑251 (MQ=255)                           gTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAgg                                                                                                                                                                                                                        <  2:154754/251‑1 (MQ=255)                                     aTCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTataca                                                                                                                                                                                                              <  2:93561/251‑1 (MQ=255)                                                         aGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCg                                                                                                                                                                                          <  2:152727/251‑1 (MQ=255)                                                                   tCAGCCAGAGAACGCTGATCGTTTACTTCCAGAATCAGTGGCAGGGATGGGGCAACCAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCAATGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGAACTCAGGAGCTTATACATAACGGAAGTAAATATCTCGGCAAtcctcc                                                                                                                                                                                >  1:334069/1‑251 (MQ=255)                                                                      gCCAGAGAACTCTGATCGTTTACTTCCAGCATCAGTGGCCAGGCTGGGCCAAACAACCCGAAGCGGCGAACATTGACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGGCGCGCCATTGCCATCAGCGGACAGACCACATTACCTCCGGGCAACCCACACTCACGATGCGCACATTGGTAATCAGTCGCTTCTACCGACCGCAAGTAAATCTCTAGGCAATCATCCGtt                                                                                                                                                                             >  2:403628/1‑251 (MQ=255)                                                                       ccAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTg                                                                                                                                                                            >  1:85418/1‑251 (MQ=255)                                                                              aCGCTGATCGTTTACTGCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAGCCCACCTTCACGCTGCGCACATTGGTAATCAGTCGCGTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTa                                                                                                                                                                     <  1:275752/251‑1 (MQ=255)                                                                                      cGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATc                                                                                                                                                             >  2:129395/1‑251 (MQ=255)                                                                                           aCTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAAGGCGATGGTGCGACCGCGTTTTGTCGCGCCAGTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAAGATCCGGc                                                                                                                                                         <  1:371278/250‑1 (MQ=255)                                                                                                           gCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATg                                                                                                                                        <  2:306780/251‑1 (MQ=255)                                                                                                             aaGGCTGGGGCAAACAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGc                                                                                                                                      <  2:151302/251‑1 (MQ=255)                                                                                                                    gggCAAACAGCCCGAAGCGGCAAAACTTAACGCGAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACg                                                                                                                               <  1:307353/251‑1 (MQ=255)                                                                                                                           cAGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGct                                                                                                                        <  2:88205/251‑2 (MQ=255)                                                                                                                            aGCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGctt                                                                                                                       <  1:178189/251‑3 (MQ=255)                                                                                                                             gCCCGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTGCACGCTGCGCACATGGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGGATGGCAAACTCGCAATGcttc                                                                                                                      <  2:426992/251‑4 (MQ=255)                                                                                                                                cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATCCATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATcc                                                                                                                                                            >  1:86196/1‑210 (MQ=255)                                                                                                                                cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATcc                                                                                                                                                            <  2:86196/210‑1 (MQ=255)                                                                                                                                cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCaa                                                                                                                            <  1:123912/242‑1 (MQ=255)                                                                                                                                cGAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCaa                                                                                                                            >  2:123912/1‑242 (MQ=255)                                                                                                                                 gAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGc                                                                                                                                      >  2:268357/1‑231 (MQ=255)                                                                                                                                 gAAGCGGCGAAACTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGc                                                                                                                                      <  1:268357/231‑1 (MQ=255)                                                                                                                                           aacttAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCt                                                                                                        <  1:16008/248‑1 (MQ=255)                                                                                                                                                                aaaGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGt                                                                                                                                                                      >  2:72389/1‑168 (MQ=255)                                                                                                                                                                aaaGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGt                                                                                                                                                                      <  1:72389/168‑1 (MQ=255)                                                                                                                                                                aaaGGCGCTACGCGTAAAGCGATGGTGCGACCGCGTGTTGTCGCGCCATTGCCAGAAGCGGAAAGACCCCAGTAACGCCTGGTAACCCCCCGTCACGCCGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTTTTTGCGCACCGCATGGCAAACACGCAAGGCTTAGTAAGTGATTTGCTGATGAAATTCCGGATAAATc                                                                                    <  2:240472/250‑1 (MQ=255)                                                                                                                                                                 aaGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGTATAAATcgt                                                                                  <  1:11314/251‑3 (MQ=255)                                                                                                                                                                      gcTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCAACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGATCTTTGCGCAACGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAAttt                                                                             >  1:287523/1‑251 (MQ=255)                                                                                                                                                                         aCGCGTAATGCGATGGTGCGACAGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCTCCGCATGGCAAACACGCAATGCTTCGTAAGTCATTTGCTTATTAAATTCCGGATAATTCGTAATTTTg                                                                           <  2:157827/250‑1 (MQ=255)                                                                                                                                                                          cgcgTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTc                                                                         <  1:252910/251‑1 (MQ=255)                                                                                                                                                                                     aTGGGGCGACCGCGTTTTGTCGCGCCTTTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCAACGCATGGCAAACACGCAATGATTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTCGGCATAATTTc                                                              >  2:114158/1‑251 (MQ=255)                                                                                                                                                                                       ggTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATATTGTCGGCATAATTTCg                                                             <  2:255399/250‑1 (MQ=255)                                                                                                                                                                                           cGACCGCGTTTTGTCGCGCCATTGACATAAGCGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATGTTGTCGGCATAATTTCGATTTa                                                        <  1:129395/251‑1 (MQ=255)                                                                                                                                                                                                                          cGGAAAGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTCGGCATAATTTCGATTTAATTTCTCCAGCAGGCGGGTAATCTTAatat                          >  1:52115/1‑250 (MQ=255)                                                                                                                                                                                                                             aaaGACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTCGGCATAATTTCGATTTAAGTTCTCCAGCAGGCGGGTAATCTTAATATAGTg                      <  1:311199/251‑1 (MQ=255)                                                                                                                                                                                                                                gACCACAGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTCGGCATAATTTCGATTTAATTTCTCCAGCAGGCGGGTAATCTTAATATAGTGACg                   <  1:136369/251‑1 (MQ=255)                                                                                                                                                                                                                                      aGTAACGCCTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTCGGCATAATTTCGATTTAATTTCTCCAGCAGGCGGGTAATCTTGATATAGTGACGTTGATa             <  2:85418/251‑1 (MQ=255)                                                                                                                                                                                                                                             ccTGGTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTCGGCATAATTTCGATTTAATTTCTCCAGCAGGCGGGTAATCTTAATATAGTGACGTTGATAATTAAgg      <  2:120997/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gTAACCCACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCTTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCGTAATTTTGTCGGCATAATTTCGATTTAATTTCTCCAGCAGGCGGGTAATCTTAATATAGTGACGTTGATAATTAAGGtttt  <  1:312198/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCGCTTATACATACCGGAAGTAAATATCTCGGCAATCATCCGTTGGCGGGTAAGAATATCCGGCTCTTTGCGCACCGCATGGCAAACACGCAATGCCTCGTAAGTTATTTGCTGATGAAATTCCGGATAAATCATAATTTTGTCGGCATAATTTCGATTTAATTTCTCCAGCAGGCGGGTAATCTTAATATAGTGACGTTGATAATTAAGGTTTT  >  NC_000913/3788900‑3789390

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