Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,705	T→G	68.8%	K173T (AAA→ACA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,705	0	T	G	68.8%	176.0 / 92.1	99	K173T (AAA→ACA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (18/13);  new base G (36/32);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATA  >  NC_000913/4116561‑4116851                                                                                                                                                 |                                                                                                                                                   gaaCGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCa                                                                                                                                                <  1:686322/149‑1 (MQ=255)   acgacgCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGa                                                                                                                                              <  3:325069/149‑1 (MQ=255)        gCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGc                                                                                                                                         >  2:21918/1‑149 (MQ=255)        gCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGc                                                                                                                                         >  2:71080/1‑149 (MQ=255)          aCTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCa                                                                                                                                       >  2:217852/1‑149 (MQ=255)            ttCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCAATCTGTCGACCCAGTCCAGGGTATGGATCTTGAACAACATGGTACGAGAGGTGTTGGTGTAATCGGTCACATGCACATGGCCCTGAGTCATTGTCCGGATAAGCcgcg                                                                                                                                     >  4:544432/1‑149 (MQ=255)               tGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTAt                                                                                                                                  >  1:95379/1‑149 (MQ=255)                ggCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATc                                                                                                                                 >  4:654019/1‑149 (MQ=255)                 gcagcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCa                                                                                                                                >  3:559317/1‑149 (MQ=255)                 gcagcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCa                                                                                                                                >  3:544725/1‑149 (MQ=255)                 gcagcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCa                                                                                                                                >  1:489202/1‑149 (MQ=255)                   agcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAAc                                                                                                                              <  4:753816/149‑1 (MQ=255)                    gcaTCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                                                                                                             >  2:590629/1‑149 (MQ=255)                       tctcGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTa                                                                                                                          <  4:288845/149‑1 (MQ=255)                         tcGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTAcc                                                                                                                        <  2:550034/149‑1 (MQ=255)                           gcgcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAg                                                                                                                                               <  2:717395/124‑1 (MQ=255)                           gcgcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAg                                                                                                                                               >  1:717395/1‑124 (MQ=255)                             gcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAAc                                                                                                                    >  1:754169/1‑149 (MQ=255)                             gcgGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAAc                                                                                                                    >  1:40510/1‑149 (MQ=255)                                  atatCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGTGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                                                                                               <  1:707956/149‑1 (MQ=255)                                     tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGAt                                                                                                                                             >  2:349974/1‑116 (MQ=255)                                     tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGAt                                                                                                                                             <  1:349974/116‑1 (MQ=255)                                     tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTc                                                                                                            >  2:7419/1‑149 (MQ=255)                                     tCCAGCACTTCCAGCATTTTGTCGTCCCAGGCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACC‑TACCAAACAGCAATTca                                                                                                           >  1:659249/1‑149 (MQ=255)                                       cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcac                                                                                                          >  1:530469/1‑149 (MQ=255)                                           aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacg                                                                                                      >  3:216122/1‑149 (MQ=255)                                                 aGCATTTTGTTGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCAAGACGTGc                                                                                                <  2:480679/149‑1 (MQ=255)                                                  gCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCa                                                                                               >  1:102288/1‑149 (MQ=255)                                                     ttttGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGc                                                                                            >  2:549823/1‑149 (MQ=255)                                                     ttttGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGc                                                                                            >  4:502645/1‑149 (MQ=255)                                                         gtcgtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgc                                                                                        <  2:40510/149‑1 (MQ=255)                                                          tcgtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgcg                                                                                       >  1:819623/1‑149 (MQ=255)                                                            gtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgag                                                                                     >  3:278244/1‑149 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                    >  2:516477/1‑149 (MQ=255)                                                             tcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCAATCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                    <  4:9120/149‑1 (MQ=255)                                                              cccAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGagag                                                                                   >  2:315256/1‑149 (MQ=255)                                                                cAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGcc                                                                                 >  2:601790/1‑149 (MQ=255)                                                                cAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGcc                                                                                 <  2:14652/149‑1 (MQ=255)                                                                       gggTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCcaca                                                                          >  2:431727/1‑149 (MQ=255)                                                                             ggATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCg                                                                    >  3:648671/1‑149 (MQ=255)                                                                             ggATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCg                                                                    >  4:443796/1‑149 (MQ=255)                                                                              caTGTTGTTCAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                    >  2:353313/2‑132 (MQ=255)                                                                              caTGTTGTTCAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCgaga                                                                                    <  1:353313/131‑1 (MQ=255)                                                                               aTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                                  >  3:315460/1‑149 (MQ=255)                                                                               aTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                                  <  4:757495/149‑1 (MQ=255)                                                                               aTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                                  <  2:777198/149‑1 (MQ=255)                                                                               aTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                                  <  2:102288/149‑1 (MQ=255)                                                                                tGTTGAACAACATGGTACGAGTGGCGTTGGTGTAATCGGACACATGGACACGGCGCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCATACTGCAATTCACAACGACGTGCACGCTCGCGAGACCCTTCCACATGGCCGAga                                                                 >  4:154917/1‑148 (MQ=255)                                                                                 gTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTc                                                                     >  2:629651/1‑144 (MQ=255)                                                                                 gTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTc                                                                     <  1:629651/144‑1 (MQ=255)                                                                                 gTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCaaa                                                                                                                     >  3:844135/1‑96 (MQ=255)                                                                                 gTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCaaa                                                                                                                     <  4:844135/96‑1 (MQ=255)                                                                                 gTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGa                                                                <  4:802655/149‑1 (MQ=255)                                                                                   tGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                              <  4:315460/149‑1 (MQ=255)                                                                                   tGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                              <  4:544725/149‑1 (MQ=255)                                                                                      acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                           <  1:43909/149‑1 (MQ=255)                                                                                      acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGACGAGGATCcac                                                           <  1:467971/149‑1 (MQ=255)                                                                                        aacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcactt                                                         <  2:866968/149‑1 (MQ=255)                                                                                         acaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCca                                                                            >  1:206956/1‑129 (MQ=255)                                                                                         acaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCca                                                                            <  2:206956/129‑1 (MQ=255)                                                                                           aTGGTACGAGAGGCGTTGGTGTAATCGGTGATATGGACACGACCCTGAGTCATTGTCCAGATGAGCCACGAATCAACCGTACCAAAGAGCAATTCACCACGACGTACACACTCGCGATAGCCTTCCACATGGTCGAGGATCcacttcac                                                      <  3:544432/149‑1 (MQ=255)                                                                                           aTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcac                                                      <  1:376708/149‑1 (MQ=255)                                                                                           aTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACAGCCCTGACTCATTGTCCAGATAAGCCACGTATCAACCGTACCCAACAGCAATTCACCACGACGTGCACGCTCGCGAGCGCCTACCACATGGTCGAGGATCcacttcac                                                      >  2:15803/1‑149 (MQ=255)                                                                                              gTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACttt                                                   <  3:593950/149‑1 (MQ=255)                                                                                                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTgg                                                 <  3:560700/149‑1 (MQ=255)                                                                                                 cGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGt                                                <  4:689902/149‑1 (MQ=255)                                                                                                    gagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGcc                                             <  4:150833/149‑1 (MQ=255)                                                                                                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCa                                            <  4:216122/149‑1 (MQ=255)                                                                                                          ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATACACTTCACTTTGGTGCCAGaaaa                                       <  2:371905/149‑1 (MQ=255)                                                                                                              tgtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                           >  2:423337/1‑125 (MQ=255)                                                                                                              tgtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                           <  1:423337/125‑1 (MQ=255)                                                                                                               gtAATCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACg                                  <  3:164750/149‑1 (MQ=255)                                                                                                                tAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                                 <  2:469453/149‑1 (MQ=255)                                                                                                                tAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                                 <  3:584898/149‑1 (MQ=255)                                                                                                                tAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                                 <  2:183468/149‑1 (MQ=255)                                                                                                                 aaTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACggg                                <  4:648671/149‑1 (MQ=255)                                                                                                                  aTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGt                               >  1:405360/1‑149 (MQ=255)                                                                                                                  aTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGt                               <  1:820601/149‑1 (MQ=255)                                                                                                                   tCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTc                              >  3:746006/1‑149 (MQ=255)                                                                                                                   tCGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTc                              >  3:727340/1‑149 (MQ=255)                                                                                                                    cGGTCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCa                             <  1:67712/149‑1 (MQ=255)                                                                                                                       tCACATGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATc                          >  2:677462/1‑149 (MQ=255)                                                                                                                         acaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAc                        <  3:502645/149‑1 (MQ=255)                                                                                                                         acaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAc                        <  1:549823/149‑1 (MQ=255)                                                                                                                           aTGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa                      >  1:95973/1‑149 (MQ=255)                                                                                                                           aTGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa                      >  3:367951/1‑149 (MQ=255)                                                                                                                           aTGGACACGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa                      <  2:95379/149‑1 (MQ=255)                                                                                                                                 acGGCCCTGGGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCgg                >  2:675728/1‑149 (MQ=255)                                                                                                                                  cGGCCCTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGt               >  4:387248/1‑149 (MQ=255)                                                                                                                                      ccTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtg           <  4:559317/149‑1 (MQ=255)                                                                                                                                      ccTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtg           >  2:680119/1‑149 (MQ=255)                                                                                                                                      ccTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtg           >  4:412477/1‑149 (MQ=255)                                                                                                                                       cTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgc          >  3:226528/1‑149 (MQ=255)                                                                                                                                       cTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgc          >  3:159079/1‑149 (MQ=255)                                                                                                                                       cTGAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgc          >  2:89980/1‑149 (MQ=255)                                                                                                                                         gAGTCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctg        >  2:839519/1‑149 (MQ=255)                                                                                                                                           gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCg      >  1:109699/1‑149 (MQ=255)                                                                                                                                            tCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCgg     >  1:540527/1‑149 (MQ=255)                                                                                                                                            tCATTGTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCgg     >  4:160859/1‑149 (MQ=255)                                                                                                                                               ttttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGata  >  3:197238/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GAACGACGCACTTCTGGCAGCATCTCGCGCGGAATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATA  >  NC_000913/4116561‑4116851

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