Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,322,707	G→A	intergenic (‑202/+26)	secG ← / ← glmM	Sec translocon subunit SecG/phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,322,707	0	G	A	98.2%	146.1 / NA	55	intergenic (‑202/+26)	secG/glmM	Sec translocon subunit SecG/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/13);  minor base C (0/1);  total (40/15)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATTGCCACAATAAGGAAAACTACTAAAAGAGCTTCATACATAATCAACCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCTGCGGTGTAACGTACGTTAACC  >  NC_000913/3322463‑3322952                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      cAATTGCCACAATAAGGAAAACTACTAAAAGAGCTTCATACATAATCAACCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTAGCTGAGAAAAAAGGCGACAATCGCCACCtttt                                                                                                                                                                                                                                                 >  1:290148/1‑251 (MQ=255)     tGCCACAATAAGGAAAACTACTAAAAGAGCTTCATACATAATCAACCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGc                                                                                                                                                                                                                                             >  1:357233/1‑251 (MQ=255)              aaGGAAAACTACTAAAAGAGCTTCATACATAATCATCCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTGAGCCAGTTATc                                                                                                                                                                                                                                     <  1:307994/250‑1 (MQ=255)                          taAAAGAGCTTCATACATAATCAACCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                                                                                                                                                                                        >  1:267086/1‑251 (MQ=255)                                     catacataATCAACCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTGTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACt                                                                                                                                                                                                              <  1:284383/250‑1 (MQ=255)                                                         cTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                                                                                                                                                                                          >  1:174565/1‑250 (MQ=255)                                                             cGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAtt                                                                                                                                                                                     <  2:257919/251‑1 (MQ=255)                                                              ggAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGAGCGTGAATACTACCCGACGCGAATGACCTTCGCAAGGGCAA‑TTTTATGGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGAGTTGGCTGAGAAAAAAGGCGAAAATCGCAACCTTTTTAGCAAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGCTGCGGAGTGCAAATTc                                                                                                                                                                                    >  1:379332/1‑251 (MQ=255)                                                                 gtttccGCGTACCAATGCTTCAACCAATTAAGCTGGATTTATTAATCTTTACCACTGAAGCGTGTGTGAAAACTAACCAAAGCTAATGACCATCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAACCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCAACTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCCGTGTGCAAATTCAGt                                                                                                                                                                                 <  2:9198/247‑1 (MQ=255)                                                                 gttgcagcgacccAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGt                                                                                                                                                                                 >  1:320143/12‑251 (MQ=255)                                                                         gTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGGATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTgcg                                                                                                                                                                         >  1:331680/1‑251 (MQ=255)                                                                                      aaCCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc                                                                                                                                                            >  1:248539/1‑251 (MQ=255)                                                                                           aTAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTcc                                                                                                                                                        >  2:146509/1‑250 (MQ=255)                                                                                                      gTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCAcgcg                                                                                                                                             >  2:300609/1‑250 (MQ=255)                                                                                                          ttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAAtt                                                                                                                                         <  1:98768/250‑1 (MQ=255)                                                                                                            aaTCTTTCCCACTGACGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAAc                                                                                                                                      >  1:145276/1‑251 (MQ=255)                                                                                                                  tCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCCACTGCGGAAAAAAACAGCCAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCg                                                                                                                                                                                                        >  1:42866/1‑179 (MQ=255)                                                                                                                  tCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCATCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCGTGTTAGCTAGTTATCTAACGCTTTAAACGGCTGTTACTGCATCg                                                                                                                                                                                                        <  2:42866/179‑1 (MQ=255)                                                                                                                             cGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGAttt                                                                                                                     >  1:391120/1‑251 (MQ=255)                                                                                                                               ggTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGTATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAc                                                                                                                   >  2:124783/1‑251 (MQ=255)                                                                                                                                gtgtGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTCACGGTTCGGTGCCGGATTTACg                                                                                                                  >  1:134376/1‑251 (MQ=255)                                                                                                                                  gtgAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgc                                                                                                                 >  1:224718/1‑250 (MQ=255)                                                                                                                                  gtgAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCGTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGTCTTCCACCAGCACGCGAATTAACGGTTCGGTGCCTGATTTACgc                                                                                                                 <  2:389016/250‑1 (MQ=255)                                                                                                                                              cAAAGCGAATGACCTTCGCAAGGGCAA‑GTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAAGTAACGTTTCGGTGCCGGATTTACGCAGCAACACGCGTc                                                                                                    <  2:267086/251‑1 (MQ=255)                                                                                                                                                aagcgacagaactccGCAAGGGCAA‑TTGTATCGCATTGTATCAACTGCGCAAAACAACACCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGACACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCCGCAACACGCGTcc                                                                                                   >  2:383623/15‑250 (MQ=255)                                                                                                                                                   cGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACgg                                                                                               >  2:230285/1‑251 (MQ=255)                                                                                                                                                    gAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGCAAAAACCAGCAAAACCGATGTGTTCGATGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGt                                                                                              >  2:195568/1‑251 (MQ=255)                                                                                                                                                    gAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGTTTTACGCAGCAACACGCGTCCACGGt                                                                                              >  2:123983/1‑251 (MQ=255)                                                                                                                                                     aaTGACCTTCGCAAGGGCCA‑TTTTATCGCATTGTCTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTAGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTCCGCAGCAACACGCGTCCACGGtt                                                                                             >  2:53082/1‑251 (MQ=255)                                                                                                                                                        gACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGCGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGcc                                                                                          >  2:295393/1‑251 (MQ=255)                                                                                                                                                         aCCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAACCAGCAAAACCGGTGTGTTCGCTGAGAAAAAAGGCGACCATCCCCACCTTTTTAGCCCGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCCGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGccc                                                                                         >  1:23292/1‑251 (MQ=255)                                                                                                                                                                cAAGGGCAAATTTTA‑CGCATTGTATCAACTGCGGAAAAAAACAGCAAAACAGATGTGTTAGCTGAGAAAATAGGAGACAATCGCCACCTATATAGCCAGTTCTCTAACGCTATAAACGGCTTTTACTACATCAGCGAAGAGGGGTGCCAGTTCATTCACCTGCTCTTAGTCTTCGCCTTCCAACATCACGCGAATTAACAGTTCGGTGCCGGGTTTACGCAGCAacccacatcaacgggtgccaagcgc                                                                                    >  2:236379/1‑226 (MQ=255)                                                                                                                                                                cAAGGGCAA‑TTTTATCGCATTGTATCACCTGCGGCAAAAAACCGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGACAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCCCGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAgcgc                                                                                    >  2:181825/1‑249 (MQ=255)                                                                                                                                                                            ttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTGCGCTGAGAAAAAAGGCTACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACctc                                                                        <  1:53082/250‑1 (MQ=255)                                                                                                                                                                             tATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACTGTTCCCCAGCGCAGCTTCAACctct                                                                       <  2:19511/250‑1 (MQ=255)                                                                                                                                                                              aTCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACTCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGc                                                                     >  2:5678/1‑251 (MQ=255)                                                                                                                                                                                 gCATTGTATCAACTGCGGAAAAAAACAGCAAAACCTATGTGGTCGCTGAGAAAAAAGGCGACCATCGCCTCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCTCAGCTGCAACCTCTGCGGt                                                                  <  2:250323/251‑1 (MQ=255)                                                                                                                                                                                     tGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACg                                                              >  2:255308/1‑251 (MQ=255)                                                                                                                                                                                      gTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTACCGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACgg                                                             >  1:280186/1‑251 (MQ=255)                                                                                                                                                                                             cTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCCGCTTCAACCTCTGCGGTCACGGCTTTAAc                                                      >  1:296365/1‑251 (MQ=255)                                                                                                                                                                                                   aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTCTCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCCTCAACCTCTGCGGTCACGCCTTTAACTGACTc                                                >  2:157166/1‑251 (MQ=255)                                                                                                                                                                                                           agcaaatacgttgttttCGCCGAGGAAAAAGGCGACAAACGCCCCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCTCGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCaa                                        <  2:81585/236‑1 (MQ=255)                                                                                                                                                                                                              aaaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACATTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTACATCGGCGATACGAAGAACAAAATCAGTCAAATAAGAATAGTCTTCACAATAAAACATCACGAGCATTAACGGTTCAGTGCCGGGTTGGCGCAGCAGCAGGGGGGCGCGGGTGCCCAGCGCAGCTGCAACCTCTGCGGTCACGGcggtaggtgactcaggctgaaggg                                     >  1:246209/1‑228 (MQ=255)                                                                                                                                                                                                                aaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTgg                                    >  2:34154/1‑250 (MQ=255)                                                                                                                                                                                                                     aTGTGTTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATcgcc                              <  2:186588/251‑1 (MQ=255)                                                                                                                                                                                                                       gtgtTCGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATcgccgc                            >  1:402763/1‑251 (MQ=255)                                                                                                                                                                                                                           tCGCTTAGAAAAAAGGCGTCAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCGTCCACCATCACGCGAATGAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTAc                         <  1:213734/250‑1 (MQ=255)                                                                                                                                                                                                                            cGCTGAGAAAAAAGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCt                       <  1:295393/251‑1 (MQ=255)                                                                                                                                                                                                                             gCTGAGAAAAAAGGCGACACTCGCCACCTTTTTCGCCAGTTCTCTAAAGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGAAAATTCAGTCACCTGCGCTTCGTCTTCACCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGACGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGCCTCCTGCTAAAGTGGATCGCCGCTACCTg                      >  2:384361/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaaaaaGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCTGCGgtgt                >  1:410810/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaaaaGGCGACAATCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCTGCGgtgt                >  2:59418/1‑250 (MQ=255)                                                                                                                                                                                                                                      aaaGGCGACAATCGACACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGAGAATTAACGGTTAGGTGCCGGATTTACGCAGCAACAAGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCAAGGCTTTAACTGACTAATGCTCAAGTGGATCGCCGCTACATGCGGGGTAac             >  1:200684/1‑251 (MQ=255)                                                                                                                                                                                                                                         ggCGACAATCGCCAACCCTTTAACCAGTTCTCTAACTCTTTAAACCACTTTTCCCGCACCGCCGAAGCCATGTCCAAACTCAGTCCCCTGCGCTTCGTCTCCCCCTTCCACCTTCACGCGCATTAACGGATCGGTGCCGGATTTACGCAGCAACACGCGTCCCCGGTTGCCTAAAGCAGCTTCAACCTCTGCTGTCACTGCCTTAACTGACTCATGCTAACGTGGACCGCCGTTGCCTGCGGTGCAacgt           >  2:187206/1‑250 (MQ=255)                                                                                                                                                                                                                                                agcgccacCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCGTTTACTGCATCGGCTATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCGGCGGTGTAACGTACGTTAAc   <  1:128561/249‑1 (MQ=255)                                                                                                                                                                                                                                                aTCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCTGCGGTGTACCGTACGTTAAc   >  1:168414/1‑251 (MQ=255)                                                                                                                                                                                                                                                aTCGCCACCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCTGCGGTGTAACGTACGTTAAc   >  1:101738/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cgccacCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCTGCGGTGTAACGTACGTTAAcc  >  1:215743/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CAATTGCCACAATAAGGAAAACTACTAAAAGAGCTTCATACATAATCAACCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAA‑TTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACACGCGTCCACGGTTGCCCAGCGCAGCTTCAACCTCTGCGGTCACGGCTTTAACTGACTCATGCTCAAGTGGATCGCCGCTACCTGCGGTGTAACGTACGTTAACC  >  NC_000913/3322463‑3322952

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