Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,905,761	G→A	G25D (GGT→GAT)	mntP →	Mn(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,905,761	0	G	A	97.8%	143.9 / NA	46	G25D (GGT→GAT)	mntP	Mn(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/24);  minor base C (1/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CTCGTTGCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTCT  >  NC_000913/1905515‑1906003                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   cTCGTTGAAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGcc                                                                                                                                                                                                                                                >  1:134616/1‑251 (MQ=255)             cGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATaaa                                                                                                                                                                                                                                    <  1:348624/251‑1 (MQ=255)              gTGGCACGTACTGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCATGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTATGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAACGATGCCACCCTCCATAAAc                                                                                                                                                                                                                                   >  2:149057/1‑251 (MQ=255)                      tACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATTTCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAAttt                                                                                                                                                                                                                           >  1:255971/1‑251 (MQ=255)                                             ggCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCt                                                                                                                                                                                                    <  2:343377/251‑1 (MQ=255)                                              gCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCtt                                                                                                                                                                                                   <  2:255971/251‑1 (MQ=255)                                                gagaCCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAt                                                                                                                                                                                                 >  1:33165/1‑251 (MQ=255)                                                                    cGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTTTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACGGTTCTTCTTGCGTTTGGTATGTCGATGGATGCAGTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAAcc                                                                                                                                                                              <  2:436004/250‑1 (MQ=255)                                                                        aTGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGAc                                                                                                                                                                         >  1:73129/1‑251 (MQ=255)                                                                             tGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTCTGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGCCATTGCCATGACTATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGACTCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGCGCCGTCGAAACCCTGAcgccgc                                                                                                                                                                    >  1:280257/1‑251 (MQ=255)                                                                                aTGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTGCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGa                                                                                                                                                                 <  2:151256/251‑1 (MQ=255)                                                                                      gcgcacaccatccgggCTTGGTGTTCCTGTATGCCTGTTTGCTGGGG‑CGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCTTGACTATCCCTGCTAATGTTCTTCTTGCGTTTGGTATGTCGAGGGATGCATTTGCTGCATCAATCGGTAAAGCTGCCACCCTCCCTACACCGAAATTTTCTGCAGCATTGCGAACCGGCCTTCTTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTg                                                                                                                                                          >  2:65271/15‑251 (MQ=255)                                                                                      gcgAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                                                                                                                                                                                                         >  2:153299/1‑173 (MQ=255)                                                                                      gcgAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTcc                                                                                                                                                                                                                                         <  1:153299/173‑1 (MQ=255)                                                                                          aaaCCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCCTCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTgggg                                                                                                                                                       >  2:191890/1‑251 (MQ=255)                                                                                           aaCCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGCGGCAGATGATGTGTTTTATGGATACCCCGGTCAGGACACTGTCATGAATCTCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTgggg                                                                                                                                                       >  1:313312/1‑250 (MQ=255)                                                                                                   cGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCAt                                                                                                                                              <  2:228290/251‑1 (MQ=255)                                                                                                       ttgttgTTCATTTATGCAGGTTTGCTGGGGGCGATGATGGGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATAGCACGGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTTCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCGTATTGGTGGTGCCGTCGAAACCCTGACGCCTCTGATCGGCTGGGGCATGGGCATGtg                                                                                                                                           <  2:61585/250‑2 (MQ=255)                                                                                                        tgttgtTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCTAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGtt                                                                                                                                           <  1:191890/249‑1 (MQ=255)                                                                                                                              gCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAAt                                                                                                                    <  1:253468/250‑1 (MQ=255)                                                                                                                              gCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATg                                                                                                                   <  2:141699/251‑1 (MQ=255)                                                                                                                                tGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATgg                                                                                                                  >  1:187029/1‑250 (MQ=255)                                                                                                                                      ggatgatGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTCCTGTTCTTCGTGCGTTTGGTATGTCGAGGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCGGATCGGCTGGGGAATTGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACt                                                                                                           <  2:282227/250‑1 (MQ=255)                                                                                                                                           atgTGTTTTATGGATACTCCGGTCAGGACCTTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGCTGCATTTGCTGACTCAATCGGTACAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGCATGGCACCACTGGAtt                                                                                                      >  2:218307/1‑251 (MQ=255)                                                                                                                                                    aTGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTtgtg                                                                                             >  1:69443/1‑251 (MQ=255)                                                                                                                                                                   aGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTc                                                                              >  2:14196/1‑251 (MQ=255)                                                                                                                                                                    ggACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcg                                                                             <  1:226089/251‑1 (MQ=255)                                                                                                                                                                       cATTGTCATGAATATCACTGCTCCTGTTCTTCTTGCGGTTGGTACTTCGATGGATGCATTTGCTGCATCAAGCGGTAAAGATGCCGCCCTCCATAAACCGAAACTTTCTGAAGCAGTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAGCCCATACGCCGCGGATCGGCAGTGGAATGGGCATGTTagcc                                                                                                                                       <  2:431819/190‑1 (MQ=255)                                                                                                                                                                       cATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGCTTGGTATGTCGCTGGATGCATTTGCTGCATCAATAGGTAAAGATGCCACCCTCCATCAACAGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTagcc                                                                                                                                       >  1:431819/1‑190 (MQ=255)                                                                                                                                                                        aTTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAt                                                                                                                                                                                                 >  1:448030/1‑131 (MQ=255)                                                                                                                                                                        aTTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTGGGTATGTCGATTTATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTAt                                                                                                                                                                                                 <  2:448030/131‑1 (MQ=255)                                                                                                                                                                           gTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCggcgggc                                                                       <  2:379975/250‑1 (MQ=255)                                                                                                                                                                                 aaTATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGt                                                                                                                                                                                          <  1:132841/129‑1 (MQ=255)                                                                                                                                                                                 aaTATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGt                                                                                                                                                                                          >  2:132841/1‑129 (MQ=255)                                                                                                                                                                                      cACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGattatt                                                            <  2:191672/250‑1 (MQ=255)                                                                                                                                                                                                tttcgtcttGCGTTTGGTATGTCGATGGATGCATTGGCTGCATCAATCGGTAAAGATGCCACCCTCCACAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGGGCCGTCGAAACCCTGCCGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGGCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATCATTATTGAGGGtttt                                                   <  2:280257/244‑1 (MQ=255)                                                                                                                                                                                                gTTCTTCTTTCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTc                                                  <  1:266114/250‑1 (MQ=255)                                                                                                                                                                                                gTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTc                                                  <  1:344637/250‑1 (MQ=255)                                                                                                                                                                                                     tcttGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcgg                                                                            >  2:196656/1‑219 (MQ=255)                                                                                                                                                                                                     tcttGCGTTTGGTATGTCGATGGATGCATGTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGGGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTcgg                                                                            <  1:196656/219‑1 (MQ=255)                                                                                                                                                                                                        tGCGTTTTGTATGTCGATGGATGGATTTGCTGCATCAATCGGTAAAGAAGGCCCCCTCCATAAACCGAAATTTTCTGAAGCCTTGGCAACCGGGCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTAGTGAGGGTTTTCGTGgcg                                            <  1:149057/248‑1 (MQ=255)                                                                                                                                                                                                            tttGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTTTTGAGGGGTTTCGTGGCGCAgatgg                                     >  1:305960/1‑250 (MQ=255)                                                                                                                                                                                                                  aTGTCGAGGGATGCATTTGCTGCATCAATCGTTAAAGATGCCACCCTCCATAAACCGAAAGTTTCTGAAGCATTGCGAACCGGCCTTATGTTTGGTGCCGTCGAAACCCTGACGCCGCTGTTCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCGTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTGTTCTTGGCGCAGATgaggaaga                               <  2:419459/251‑1 (MQ=255)                                                                                                                                                                                                                    gTCGATGGATGCATTTGCTGCATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGCACCCCTGGATTGCGTTTGTGCTGCTGATATTCATCGGCGGGGGCATGGTTATTTAGGGGTTTTTTGGCGCacaagttgtcgg                               >  2:393360/1‑238 (MQ=255)                                                                                                                                                                                                                                     ctgcATCAATCGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGGATGATTATTGAGGGGTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCgc            >  1:361123/1‑250 (MQ=255)                                                                                                                                                                                                                                               cGGTAAAGATGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGACTGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGGATGGTTTTTGTGGGGTTTCGTGGCGGAGGTGATGGAGGTGAAGAGCCGCGCCggcggccagggttct  >  2:180621/1‑237 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CTCGTTGCAAAACGTGGCACGTACGGACTGAATACTTTCAGTCAGGCGAGACCATATGCACATCAATCGCTATGCCTGCATGAACGCGAAAACCATCACGGCTTGTTGTTCATGTATGCATGTTTGCTGGGGGCGATGATGTGTTTTATGGATACCCCGGTCAGGACATTGTCATGAATATCACTGCTACTGTTCTTCTTGCGTTTGGTATGTCGATGGATGCATTTGCTGCATCAATCGGTAAAGGTGCCACCCTCCATAAACCGAAATTTTCTGAAGCATTGCGAACCGGCCTTATTTTTGGTGCCGTCGAAACCCTGACGCCGCTGATCGGCTGGGGAATGGGCATGTTAGCCAGCCGGTTTGTCCTTGAATGGAACCACTGGATTGCGTTTGTGCTGCTGATATTCCTCGGCGGGCGAATGATTATTGAGGGTTTTCGTGGCGCAGATGATGAAGATGAAGAGCCGCGCCGTCGACACGGTTTCT  >  NC_000913/1905515‑1906003

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