Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,116,656	C→A	R189L (CGT→CTT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,656	0	C	A	98.7%	228.4 / NA	80	R189L (CGT→CTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (36/42);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCAC  >  NC_000913/4116511‑4116798                                                                                                                                                  |                                                                                                                                               aCGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAagag                                                                                                                                             >  2:340255/1‑149 (MQ=255)  cgcgTGCCGCCTTTGGCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAACATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAgaga                                                                                                                                            <  1:359767/149‑1 (MQ=255)  cgcgTCCCGCCTTCGCCGCCAATGTCAGTCTGCGCGTATACTCCGTAACAACGACGCACTTCTGGCCGCATCTCGCCCGTAATATCCAGCACTCCTAGCTTTTTGTCACCCCACGCCACACTATGGATGTTGAACAACATGGTAAgaga                                                                                                                                            <  2:555196/149‑1 (MQ=255)       gccgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGCTGTTGAACAACATGGAAAGAGAGGCGt                                                                                                                                       >  1:620736/1‑149 (MQ=255)         cgccTTTGCCGCCAATGTTAGTCTGACCGTCTATAATGAATTGACTACTCTTTGCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAagagagcagtta                                                                                                                                     <  2:227535/149‑7 (MQ=255)         cgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTg                                                                                                                                     >  1:52264/1‑149 (MQ=255)         cgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTg                                                                                                                                     <  2:291277/149‑1 (MQ=255)         ccgcTTTGCCGCCAAGGTTAGTCTGACCGTATACTGCGGAAGAACGACGCACTTCGGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTg                                                                                                                                     <  2:68205/146‑1 (MQ=255)          gccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGAAGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCACTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTgg                                                                                                                                    >  2:83822/1‑149 (MQ=255)           cccttGCCGCCAAGGTTAGTCTGCCCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGgt                                                                                                                                   <  2:668919/146‑1 (MQ=255)           acTTGGCCGCCACTGTTAGTCTGACCGTATGCCTCGGCCGACCGGAGCACTTCGGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAGGGATGTTGAACAACGTGGCAAGAGAGGGggtggt                                                                                                                                   >  2:772702/2‑149 (MQ=255)            cTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCCCTTCTGCCACCATCTCGCGCGGAATATCCAGCACTTCCAGACTTTTGTCGCCCCAGTCCCGGGCATGGATGTTGCACAACCTGGTACAGGAGGCGTTggcg                                                                                                                                  >  2:393077/1‑147 (MQ=255)                 ccgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATc                                                                                                                             <  1:195304/149‑1 (MQ=255)                    gcAATGTTCGTCTGACCGTATACTTCGGAAGAACGACGCACCTCTGGCAGCATCTCGCTCGGAATATCCAGCACTTCCAGCATTGTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGt                                                                                                                          <  2:662446/148‑1 (MQ=255)                     caaagctAGTCTGACCGTATACTTCCGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCCCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTc                                                                                                                         <  2:309936/143‑1 (MQ=255)                         gTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACAt                                                                                                                     <  1:3364/149‑1 (MQ=255)                         gTTAGTCTGACCGTATACTACGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGACATGTTTGTTGAACAGCATGGTGAGGGAGGCGTTGGTGTGATCGGTCACAt                                                                                                                     >  1:791552/1‑149 (MQ=255)                          ttAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATg                                                                                                                    <  2:613006/149‑1 (MQ=255)                          ttAGTCTGACCGTATACTTCGGAAGAACCACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATg                                                                                                                    <  2:335624/149‑1 (MQ=255)                           tAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGCTGGTGTAATCGGTCCCATgg                                                                                                                   >  2:1025619/1‑149 (MQ=255)                              tCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCAAGCACTTCCAGCATTTTGTAATCCCATTCCAGGGAATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAGTCGGTCACATGCacc                                                                                                                >  1:1072740/1‑148 (MQ=255)                                 gACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTAGGGATGTTCAACCACATGGCAAGAGAGGCGTTGGTGTAATCGGTAACATGGGCGCgg                                                                                                             >  2:335377/1‑149 (MQ=255)                                  aCCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTATCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                                                                                            >  1:35766/1‑149 (MQ=255)                                  aCCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATACAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                                                                                            >  1:579314/1‑149 (MQ=255)                                   ccGTATACTTCGGAAGAACGACGAACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGcc                                                                                                           >  1:574565/1‑149 (MQ=255)                                    cGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGgcc                                                                                                          >  1:454750/1‑149 (MQ=255)                                    cGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCAAGCATTTTTTCGTCCCAGTCCAGGGTAGGGATGTTGAACAACATGGTAAGAGAGGCGGTGGTGTGATCGATCACATGGCCACGGccc                                                                                                          >  2:459969/1‑149 (MQ=255)                                       atccTTCGGAAGAACGACGCCCTTTGGGCAGCATCTCGCGCGGAATATCCAGAACGTGCATAGCTTCGTGGTCCGAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                                                                                                       <  2:389960/146‑1 (MQ=255)                                       ataccgcGGAAGCACGACGCCCTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                                                                                                       <  1:440782/143‑1 (MQ=255)                                       ataCCTCTGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                                                                                                       <  1:837132/149‑1 (MQ=255)                                        taCTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                                                                                                      <  1:491188/149‑1 (MQ=255)                                           ttCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGCTGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACCTGGACACGGCCCCGAGGCa                                                                                                   >  2:312632/1‑149 (MQ=255)                                              ggAAGAACGACGCACTTCGGGCAGGATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                                                                                                <  2:384721/149‑1 (MQ=255)                                                aagaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTc                                                                                              >  1:1004060/1‑149 (MQ=255)                                                 agaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATTGACACGGCCCTGAGTCATTTTcc                                                                                             <  1:574238/149‑1 (MQ=255)                                                  gagcgacgCACTTCTAGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCa                                                                                            <  2:881466/146‑1 (MQ=255)                                                  gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCa                                                                                            <  2:493813/149‑1 (MQ=255)                                                  gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCa                                                                                            <  1:273181/149‑1 (MQ=255)                                                  gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTCCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCa                                                                                            <  1:508197/149‑1 (MQ=255)                                                     cgacgCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGAt                                                                                         >  1:107187/1‑149 (MQ=255)                                                      gacgCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGGTGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGGCACATGGACACGGGCCTGAGGCATGTTCCAGATa                                                                                        >  2:717388/1‑149 (MQ=255)                                                       aggCACATCTGGCAGAATCTCGCAAGGATTATCCCGCCCTACCAGCACTTTGTCGTCACAGTTAAGGGTATGGATGTTGAACAACATTGTAGGAGTGGCTATGGTGGAATCGGTCACATCGACACGGCCTTGAGTCATTTTCCAGAtga                                                                                       <  2:77394/147‑3 (MQ=255)                                                       acgCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCCTTTTCCCGATaa                                                                                       >  2:745636/1‑149 (MQ=255)                                                            cTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGCGTCATTTTCCAGATAAGCCAc                                                                                  >  2:267704/1‑149 (MQ=255)                                                                 ggCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGACGATGTAGTAATCAGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATc                                                                             >  2:280909/1‑149 (MQ=255)                                                                   cagcaCCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTGTAGGGTATGGATGTTGAACAACATGGTAAGTGTGGTGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCaa                                                                           <  1:281596/149‑1 (MQ=255)                                                                    agcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAc                                                                          >  1:561332/1‑149 (MQ=255)                                                                    agcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGCATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAATCCACGTATCCAc                                                                          >  1:591600/1‑149 (MQ=255)                                                                           cgcgcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGAGAAGCTACGTATCAACCTTACCa                                                                   <  2:371031/149‑1 (MQ=255)                                                                            gcgcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCaa                                                                  <  1:340255/149‑1 (MQ=255)                                                                            gcgcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCaa                                                                  >  2:979836/1‑149 (MQ=255)                                                                             cgcgGAATATCCACCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCaaa                                                                 <  1:572845/149‑1 (MQ=255)                                                                              gcgGAATATCCAGCACTTCCAGCATTTTGGCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAGTCCGTCACAAGTGCACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAAc                                                                <  2:123476/149‑1 (MQ=255)                                                                               cgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACa                                                               <  1:380702/149‑1 (MQ=255)                                                                                 gAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGc                                                             >  1:1016473/1‑149 (MQ=255)                                                                                   attgccAGCCCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                                           <  2:574565/145‑1 (MQ=255)                                                                                      tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTc                                                        >  1:190009/1‑149 (MQ=255)                                                                                        cacctcTTCCAGCATGTGGTCGGCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGACTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcac                                                      <  2:981023/144‑1 (MQ=255)                                                                                        cAGCACTTCCAGCATTTTGTCGTACCAGTCAAGAGTATGGACGTTGAACAACAAGGTAAGAGAGGCGTTGGGGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAACCAGCCATacac                                                      >  2:694581/1‑149 (MQ=255)                                                                                          gCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACGTGGTAATAGAGGCGTTGGTGTAATCGGTCACATGTACACGGCCCTGAGTCATTTTTCAAATAAGCCACGTATCAACCGTACCAAACAGCAATTcacca                                                    <  2:724400/149‑1 (MQ=255)                                                                                          gCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACCTGGTAATAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTGCCGATAAGCCACGTATCAACCGTACCAAACAGCAATTcacca                                                    <  2:724370/149‑1 (MQ=255)                                                                                                ccAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGCCGt                                              >  2:364677/1‑149 (MQ=255)                                                                                                 cAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTg                                             >  1:1088895/1‑149 (MQ=255)                                                                                                    cATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTAACATGGAAACGGAATTGCGTCATTTATCATATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCAc                                          <  2:395890/149‑1 (MQ=255)                                                                                                     aTTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCCCg                                         >  2:263682/1‑149 (MQ=255)                                                                                                          gtcgtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgc                                    >  2:839294/1‑149 (MQ=255)                                                                                                              tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                >  1:1037690/1‑149 (MQ=255)                                                                                                              tcCCAGTCCAGGGTATGGAGGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACAGGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                <  2:876700/149‑1 (MQ=255)                                                                                                              acccGGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                <  2:437022/148‑1 (MQ=255)                                                                                                                ccAGTCCAGGGTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGGGc                              >  1:398683/1‑149 (MQ=255)                                                                                                                   gTCCAGGTTAGGCACGTTGAACAACATGGTAAGCGCGGCGTTGGTGTCATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCtt                           <  2:337789/149‑1 (MQ=255)                                                                                                                      cAGGGTATGGATGTTGAACAACATGGTAAGGGAGGCGTTGGTGGAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACCTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCca                        <  2:503327/149‑1 (MQ=255)                                                                                                                        gggTATGGATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCcaca                      <  2:642920/149‑1 (MQ=255)                                                                                                                             tGGATGTTGAACAACATGGTAAGAGAGGCGTTGGGGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTc                 <  2:673828/149‑1 (MQ=255)                                                                                                                               gATGTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACCTGGTCGa               >  2:109848/1‑149 (MQ=255)                                                                                                                                  gTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGa            <  1:312632/149‑1 (MQ=255)                                                                                                                                  gTTGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTAACCACGACGTGGACGCTCGCGAGAGCCTTCCACATGGTCGAGGa            >  1:742542/1‑149 (MQ=255)                                                                                                                                    tGAACAACATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc          <  2:1016473/149‑1 (MQ=255)                                                                                                                                      aacaacATGGTAAGGGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCca        <  2:758829/149‑1 (MQ=255)                                                                                                                                      aacaacATGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCCATTCACCCCGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATccc        >  1:830746/1‑148 (MQ=255)                                                                                                                                       acaacaTGGTAAGAGAGGCGTTGGCGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac       <  2:571429/149‑1 (MQ=255)                                                                                                                                            aTGGTAAGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGAATCAACCGTACCACACAGCAATTCACCACGACGTGCACAACCGACAAAGCCTTCCACATGGTCGAGGGTCcacctcac  >  2:7654/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                               ACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCAC  >  NC_000913/4116511‑4116798

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑