Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	851,631	A→G	Y323H (TAC→CAC)	opgE ←	phosphoethanolamine transferase OpgE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	851,631	0	A	G	98.0%	173.1 / NA	50	Y323H (TAC→CAC)	opgE	phosphoethanolamine transferase OpgE
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (28/21);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGCTTGCGCAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑ACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGG  >  NC_000913/851387‑851873                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   agcatGCGCAGGAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTAGTTTCCTTGTGTCCTGTCGCAGGCCGGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCTGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑gcacatc                                                                                                                                                                                                                                              <  2:119271/247‑4 (MQ=255)       cgcAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGACCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGc                                                                                                                                                                                                                                        >  2:108627/1‑251 (MQ=255)               aTCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATTAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATGACCTTCTTTCAGTAAGT‑GCACTTCATCTGCTCGTTTgg                                                                                                                                                                                                                                <  1:278602/251‑1 (MQ=255)                 cGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCGGCTCGTTTGGc                                                                                                                                                                                                                               <  2:254024/250‑1 (MQ=255)                      gTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGCTGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCt                                                                                                                                                                                                                         >  2:74550/1‑251 (MQ=255)                         tGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTgg                                                                                                                                                                                                                       <  2:285651/250‑1 (MQ=255)                               aTGGTTTAGATATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGACATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCGTCAAAATGACCTTCTTTCATGAAGT‑GCACTTCATCTACTCGTTTGGCGATGCTGGCGATAGc                                                                                                                                                                                                                <  2:187199/251‑1 (MQ=255)                                      agagATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTCTCg                                                                                                                                                                                                         >  1:324756/1‑251 (MQ=255)                                          aTAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGTTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTAtt                                                                                                                                                                                                     <  1:91978/251‑1 (MQ=255)                                                cgacgTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCg                                                                                                                                                                                                >  2:135469/1‑250 (MQ=255)                                                       tctttcGATAGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCCGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGACTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGAcc                                                                                                                                                                                        >  2:185610/1‑251 (MQ=255)                                                         tttcGATTGCACACAGGTTACGTCTTTTCCTTGCGTCCTGTCGCAGGCCTGCGGATGTGCGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCCGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCt                                                                                                                                                                                      >  2:208991/1‑251 (MQ=255)                                                             gATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGtt                                                                                                                                                                                  >  2:17281/1‑251 (MQ=255)                                                                 gCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTGCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGa                                                                                                                                                                               <  1:289335/250‑1 (MQ=255)                                                                     aaaGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAacc                                                                                                                                                                          >  1:10434/1‑251 (MQ=255)                                                                      aaGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAacc                                                                                                                                                                          >  2:94980/1‑250 (MQ=255)                                                                            tCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGCTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGCAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTTACCCTGTTTGGAAAACCACCAGGt                                                                                                                                                                   >  2:143678/1‑251 (MQ=255)                                                                             cGTATTTTCCTTGGGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGt                                                                                                                                                                   <  2:110740/250‑1 (MQ=255)                                                                                      cgtgtgTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGAGGTAGAACAATCAGCTGCGGTGGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGGAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCGGGCGATAGCGGGATCGTAGTCGCCGAATTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGccc                                                                                                                                                         <  2:212608/249‑1 (MQ=255)                                                                                         gtgtCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCcgcg                                                                                                                                                      >  2:168329/1‑251 (MQ=255)                                                                                              cTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGAtt                                                                                                                                                 >  2:50095/1‑251 (MQ=255)                                                                                                gTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTGTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGAttt                                                                                                                                                <  1:143134/250‑1 (MQ=255)                                                                                                   gCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGcc                                                                                                                                             >  1:121928/1‑250 (MQ=255)                                                                                                            cGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACa                                                                                                                                    >  1:403083/1‑250 (MQ=255)                                                                                                                 gtgAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCCGGGTGACAAAGttt                                                                                                                              >  2:16249/1‑250 (MQ=255)                                                                                                                 gtgAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGtt                                                                                                                               >  1:312584/1‑250 (MQ=255)                                                                                                                                   gAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGAACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTg                                                                                                            >  1:196091/1‑251 (MQ=255)                                                                                                                                      cAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACATGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGcc                                                                                                          >  1:172591/1‑250 (MQ=255)                                                                                                                                           aGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGAc                                                                                                    >  2:140272/1‑251 (MQ=255)                                                                                                                                                       gAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTg                                                                                          >  1:82062/1‑249 (MQ=255)                                                                                                                                                         gTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCg                                                                                      <  1:135469/251‑1 (MQ=255)                                                                                                                                                                 gcgcCAGCACTTGAGCGGTCATAGCCAGTAACGCTGCGTCTTTGGTGTTTTTATCTGCTTCAAAATTTCCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTTGAGAACCGCCAGGTCTGGAGGCGCTCGCGATTTGCCAGGGTGACAAAGGTAGCCTGAAGCTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGcc                                                                               <  1:181872/250‑1 (MQ=255)                                                                                                                                                                  cggcAGCACTTGAGCGGTCAATTCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCGTCAAAATTACCTGCTTTCAGGAAGT‑GCACTTCAGCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCaa                                                                             <  2:104158/248‑1 (MQ=255)                                                                                                                                                                  cgcCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCaa                                                                             <  2:242251/251‑1 (MQ=255)                                                                                                                                                                             gAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGc                                                                  <  1:248885/251‑1 (MQ=255)                                                                                                                                                                              aGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCg                                                                 >  1:349257/1‑251 (MQ=255)                                                                                                                                                                                     atatCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGcc                                                          <  1:168329/251‑1 (MQ=255)                                                                                                                                                                                                 gtttcgtcttggggtgttttCTCTGCTTTAAAATTACCTGCTTTCAGGAAGTGACA‑TTCAGCTGCTCGTTTGGCGAGGCTGGCGATAGCGGTATCGTATTCGCCGAGTTGACCCTGGTTGGACACCCACCAGGTCTTGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTGATCCTGAAACTGTGGTTTGCCATCGACAACGTGATTGAGCGGTAAGCCAAGCGATTTCTGCGTGGAGCCCCTCGCGGCAAt                                              <  2:130523/235‑1 (MQ=255)                                                                                                                                                                                                  cTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGAGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATg                                             >  2:267831/1‑251 (MQ=255)                                                                                                                                                                                                    tCGTCTTTGGAGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGt                                            >  1:341734/1‑250 (MQ=255)                                                                                                                                                                                                       tCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTACGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCACTGTAGTc                                        >  2:249293/1‑251 (MQ=255)                                                                                                                                                                                                           tGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGc                                     <  1:247705/250‑1 (MQ=255)                                                                                                                                                                                                                             cAAAATTGCCTTCTTGCAGAAAGT‑GCGCTTCATCTGCTCGTTTGGCGATGCTGGCTATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGTCAAAGTTATCCGGAAGCTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGAATTCTGCGGGGAGCCACTGGCGGCAATGTACTCAGCAAATATCAAACCGTTAAc                   <  2:294790/250‑1 (MQ=255)                                                                                                                                                                                                                             cAAAATTACCTTCTTTCAGGAAGT‑GCCCTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGTCCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCACACCGTTAAc                   <  2:312584/250‑1 (MQ=255)                                                                                                                                                                                                                              aaaaTTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGACAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACgc                 >  2:81649/1‑251 (MQ=255)                                                                                                                                                                                                                              aaaaTTACCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGAGATGCTGGCGATAGCGGTATCGTATCCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGGTTGCCATCGACAACGCGATTGAGCTTTAAGCCAAGGGCATTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATAAAACCGGTAACgc                 >  1:283639/1‑251 (MQ=255)                                                                                                                                                                                                                               aaaTTACCTTCTTTCAGGCAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACgct                >  1:136719/1‑251 (MQ=255)                                                                                                                                                                                                                                    aCCTTCTTTCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGctgctg            <  2:82062/250‑1 (MQ=255)                                                                                                                                                                                                                                          tttCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGAGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGCAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGCAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGACGCTGGCAAACg     >  2:39787/1‑251 (MQ=255)                                                                                                                                                                                                                                           ttCAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCGGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTAGCCTGAAACTGTGGTTTTCCATCGACAACTCGATTGAGCTTTAAGCCAATCGATTTCTGCGTGGAGCCACTGGCGGGAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCaaa       <  1:377505/248‑1 (MQ=255)                                                                                                                                                                                                                                             cAGGAAGT‑GCACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACgggg  <  1:74550/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   AGCTTGCGCAGTAAATCGTCCGTTTGCGTCATGGTATAGAGATAGCACGACGTTTCTTTCGATTGCACAAAGGTTTCGTATTTTCCTTGTGTCCTGTCGCAGGCCTGCGGATGTGAGCCCATCAGATGTAGAACAATCAGCTGCGGTTGCGAGTGCTCTTGCGCCAGCACTTGAGCGGTCATATCCAGTAACGCTTCGTCTTTGGTGTTTTTATCTGCTTCAAAATTACCTTCTTTCAGGAAGT‑ACACTTCATCTGCTCGTTTGGCGATGCTGGCGATAGCGGTATCGTATTCGCCGATTTGACCCTGGTTGGAAAACCACCAGGTCTGGAAGCCCGCGCGATTTGCCAGGGTGACAAAGTTATCCTGAAACTGTGGTTTGCCATCGACAACGCGATTGAGCGTTAAGCCAAGCGATTTCTGCGTGGAGCCACTGGCGGCAATGTAGTCAGCAAATATCAAACCGTTAACGCTGCTGGCAAACGGGG  >  NC_000913/851387‑851873

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: ‑