Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,050,472	G→A	L66L (CTG→TTG)	gcvT ←	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,050,472	0	G	A	98.6%	238.4 / NA	70	L66L (CTG→TTG)	gcvT	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/32);  minor base C (1/0);  total (37/33)
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

CAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAGAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATT  >  NC_000913/3050232‑3050716                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     cAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCGGCAGATAACGCAAAAACTCCCgg                                                                                                                                                                                                                                            >  2:609018/1‑251 (MQ=255) cAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCgg                                                                                                                                                                                                                                            >  1:498137/1‑251 (MQ=255)     cATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGc                                                                                                                                                                                                                                        <  1:227319/251‑1 (MQ=255)      aTGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCg                                                                                                                                                                                                                                       <  2:190121/251‑1 (MQ=255)        ggAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGc                                                                                                                                                                                                                                     <  2:225508/251‑1 (MQ=255)                catcaCGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCcgcg                                                                                                                                                                                                                              >  1:409000/1‑250 (MQ=255)                     cGAACGGTAATTTCGCTGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGa                                                                                                                                                                                                                             >  2:108013/1‑246 (MQ=255)                     cGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGa                                                                                                                                                                                                                             >  2:630419/1‑246 (MQ=255)                     cGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGa                                                                                                                                                                                                                             >  1:223421/1‑246 (MQ=255)                      gAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGACAACGCAAAAACTCCCGGGTGCGGCTGCCGCGa                                                                                                                                                                                                                             >  2:121013/1‑245 (MQ=255)                        aCGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGTTAATCCAGGAGAGGTCGTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGAGAACGCAAAAACGCCCGGGTGCGTCTGCCGCGAAGATcgac                                                                                                                                                                                                                     <  2:702876/251‑1 (MQ=255)                              aTTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGt                                                                                                                                                                                                               >  1:204565/1‑251 (MQ=255)                                ttCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGTTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTc                                                                                                                                                                                                              <  1:364042/250‑1 (MQ=255)                                   gATGCCGAAAGGATCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGGGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCatat                                                                                                                                                                                                          >  1:627912/1‑251 (MQ=255)                                    aTGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGGCATAtg                                                                                                                                                                                                         <  1:599377/251‑1 (MQ=255)                                    aTGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGGAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATAtg                                                                                                                                                                                                         >  2:229272/1‑251 (MQ=255)                                       ccGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCTGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGa                                                                                                                                                                                                      <  1:22400/251‑1 (MQ=255)                                                          gggTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGGAGACGATGAGGTCATCTTTCACACCGCCAGAGGCATGCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCGTCGCCACATCGTTCGCCTGCAGATAACGCAAAAACCCCCGGGTGCGGCTGCCGCGAAGATCGACGAGGGTCATATGTGACACATCAAACATTCCGGCa                                                                                                                                                                                   <  2:276322/251‑1 (MQ=255)                                                            gTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTCACAACGCGGCGGAAGAAATCTTCAGTAAAGTAGGAGACGATGAGGTCATCTCTCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAACCTCCAGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACCGCCAACATTCCGGCATc                                                                                                                                                                                 >  2:357418/1‑251 (MQ=255)                                                              aaTCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCgg                                                                                                                                                                               >  2:85737/1‑251 (MQ=255)                                                               aTCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGgt                                                                                                                                                                              >  1:562228/1‑251 (MQ=255)                                                                  cAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGgtacg                                                                                                                                                                           <  1:253769/251‑1 (MQ=255)                                                                              ttttCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCatga                                                                                                                                                                 <  1:632549/249‑1 (MQ=255)                                                                               tttCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGt                                                                                                                                                              >  1:209982/1‑251 (MQ=255)                                                                                ttCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGtt                                                                                                                                                             >  1:416104/1‑251 (MQ=255)                                                                                 tCGCGGGTGGCGGAGTTAACAACGCGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCAAACAGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTtc                                                                                                                                                            >  2:722041/1‑251 (MQ=255)                                                                                  cgcgGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTtcg                                                                                                                                                           <  1:131216/251‑1 (MQ=255)                                                                                     gggTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTtcgtcg                                                                                                                                                        <  1:85737/251‑1 (MQ=255)                                                                                         ggCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGAttt                                                                                                                                                    <  1:57842/251‑1 (MQ=255)                                                                                          gCGGGGTTAACGACGAGGCGGAAGAAATCGTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGAtttt                                                                                                                                                   <  2:168803/251‑2 (MQ=255)                                                                                          gCGGAGTTAACAACGAGGAGGAAGAAATCTTCAGTAAATTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTGCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCAGGATGGTCGGCGATTTg                                                                                                                                                   <  2:685989/251‑1 (MQ=255)                                                                                            ggAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCg                                                                                                                                                 >  1:675999/1‑251 (MQ=255)                                                                                             gAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCg                                                                                                                                                 >  1:325520/1‑250 (MQ=255)                                                                                              aGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGa                                                                                                                                                >  2:147954/1‑250 (MQ=255)                                                                                                    catcGAGGCGGAATAAATCTTCAGTAAAGTAGTAGACGATGATGTCATCTATCACCCCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACACATAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGGTCGTCGATTTGCGAACCGTaa                                                                                                                                         <  1:573971/248‑1 (MQ=255)                                                                                                        gcggCGGAAGCAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTACCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGc                                                                                                                                      >  1:679544/3‑250 (MQ=255)                                                                                                        gAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGc                                                                                                                                      >  1:679554/1‑250 (MQ=255)                                                                                                         aGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAg                                                                                                                                    <  1:147954/251‑1 (MQ=255)                                                                                                          ggCGGAAGCAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGTTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCACAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGCCACATCAAACATTCCGGCATCGGGACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGc                                                                                                                                   >  2:678344/1‑251 (MQ=255)                                                                                                               aagaaATCTTCAGTAAAGTAGAAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGcat                                                                                                                              >  2:111325/1‑251 (MQ=255)                                                                                                                      cTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAACCATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATcc                                                                                                                         >  2:518191/1‑249 (MQ=255)                                                                                                                       ttCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGc                                                                                                                      >  2:110815/1‑251 (MQ=255)                                                                                                                        tCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAg                                                                                                                       <  1:630419/249‑1 (MQ=255)                                                                                                                            tAAAGCAGTAGACGATAAGGACATCTATAACACCGCCAGAGGCATTCACCAAACCAGAGTAAACGGCTTTGACGCTTTTCGCGAGCTTAGCCACACCATTCGCCACCAGATAACGCAAAAACTAAAGGGTGCGGCTGCCGCGACGATCGACGATGCTCATATGTGTCACATCACACATTCCGGAATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTgg                                                                                                                 >  2:365200/1‑251 (MQ=255)                                                                                                                              aaggagtagACGCTGAGGTAGTCTATCTCACCGCCAGAGGCAGTCAACATCCCCGAGTAAAGGGCTTTGCCGCTGGTGGTGATCTCCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGGCACATCAAACAGTCCGTCATCGGTACGTACCGCATGATGTGCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGaa                                                                                                               <  1:454648/247‑1 (MQ=255)                                                                                                                               aggagtagACGAAGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGaaa                                                                                                              <  2:498137/248‑1 (MQ=255)                                                                                                                                         aaTGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAACTCCACCATGc                                                                                                    >  2:118190/2‑251 (MQ=255)                                                                                                                                            gAGGTCATCTACCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTGGCGAACCGTAATGCAGCGGCAACATCCAGCCGGGGAAATCCACGATGCGAg                                                                                                 <  2:325520/251‑1 (MQ=255)                                                                                                                                                     tatcacaactcaagagGCATACAACATCTCCGAGTCAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGACAACGCAAAAACTCCCGGGTCCGCCTGCCGCGAAGATCGACGATGGTCCTATGAGAAACATCACACCTTCCGGCCTCGGTACGTACCGCATGATGTTCGTCGTATTGCGAACCGTCATGCAGCGGCATCATCCAGCCGTGGCAATCAACAATGCCAgccaagaca                                                                                         >  2:716534/13‑243 (MQ=255)                                                                                                                                                      aTCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGTGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATc                                                                                                                                                                                 <  2:66481/161‑1 (MQ=255)                                                                                                                                                      aTCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATc                                                                                                                                                                                 >  1:66481/1‑161 (MQ=255)                                                                                                                                                        cacaCCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCg                                                                                     >  2:730209/1‑251 (MQ=255)                                                                                                                                                           acCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGCGcgccg                                                                                            >  2:342618/1‑241 (MQ=255)                                                                                                                                                           acCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCAGCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGcgccg                                                                                            <  1:342618/241‑1 (MQ=255)                                                                                                                                                             cGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCgtgtgt                                                                                >  1:626952/1‑251 (MQ=255)                                                                                                                                                                   agGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCg                                                                          >  1:566485/1‑251 (MQ=255)                                                                                                                                                                                    gAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCtgttg                                                         <  1:156831/251‑1 (MQ=255)                                                                                                                                                                                          aGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGGGTTGTTCGTACAAAGGAGTCTGTTGTGCCAt                                                   <  1:538881/251‑1 (MQ=255)                                                                                                                                                                                                      cTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCAtt                                       >  2:454343/1‑251 (MQ=255)                                                                                                                                                                                                      cTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCAtt                                       >  1:211223/1‑251 (MQ=255)                                                                                                                                                                                                           ggTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCGCCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAAt                                   <  1:661965/250‑1 (MQ=255)                                                                                                                                                                                                            gTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATGCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATGGAATaa                                 <  2:596454/251‑1 (MQ=255)                                                                                                                                                                                                                    cGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCt                         <  1:518191/251‑1 (MQ=255)                                                                                                                                                                                                                      ccACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGa                       <  1:99033/251‑1 (MQ=255)                                                                                                                                                                                                                        acaTCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACa                     <  1:117263/251‑1 (MQ=255)                                                                                                                                                                                                                                cGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACttttt              <  1:229272/250‑1 (MQ=255)                                                                                                                                                                                                                                      cagATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTg       <  2:566485/251‑1 (MQ=255)                                                                                                                                                                                                                                       agATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTg       >  1:312496/1‑250 (MQ=255)                                                                                                                                                                                                                                         aTAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGaa     >  2:573811/1‑250 (MQ=255)                                                                                                                                                                                                                                           aaCGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAAtt  <  1:494145/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     CAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAGAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATT  >  NC_000913/3050232‑3050716

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