Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	793,557	G→A	L78F (CTC→TTC)	modF ←	ABC family protein ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	793,557	0	G	A	97.7%	139.7 / NA	43	L78F (CTC→TTC)	modF	ABC family protein ModF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/24);  minor base T (0/1);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCAAAATTTGCAACGATGACATTTTATATCCCAC  >  NC_000913/793311‑793800                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     caAGAATCAACAAGTCAGGCACCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGGCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAAttt                                                                                                                                                                                                                                                   <  1:262251/250‑1 (MQ=255)   aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTg                                                                                                                                                                                                                                                <  1:305026/251‑1 (MQ=255)        caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAAtt                                                                                                                                                                                                                                           <  1:215959/251‑1 (MQ=255)                cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAg                                                                                                                                                                                                                                   <  2:317880/251‑1 (MQ=255)                         aCATCAGCGCCTGAAACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACg                                                                                                                                                                                                                           >  1:326268/1‑250 (MQ=255)                          cATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTGATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCTAAGGAGAGACGAg                                                                                                                                                                                                                         <  2:198395/251‑1 (MQ=255)                                         cagcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑t                                                                                                                                                                                                         <  1:396773/251‑1 (MQ=255)                                         cagcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCAGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑t                                                                                                                                                                                                         <  1:226040/251‑1 (MQ=255)                                                   tttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATCCCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCg                                                                                                                                                                                               >  1:44655/1‑251 (MQ=255)                                                       gcgTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCa                                                                                                                                                                                           <  2:326268/251‑1 (MQ=255)                                                             cGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTc                                                                                                                                                                                     >  2:75328/1‑251 (MQ=255)                                                              gCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTc                                                                                                                                                                                     >  1:235763/1‑250 (MQ=255)                                                                            aTTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTcc                                                                                                                                                                      >  2:108609/1‑251 (MQ=255)                                                                            aTTTAAAGCGTCGGTCGAGCAGGGCGGTAATCCCGAACTCCTGTGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGCTCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGACGGAGAGCCGAGTGATGTGGGACAAC‑TGGCTTTGCCGTTCCCCTTTCACAAGCGGAAGTTcc                                                                                                                                                                      >  2:178447/1‑251 (MQ=255)                                                                              ttAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTcccc                                                                                                                                                                    >  1:203844/1‑251 (MQ=255)                                                                                       cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑Tgg                                                                                                                                                                                                       >  1:113110/1‑207 (MQ=255)                                                                                       cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑Tgg                                                                                                                                                                                                       <  2:113110/207‑1 (MQ=255)                                                                                        ggtggatgaggGCGGTAATCCCGAACTTCTGCGCCAGTTGCATGAAACGCGTTGCATCCTTGCCTTCATCCTGAATGATCTCAGCCGTAGTCCGTCCGGTGTCATCTCCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCGGAGGAATTTTGGCAATTGCGCGAAGGAGAGACGAGTGATG‑GGGAAAACCGGGCTTTGCTGTACTCCGTTCAAAAGCGAAATGTCCCCCGCCAgcgcg                                                                                                                                                          <  2:373454/244‑1 (MQ=255)                                                                                         gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAgcgcgc                                                                                                                                                         >  1:32915/1‑251 (MQ=255)                                                                                             gggaggGCGGTAATACCGAACTGCTGTGCCAGTTCGATGCACCGCGGTGCTCCCTTTCCTTCATCCTTACTGCTCTCAGCCGTGGGGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTTTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCAAAGGAGAGACGAGTGATGTGGGAAACC‑TGGCTTTGCCGTTCCCCTTTCAAATGCGGAAGTCCCCCCGCCAGCGCGCGGGc                                                                                                                                                     <  2:329229/250‑1 (MQ=255)                                                                                                    cGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGcc                                                                                                                                              >  1:16215/1‑251 (MQ=255)                                                                                                        aaTACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGAtt                                                                                                                                          <  1:75328/251‑1 (MQ=255)                                                                                                                   gctgcGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCa                                                                                                                               <  1:182592/251‑1 (MQ=255)                                                                                                                    ctgcGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGTTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCAt                                                                                                                              >  2:280049/1‑251 (MQ=255)                                                                                                                      gcgcCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCAtt                                                                                                                             <  1:139941/250‑1 (MQ=255)                                                                                                                      gcgcCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTc                                                                                                                            <  1:284422/251‑1 (MQ=255)                                                                                                                                 atgcaACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGAATATCGGTGTTATTCCGCTGCCATTCGTGGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGTTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACaaa                                                                                                                 <  1:62575/251‑1 (MQ=255)                                                                                                                                        gcgGTGCATCCTTTCCTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCaa                                                                                                          <  1:332843/251‑1 (MQ=255)                                                                                                                                         cgGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAAc                                                                                                         >  2:77375/1‑251 (MQ=255)                                                                                                                                                                gAATGATCTCGTCCGTAGTGCTTCCGGTGTCATCTGCGCCAGGGCCGAGAATTGCTGTTTTTTTCCGCTGCCATTCGTCGCTGACAAATTTTTGCAATCGCTTGAAGGAGAGACGGTTGAGGTGGGAAAAC‑TGGCTTTGCGGTTCGCCTTCGAAAAGCGGAATTGCCCCCGCCCGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGAGAAACGCCCAACTATCACCCGCGTTTAGCGTTa                                                                                    <  2:412155/249‑1 (MQ=255)                                                                                                                                                                      tctcAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGc                                                                            >  1:410509/1‑251 (MQ=255)                                                                                                                                                                             cGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGc                                                                      >  1:336794/1‑250 (MQ=255)                                                                                                                                                                                        cGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTtgtg                                                          <  1:77375/251‑1 (MQ=255)                                                                                                                                                                                         ggTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCAACTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATCCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTgtgt                                                         >  1:59551/1‑251 (MQ=255)                                                                                                                                                                                                   tCGCCAGGGCGGAGCATATCGGGTTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGGGAGCCGAGTGAGGTGGGAAAAC‑TGGCTTTGCCGTTCCCCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAg                                                  <  2:339114/248‑1 (MQ=255)                                                                                                                                                                                                   tCGCCAGGGCCGAGCATATCGGTTTTATTCCGCTGCCTTTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCCCCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAg                                                  <  2:336794/248‑1 (MQ=255)                                                                                                                                                                                                                    atCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCCGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCaaa                               >  2:212849/1‑250 (MQ=255)                                                                                                                                                                                                                    atCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑GGGCTTTGCCGTTCCCCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCaaa                               <  2:203844/250‑1 (MQ=255)                                                                                                                                                                                                                                 cGCTGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCGGAGGCAATTGCAGCGTTTTTGTTTCGCTAAGACGACACGTGCCTTGCAAAATTTTCAACGATg                  >  2:366496/1‑250 (MQ=255)                                                                                                                                                                                                                                    tGCCATTCGTCGCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCAAAATTTGCAACGATGACAt              >  2:45376/1‑251 (MQ=255)                                                                                                                                                                                                                                          tcgtcgCTGACGAATTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑GGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCAAAATTGGCAACGATGACATTTtat         <  1:108609/250‑1 (MQ=255)                                                                                                                                                                                                                                            gtcgCTGACGAAGTTTGACAATTGCGCGAAGGAGAGACGAGTGAGGTGGGAAAAC‑TGGCTTGGCCGTGCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGACTTGCAAAATTTGCAACGATGACATTTtatagc      <  2:175514/251‑3 (MQ=255)                                                                                                                                                                                                                                                  gACGATTTTTTGCAATGTCTCGAAGGAGAGGCGAGGTATGGGGAAAAAC‑GGGTTTGGCCGTTCACCTTACAAAAGCGGAAGTCCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCAAAATTTGCAACGATGACATTTTATATCCCAc  <  2:26498/249‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     CAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGGGAAAAC‑TGGCTTTGCCGTTCACCTTTCAAAAGCGGAAGTTCCCCCGCCAGCGCGCGGGCCAGGGCCGATTTCCCGCTTCCATTCGAACCGACAAACGCCCAACTATCACCCGCGTTTAACGTTAGCTGAGGCAATTGCAGCGTTTTTGTGTCGCTAAGACGAAACGTGCCTTGCAAAATTTGCAACGATGACATTTTATATCCCAC  >  NC_000913/793311‑793800

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