Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,109	T→G	Q38P (CAA→CCA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,109	0	T	G	98.6%	229.4 / ‑6.1	70	Q38P (CAA→CCA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (37/32);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAG  >  NC_000913/4116963‑4117249                                                                                                                                                   |                                                                                                                                             tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGc                                                                                                                                            <  1:897882/149‑1 (MQ=255)   gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTc                                                                                                                                          >  2:724409/1‑149 (MQ=255)        aCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAAtt                                                                                                                                     >  1:499229/1‑149 (MQ=255)        aCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAAtt                                                                                                                                     <  2:601354/149‑1 (MQ=255)        aCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAAtt                                                                                                                                     <  1:1191345/149‑1 (MQ=255)         cGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTc                                                                                                                                    >  1:341078/1‑149 (MQ=255)            tGGTTCGTAATACCGATAGCTGCAATTTGGTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTcgcg                                                                                                                                 <  2:804306/149‑1 (MQ=255)            tGGTGCGTAATAGCGATAGCTGCAATTTGATCGGAATTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGTGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTcgcg                                                                                                                                 <  1:771633/149‑1 (MQ=255)             ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTCGATGGGCTCAAATTCgcgc                                                                                                                                >  2:781783/1‑149 (MQ=255)                tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGAGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGc                                                                                                                             >  1:141690/1‑149 (MQ=255)                tCGTAATACCGATAGCCGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGc                                                                                                                             <  2:230088/149‑1 (MQ=255)                   tcaTACCGAGCGCGTCAATTTGATCGGAACTGATATCGGCTTGCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACGTGTACTTGGTTAGATTGGCTCAAATTCGCGCTGCGac                                                                                                                          <  2:1306447/147‑1 (MQ=255)                   tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGac                                                                                                                          >  2:499704/1‑149 (MQ=255)                    aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGaca                                                                                                                         <  1:1028491/149‑1 (MQ=255)                          gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTaa                                                                                                                   <  1:1028802/149‑1 (MQ=255)                             aGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGa                                                                                                                >  1:28228/1‑149 (MQ=255)                              gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGat                                                                                                               <  2:499229/149‑1 (MQ=255)                                     tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGCCCCGCTAATGATATTGGCa                                                                                                        >  1:963728/1‑149 (MQ=255)                                     tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCa                                                                                                        >  2:165063/1‑149 (MQ=255)                                         aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGGAGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGATCAAATCCGCGCTGCGACACGCTACTGATATTGGCCTCGt                                                                                                    >  1:538265/1‑149 (MQ=255)                                          tCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGGCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTAGCATCGTg                                                                                                   >  2:1284267/1‑149 (MQ=255)                                             gAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGAGTGGAGCTTTGGGTGGCCCGGATTTCCATTGGGTCTTCTTCTACCCAACCTGGTTTTGGGTAGGTTGGCTCAAATTCGGGTTGCGACACGCTAATGACATTGGCATCGTGATc                                                                                                >  2:819242/1‑149 (MQ=255)                                              aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGCTATACATTGGGTCGTGTTCTACACAACCTGGGATTGGGTAGTTTTGCTTAAATTAGCGATGCGACATGCTTATGATATTGGCATCGTGATcc                                                                                               >  1:854297/1‑149 (MQ=255)                                                 tGATATCGCCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                                                                            <  2:969463/149‑1 (MQ=255)                                                  gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTa                                                                                           <  1:901933/149‑1 (MQ=255)                                                  gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTa                                                                                           <  2:702668/149‑1 (MQ=255)                                                     atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAGCCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGa                                                                                        <  2:583868/149‑1 (MQ=255)                                                     atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGa                                                                                        <  2:341078/149‑1 (MQ=255)                                                        ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcg                                                                                     >  2:367852/1‑149 (MQ=255)                                                                ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGCCCGCGCGGGAg                                                                             >  1:614610/1‑149 (MQ=255)                                                                ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                             <  2:875491/149‑1 (MQ=255)                                                                 cAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTGCCAGGGGGTCGTGTTCTACATAAGAGGGATTTGGGTAGTTTGGCTCAAATTCGCGCTGCGACACGCTAATGATCTTGGCATCGTGATCCATTACGCCCGCGCGGGAGc                                                                            >  1:383182/1‑149 (MQ=255)                                                                  aGCACTTCTACCAGTATGAAGCTTTGGGTGGCCCAGAAATATATGGGGTAGTGTTCTACCCAACCTGGTTTTGGGGAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGCGATCCATTACGACCGCGCGGGAGCt                                                                           >  1:999465/1‑149 (MQ=255)                                                                  aGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCt                                                                           >  1:508884/1‑149 (MQ=255)                                                                   gCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCCGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtg                                                                          <  1:392223/149‑1 (MQ=255)                                                                   gCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACGCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGCTATTGGCATCGGGATCCATTACGACCGCGCGGGAGCtg                                                                          <  2:647519/149‑1 (MQ=255)                                                                      ctgcgaCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCACTGGGTCGTGTTCTCCCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                                       <  2:697919/144‑1 (MQ=255)                                                                        tctcCCAGCGTGGAGCTTTGGGTGGCCCAGTTTTCCATAGGGACGTGTTCTTCCCAACCTGGTTTTGGGCAGATTGGCTCAAATTCGCGCTGCGGCACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTggtggt                                                                     <  2:182692/149‑1 (MQ=255)                                                                           aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGcc                                                                  <  1:724409/149‑1 (MQ=255)                                                                           aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCCGCGACACGCTAATGATATTGGCACCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGcc                                                                  <  2:854297/149‑1 (MQ=255)                                                                             cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATCCGCGCTGCGACACGCTAATGATATTGGGATCGTGATCCATTACGACCGCGCGGGAGCTGGTAGCGCCCt                                                                >  2:453289/1‑149 (MQ=255)                                                                              agcgtggagcgtggGGTGGCCCAAATTTCCATTCGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCTATTACGACCGCGCGGGAGCTGGTGGTGCCCTg                                                               <  1:1321643/149‑1 (MQ=255)                                                                                 gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTc                                                            <  1:367852/149‑1 (MQ=255)                                                                                   ggAGCTTTAGGTGGCCCAGATTTCCATTGGGTCATGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCCCGCGGGGGCTGGTGGTGCCCTGGTCgg                                                          >  1:992342/1‑148 (MQ=255)                                                                                     ggCTTTGGGTGGCCCAGATTTCCATTGGGGCGTGTTCTACCCACCCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc                                                        <  2:299661/148‑1 (MQ=255)                                                                                     aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGGGCTGGTGGTGCCCTGGTCGGgc                                                        >  1:1204563/1‑149 (MQ=255)                                                                                     aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGGGTTCTACCAAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc                                                        <  2:237248/149‑1 (MQ=255)                                                                                      gCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGGGCTGGGGGTGCCCTGGTCGAgcg                                                       >  1:1021949/1‑149 (MQ=255)                                                                                       cTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACCCGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGACCTGGTCGAgcgc                                                      >  2:326853/1‑149 (MQ=255)                                                                                        tttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGAGCCATGACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCa                                                     >  1:712410/1‑149 (MQ=255)                                                                                               ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGAAAATTGGCTCAAATTCGCGCTGCGACCCGCTAATGATCTTGGCATCGTTATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatat                                              <  1:37941/149‑1 (MQ=255)                                                                                                  ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtt                                           >  1:1182416/1‑149 (MQ=255)                                                                                                   cAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTCTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTCATGATATTGGCATCGTGCTCCATTACGACCGCGCGGGAGCTGGTGGTGCCCAGGTCGAGCGCAACGACATAttt                                          >  2:275615/1‑149 (MQ=255)                                                                                                    aGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttt                                         <  1:499704/149‑1 (MQ=255)                                                                                                       tttCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGCCCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAttttttt                                      >  2:600798/1‑149 (MQ=255)                                                                                                          ccATTGTGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAACTTCGAACCACTACACGCTAATGATATTGGCAGCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                   <  1:651720/149‑1 (MQ=255)                                                                                                          ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                   >  2:1012993/1‑149 (MQ=255)                                                                                                            gttGGGTCGTGTTCTACCCCACCTGGTTGTGGGTAGATTGGCTCAACTGCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGGGCTGGTGGTGCCCTTGTCGAGCGCAACGATATATTTTTTTTCAGt                                 <  2:519094/148‑1 (MQ=255)                                                                                                            aTTGGGTCGTGTTCTACCCAACCTGGTTTTGGGGCGATTGGCTCAGATTCGCGCTGCGACACGCTACTGATATTGGCATCGTGATCCCTTACGACCGCGCGGGCGCTGGTGGTGCCCCGGTCGAGCGCAAGGGTATAttttttttagtt                                 >  2:262077/1‑145 (MQ=255)                                                                                                                     tgtTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                        >  1:1070648/1‑149 (MQ=255)                                                                                                                       gtCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAtt                      <  1:238484/148‑1 (MQ=255)                                                                                                                            cccAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTccc                 <  2:614610/149‑1 (MQ=255)                                                                                                                             ccAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCg                >  1:971475/1‑149 (MQ=255)                                                                                                                               aaCCTGGTTTTGGGTAGATTGGCTCACATTCACGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCATTCATGTTTAATTGTCCCGTa              >  2:967898/1‑149 (MQ=255)                                                                                                                                  cTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCCTGTTTAATTGGCCCgacctc           >  2:811199/1‑144 (MQ=255)                                                                                                                                   tGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCa          >  2:805542/1‑149 (MQ=255)                                                                                                                                   tGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGGCATGTTTAATTGTCCCGTAGTcc          >  1:1124015/1‑148 (MQ=255)                                                                                                                                     gTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCata        >  2:1152025/1‑149 (MQ=255)                                                                                                                                          gggTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGAAGTCATATTACa   >  2:186414/1‑149 (MQ=255)                                                                                                                                           ggTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAg  >  2:1281214/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                             TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAG  >  NC_000913/4116963‑4117249

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