Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	9,668	C→T	G121G (GGC→GGT)	mog →	molybdopterin adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	9,668	0	C	T	96.8%	180.8 / NA	64	G121G (GGC→GGT)	mog	molybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (23/38);  minor base G (2/0);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCT‑GGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCG  >  NC_000913/9426‑9906                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               gAACTGGAAACCCGCTTAATCCCCGATGAGCAGTCGATCATCGAGCAAACGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCCCGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTc                                                                                                                                                                                                                                         <  1:193141/250‑1 (MQ=255)      ggAAACCCGCTTAATCCCCGATGAGCAGGCGAGCATCGAGCAAACGTTGTGTGAGGT‑AGTGGAGGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCaaa                                                                                                                                                                                                                                    <  2:188450/250‑1 (MQ=255)                      cccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAAAGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCAAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGTTTCGCAAACAGGCGCTGATCCTTa                                                                                                                                                                                                                    >  2:384946/1‑250 (MQ=255)                        cGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGGGTGATTCGCAAACAGGCGCTGATCCTTAACt                                                                                                                                                                                                                 >  2:112648/1‑251 (MQ=255)                               cAGGCGATCATCGAGCAAACGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCgg                                                                                                                                                                                                          >  2:152323/1‑251 (MQ=255)                                 ggCGATCATCGAGCAAACGTTGTTTGAGCTAGGTG‑ATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACGGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTTCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTGACCCGGTc                                                                                                                                                                                                        <  2:214334/251‑1 (MQ=255)                                     atcatcGAGCAAACGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCCGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCCGcc                                                                                                                                                                                                    >  2:205221/1‑251 (MQ=255)                                          cGAGCAAACGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGCTGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGGTCCTTTCGCGTCAGGTTGGTGTGCTTCGCAAACAGGCGCTGCACCTTAACTTACCCGGTCAGCCGAAGt                                                                                                                                                                                               >  2:35217/1‑251 (MQ=255)                                               aaaCGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAtt                                                                                                                                                                                          >  1:172525/1‑251 (MQ=255)                                                 aCGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAg                                                                                                                                                                                                <  1:152994/243‑1 (MQ=255)                                                 aCGTTGTGTGAGCT‑AGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAg                                                                                                                                                                                                >  2:152994/1‑243 (MQ=255)                                                      gtggtgagttagtggAGGAAATAAGTTGCCACCTGGTGCTCACCACGGGCGGTACTGGCCCGGCGCGTCGTGACGTAACGGCCGATCCGACGCTGTCAGTAGAGGACCGTGAGATTCCTGTCTTTGGTGAACTGATGCGCCTGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGTGTGAGTCGCAAACAGGCGCTGATCCTTAACTTACCCGGGCAGCCCAAGTCGATTAAaga                                                                                                                                                                                     <  2:227523/239‑1 (MQ=255)                                                             gctagtggATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTgg                                                                                                                                                                             >  1:324975/5‑251 (MQ=255)                                                               tagtggCTGAAATGAGTTGCCCTCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGCTGCGACGCTGGCAGTAGCGGACCGCGCGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTAATGTTGGCGTGATTCGCAAACATGCGCTGATCCTTAACTTACACGGTAAGCATCACTCTATTAAAGAAACACTGGa                                                                                                                                                                            >  2:233595/3‑250 (MQ=255)                                                                      tGAAATGCGTTGCCATCTGGTGCTCACCACAGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGCGACGCTGGATGGtgtg                                                                                                                                                                     >  2:59818/1‑250 (MQ=255)                                                                      tGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCCACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAGGGGACGCTGGAAGgtgtt                                                                                                                                                                     >  2:248887/1‑249 (MQ=255)                                                                      tGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAACGAGACGCGGGGAGGtgtg                                                                                                                                                                     >  2:398698/1‑250 (MQ=255)                                                                                         gTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCTCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGt                                                                                                                                                 <  1:35217/251‑1 (MQ=255)                                                                                                ccacGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTa                                                                                                                                          <  1:220651/251‑1 (MQ=255)                                                                                                    gggCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGAGGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACg                                                                                                                                      <  1:28124/251‑1 (MQ=255)                                                                                                     ggCGGAACGGGCCCGGCGCGTCGGGACGTACCGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACgg                                                                                                                                     <  2:142557/251‑1 (MQ=255)                                                                                                            cTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTAttttt                                                                                                                              >  1:31307/1‑251 (MQ=255)                                                                                                              ggCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTAttttt                                                                                                                              <  1:457079/249‑1 (MQ=255)                                                                                                                       cgtcgtGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACg                                                                                                                                      >  1:292621/1‑232 (MQ=255)                                                                                                                       cgtcgtGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACg                                                                                                                                      <  2:292621/232‑1 (MQ=255)                                                                                                                            ggACGTAACGACCGATACGACGCGGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGAGAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTGTGCCAGCGTACCGTACt                                                                                                              <  2:325861/250‑1 (MQ=255)                                                                                                                                  aaCGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCAtt                                                                                                         <  1:251984/250‑1 (MQ=255)                                                                                                                                    cGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAg                                                                                                      <  1:405751/251‑1 (MQ=255)                                                                                                                                      cccGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGaa                                                                                                                                                                   <  1:67376/188‑1 (MQ=255)                                                                                                                                      cccGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGaa                                                                                                                                                                   >  2:67376/1‑188 (MQ=255)                                                                                                                                      cccGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGtt                                                                                                    >  1:429157/1‑251 (MQ=255)                                                                                                                                      cccGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGtt                                                                                                    <  1:164148/251‑1 (MQ=255)                                                                                                                                         gATGCGACGCTGGCAGTAGCGGGCCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTGGCCAGCGTACCGTACTGCATTCAGTTGCt                                                                                                 <  1:347455/251‑1 (MQ=255)                                                                                                                                          aTGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGAGCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTg                                                                                                <  2:429157/251‑1 (MQ=255)                                                                                                                                            gCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGa                                                                                              <  1:421924/251‑1 (MQ=255)                                                                                                                                                    ggCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGcc                                                                                        <  1:205221/249‑1 (MQ=255)                                                                                                                                                    ggCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCa                                                                                       <  1:112648/250‑1 (MQ=255)                                                                                                                                                    ggCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGATGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTATATTAAAGAGACGATGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGATTGAAGGGcact                                                                                      >  2:253128/1‑248 (MQ=255)                                                                                                                                                    ggCAGTAGCGGACCGCGAGATGCCTGGCGTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTTCCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGGCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTTCACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCAt                                                                                      <  2:283357/251‑1 (MQ=255)                                                                                                                                                             ggACCGCGAGATGCATGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGCACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGCGTGAAGGACGCTGCGGGTACCGTTGTGGAACACGGGATTTTTGCCAGCGTACGGTACAGCATTAAGTTGCTGGAGGGGCCAACCGTTGaa                                                                              >  2:321023/1‑250 (MQ=255)                                                                                                                                                              gACCGCGAGATGCCTGGCTTTGGTGAATAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCTATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAAc                                                                            <  2:324975/251‑1 (MQ=255)                                                                                                                                                               aCCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACg                                                                           >  1:216762/1‑251 (MQ=255)                                                                                                                                                                  gcgAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCa                                                                        <  1:153472/251‑1 (MQ=255)                                                                                                                                                                     agaTGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATCCGTTGAAACGGCAcc                                                                      >  2:101407/1‑250 (MQ=255)                                                                                                                                                                                      gAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGACGGGCCATACGTTGAACCGGCACCAGAAGTGGTTGCTGCAtt                                                    >  2:222251/1‑251 (MQ=255)                                                                                                                                                                                       aaCAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTTGTGTGATTCGCAAACAGTCGCTGAACCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTc                                                   <  2:69483/251‑1 (MQ=255)                                                                                                                                                                                            aTGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGACGTGGTTGCTGCATTCAGAcc                                              >  1:435190/1‑251 (MQ=255)                                                                                                                                                                                                       aaGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAGGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCa                                    <  1:248887/249‑1 (MQ=255)                                                                                                                                                                                                        aGCCTGCTTTTTGTACCAACTGCGATCCTTTCGCGGCATGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCGGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTAGTTGTTCCAGCGGACCGTACTGCATTCAGTCGCTGGAAAGGCCAGACGGTGAAACGGCACCAGAAGTGGTTAGTTCATTCAGACCGAAGAGTGCAAg                                  <  2:465618/251‑1 (MQ=255)                                                                                                                                                                                                             gCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAcgc                              <  1:467644/250‑1 (MQ=255)                                                                                                                                                                                                                 tttGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACg                          <  1:384946/250‑1 (MQ=255)                                                                                                                                                                                                                    gTACCAACTGCGATCCTTTCGCGTCCGGTTGGTGGGATTCGCAAACAGGCGCTGATCCTTAAATTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATGCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGTCGTTAg                      <  2:366685/251‑1 (MQ=255)                                                                                                                                                                                                                     tACCAACTGCGATTCTTTCGCGTAAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGGCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTTCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAg                      <  2:268424/250‑1 (MQ=255)                                                                                                                                                                                                                     tACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGACGGACGCTGCGGGGAACGTtg                                                                                                                                               >  2:354305/1‑129 (MQ=255)                                                                                                                                                                                                                     tACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATGCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTTTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTtg                                                                                                                                               <  1:354305/129‑1 (MQ=255)                                                                                                                                                                                                                         aaCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGCAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAAt                 >  1:275410/1‑251 (MQ=255)                                                                                                                                                                                                                              cGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaaaa            <  2:172525/251‑1 (MQ=255)                                                                                                                                                                                                                              cGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGGACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaaaa            <  2:76192/251‑1 (MQ=255)                                                                                                                                                                                                                                  ccTTTCGCGTCAGGGTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATc        <  2:464000/251‑1 (MQ=255)                                                                                                                                                                                                                                      tCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccc     <  1:101407/250‑1 (MQ=255)                                                                                                                                                                                                                                       cgcgTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccccc   <  2:4090/251‑1 (MQ=255)                                                                                                                                                                                                                                          gacAGGTTGGGGTGATTCGCCAAAAGGAGCTGCTCCTTACCTTAAATGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGGTGTACACGGAATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATCCGTTGAACCGGCACCAGAAGTGGTTGCTGCATTCCGACCGGAGAGTGCAATACGCAACGATAGCGAATAAAACAATaccccct  >  2:291488/3‑248 (MQ=255)                                                                                                                                                                                                                                           tCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTa                       <  2:101126/227‑1 (MQ=255)                                                                                                                                                                                                                                           tCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTa                       >  1:101126/1‑227 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               GAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCT‑GGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCG  >  NC_000913/9426‑9906

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