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.3%	203.8 / NA	60	L66L (CTG→TTG)	gcvT	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (28/30);  minor base C (1/0);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAGAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTAT  >  NC_000913/3050231‑3050718                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       gCAATAATAGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAACGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCACCGAGTAAATGGCTTTGCCGCTTTTGGTGAGCTTCCCCACATCGTTCGCCAGCAGCAAACGCAAAAACTCCCg                                                                                                                                                                                                                                               >  2:336805/1‑251 (MQ=255)      cATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGc                                                                                                                                                                                                                                          <  1:96429/251‑1 (MQ=255)      cATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGc                                                                                                                                                                                                                                          <  1:320317/251‑1 (MQ=255)            aaGGTCATCACGCACGGTAGTTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCACAACCTCCCGGGTGCGGCTGc                                                                                                                                                                                                                                    >  2:80277/1‑251 (MQ=255)                     ccGAACGGTAATTTCGATGCCGAAAGGTTAAGCGTGTTGGTTAAGCCAGGAGAGGTCTTTTTCGCGGGTTGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGAGGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGTGTTTGGTGAGCTTCGCCACATTTTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGag                                                                                                                                                                                                                           <  2:572219/250‑2 (MQ=255)                        aaCGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGCGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATcg                                                                                                                                                                                                                         >  1:346568/1‑250 (MQ=255)                         aCGGTAATTGCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATAACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATcgac                                                                                                                                                                                                                       <  1:362892/251‑1 (MQ=255)                                    gATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCatat                                                                                                                                                                                                            >  1:59423/1‑251 (MQ=255)                                              ggTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATc                                                                                                                                                                                                  >  1:160170/1‑251 (MQ=255)                                                   aGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGCTGGTAATATGTGACACATCAAAc                                                                                                                                                                                              >  2:406398/1‑250 (MQ=255)                                                           gggTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGc                                                                                                                                                                                      >  1:185745/1‑250 (MQ=255)                                                            ggTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGCCGATGAGGTCATCTATCCCACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATTTGACACATCAAACATTCCGGCAt                                                                                                                                                                                    <  2:250510/251‑1 (MQ=255)                                                            ggTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAACCATTCCGGCAt                                                                                                                                                                                    >  1:361785/1‑251 (MQ=255)                                                             gTAATCCAGGAGAGGTCTTGTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGGAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACAGCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATc                                                                                                                                                                                   <  1:464098/251‑1 (MQ=255)                                                              tAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCg                                                                                                                                                                                  >  1:530508/1‑251 (MQ=255)                                                                aTCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGgt                                                                                                                                                                                >  1:130434/1‑251 (MQ=255)                                                                  ccAGGAGAGTTCTTTTTCGCGGGGGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTGGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGgtac                                                                                                                                                                              <  2:464427/251‑1 (MQ=255)                                                                  ccAGGAGAGGTTTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCTGCATCGgtac                                                                                                                                                                              <  2:378670/251‑1 (MQ=255)                                                                                   cgcgGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTtcg                                                                                                                                                             >  1:100063/1‑251 (MQ=255)                                                                                    gcgGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTtcgt                                                                                                                                                            <  1:516175/251‑1 (MQ=255)                                                                                                      aacGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAAt                                                                                                                                          >  2:211782/1‑251 (MQ=255)                                                                                                      aacGAGGCGGAAGAAATCTTAAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCCGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGGGAACCGTAAt                                                                                                                                          >  2:328636/1‑251 (MQ=255)                                                                                                       acGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATg                                                                                                                                         >  1:573302/1‑251 (MQ=255)                                                                                                         gAGGCGGAAGAAATCTTCAGTAAAGGAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCa                                                                                                                                       <  1:567653/251‑1 (MQ=255)                                                                                                                aagaaATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGAATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGAAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGCACCGTAATGCACCGGcat                                                                                                                                >  1:340487/1‑251 (MQ=255)                                                                                                                 agaaATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGcatc                                                                                                                               <  1:427000/251‑1 (MQ=255)                                                                                                                       cTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATcc                                                                                                                           >  1:424086/1‑249 (MQ=255)                                                                                                                             aaaaGGAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACAGCCCCGAGTAAAGGGCTTTGCCGCTTTGGGTGAGCTTCGCCACAGCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGAGCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGttg                                                                                                                   <  2:353712/250‑3 (MQ=255)                                                                                                                                    gtagACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATcc                                                                                                             >  2:66743/1‑250 (MQ=255)                                                                                                                                         cGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCAtc                                                                                                       >  2:62811/1‑250 (MQ=255)                                                                                                                                          gATGTGGTCATCTATCTCCCCGCCAGAGGCATTCAACACCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGc                                                                                                      <  2:480251/251‑1 (MQ=255)                                                                                                                                           aTGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGCTGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATTCACCAAGCg                                                                                                     >  2:149887/1‑251 (MQ=255)                                                                                                                                             gAGGTCATCTATCACACCGCCAGAGGCATTCAACACCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGa                                                                                                    <  1:437657/250‑1 (MQ=255)                                                                                                                                              gggccacctatccccccGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGGGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAgt                                                                                                  <  2:496004/236‑2 (MQ=255)                                                                                                                                                         caccccGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCg                                                                                       <  1:310433/247‑1 (MQ=255)                                                                                                                                                               gCCAGGGGCATTCAACATCCCCGAGTAAAGGTCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGtgtt                                                                                 <  2:504949/251‑1 (MQ=255)                                                                                                                                                                 cAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCAGCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGtgttgt                                                                               <  1:328636/251‑1 (MQ=255)                                                                                                                                                                   gatgcagccaccacccccGAGTTAAGTGATTAGCCGCTTGTGGTGAGCTTCCCCACATCGTTGGCCAGCAGAAAACGCAAAAACTCCCGGGGGCGGCTGCCGCGAAGATCGACGATGGTTATTTGTGACACAGCAAACATTCCGGCACCGGTACGTACCGCATGATGTTCGACGAGTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAAACCTCCATGCGCGCGCCGCAAAGCGTGTGTTGTTc                                                                             <  2:515878/237‑1 (MQ=255)                                                                                                                                                                   gaggcaacaaccATCCCCGATTAAAGGGCTTTGCCGCTTTTGGTGAGCTGCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTc                                                                             <  2:515880/240‑1 (MQ=255)                                                                                                                                                                            aaCATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAgg                                                                    <  2:160170/251‑1 (MQ=255)                                                                                                                                                                              cATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAg                                                                  <  1:406398/251‑1 (MQ=255)                                                                                                                                                                                 ccccGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAGAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCt                                                               <  1:312950/251‑1 (MQ=255)                                                                                                                                                                                      aGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCtgttgt                                                          >  1:257833/1‑251 (MQ=255)                                                                                                                                                                                            gggCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATc                                                    >  1:286368/1‑251 (MQ=255)                                                                                                                                                                                                 ttGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTc                                               >  1:494251/1‑251 (MQ=255)                                                                                                                                                                                                    ccGCTTTTGGGGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTc                                            <  2:494251/251‑1 (MQ=255)                                                                                                                                                                                                     cGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGCGTGCGGCTGCCGCGAAGATCGACGATGGACATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTc                                            >  1:47179/1‑250 (MQ=255)                                                                                                                                                                                                          ttGGTGAGCTTCTCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATTGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATTATGTTCGTCGATTTGCGAACCGTAATGCAGCTGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGa                                       <  1:80335/250‑1 (MQ=255)                                                                                                                                                                                                            ggggAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAAt                                     <  2:100240/247‑1 (MQ=255)                                                                                                                                                                                                                aGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCg                                <  1:342496/251‑1 (MQ=255)                                                                                                                                                                                                                 gCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCgg                               >  2:144977/1‑251 (MQ=255)                                                                                                                                                                                                                  cTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCa                                                                                                                          <  1:303834/159‑1 (MQ=255)                                                                                                                                                                                                                  cTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCa                                                                                                                          >  2:303834/1‑159 (MQ=255)                                                                                                                                                                                                                     cGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCt                           <  2:257833/251‑1 (MQ=255)                                                                                                                                                                                                                      gCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTg                          >  1:437459/1‑251 (MQ=255)                                                                                                                                                                                                                        cacaTCGTTGGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGat                        <  2:130434/251‑2 (MQ=255)                                                                                                                                                                                                                                 cGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTc               >  2:439105/1‑251 (MQ=255)                                                                                                                                                                                                                                     agcagATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATgg           >  2:287493/1‑251 (MQ=255)                                                                                                                                                                                                                                             aCGCAAAAACTCCCGGGGGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGATCCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTACGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTa   <  2:47179/251‑1 (MQ=255)                                                                                                                                                                                                                                              cGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGAAAAAGGTAAAAAGAGACACAACAAACAGGCAGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGAAACCCTGTCCTCATTGAATAAGCGGGGCTGACAACATTTTCATGGTGAAATTAt  >  1:321524/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       GCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCATCTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAGAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTAT  >  NC_000913/3050231‑3050718

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