Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,552,845	G→A	10.5%	intergenic (+25/‑56)	uxuA → / → uxuB	mannonate hydrolase/D‑mannonate oxidoreductase, NAD‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,552,845	0	G	A	10.5%	225.6 / 14.6	76	intergenic (+25/‑56)	uxuA/uxuB	mannonate hydrolase/D‑mannonate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base A (5/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

TCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCATCTCGGTTGCGGG  >  NC_000913/4552697‑4552993                                                                                                                                                      |                                                                                                                                                     tcaGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCg                                                                                                                                                      >  4:206904/1‑149 (MQ=255) tcaGATGCTGGACGACCTGAAGAAGAAAACAAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCg                                                                                                                                                      >  1:377320/1‑149 (MQ=255)  caGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGt                                                                                                                                                     >  1:408329/1‑149 (MQ=255)  caGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGt                                                                                                                                                     <  3:61487/149‑1 (MQ=255)       gCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCgg                                                                                                                                                >  2:477014/1‑149 (MQ=255)       gCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCACGGC‑GCGCAGCCATGCCgg                                                                                                                                                <  1:81046/149‑1 (MQ=255)          ggACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTc                                                                                                                                             >  4:448697/1‑149 (MQ=255)          ggACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTc                                                                                                                                             <  3:469461/149‑1 (MQ=255)              gacCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCcttc                                                                                                                                         <  1:178150/149‑1 (MQ=255)               acCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCcttct                                                                                                                                        <  4:348739/149‑1 (MQ=255)                 cTGAAGAAGAAAACCAACCCAGGTTACTCCGCGATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGACTGCGC‑GCCGTGCCGGTTCCttcttc                                                                                                                                      >  2:91454/1‑149 (MQ=255)                   gaagaagaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGTGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCt                                                                                                                                    >  4:310059/1‑149 (MQ=255)                   gaagaagaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGc                                                                                                                                                   <  2:334145/134‑1 (MQ=255)                   gaagaagaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGc                                                                                                                                                   >  1:334145/1‑134 (MQ=255)                    aagaagaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCtt                                                                                                                                   >  2:188282/1‑149 (MQ=255)                        agaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCg                                                                                                                               >  3:481653/1‑149 (MQ=255)                         gaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGt                                                                                                                              >  4:194426/1‑149 (MQ=255)                         gaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCACGGC‑GCGCAGCCATGCCGGtggat                                                                                                                                           <  3:261282/136‑5 (MQ=255)                         gaaAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCACGGC‑GCGCAGCCATGCCGGtggat                                                                                                                                           >  4:261282/1‑132 (MQ=255)                               caaccCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCActct                                                                                                                        >  3:452070/1‑149 (MQ=255)                                aaccCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCtt                                                                                                                       <  3:111497/149‑1 (MQ=255)                                    cAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGaa                                                                                                                   >  1:85015/1‑149 (MQ=255)                                       gTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGAc                                                                                                                <  3:299733/149‑1 (MQ=255)                                       gTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGAc                                                                                                                <  3:206904/149‑1 (MQ=255)                                            tCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGAtt                                                                                                           >  2:99990/1‑149 (MQ=255)                                            tCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGAtt                                                                                                           >  4:423117/1‑149 (MQ=255)                                             ccGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTTGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTc                                                                                                          >  4:99045/1‑149 (MQ=255)                                              cGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGAc                                                                                                                <  2:282271/142‑1 (MQ=255)                                              cGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGAc                                                                                                                >  1:282271/1‑142 (MQ=255)                                              cGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCt                                                                                                         >  3:460603/1‑149 (MQ=255)                                               gCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGTGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTg                                                                                                        <  3:310059/149‑1 (MQ=255)                                               gCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTg                                                                                                        <  1:188282/149‑1 (MQ=255)                                                 aaTTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGa                                                                                                      <  2:408329/149‑1 (MQ=255)                                                 aaTTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGa                                                                                                      <  3:448697/149‑1 (MQ=255)                                                      gtcgtcTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTa                                                                                                 >  4:562423/1‑149 (MQ=255)                                                       tcgtcTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTAc                                                                                                >  4:472557/1‑149 (MQ=255)                                                            tGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGa                                                                                           <  4:452070/149‑1 (MQ=255)                                                             gAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGac                                                                                          >  2:161244/1‑149 (MQ=255)                                                             gAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGac                                                                                          <  4:481653/149‑1 (MQ=255)                                                              aaaGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGact                                                                                         >  1:169108/1‑149 (MQ=255)                                                               aaGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGacta                                                                                        <  1:99990/149‑1 (MQ=255)                                                                aGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGactac                                                                                       <  3:562423/149‑1 (MQ=255)                                                                aGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGactac                                                                                       >  3:207722/1‑149 (MQ=255)                                                                      ggCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTAttgt                                                                                 >  1:437041/1‑149 (MQ=255)                                                                         cGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGa                                                                              <  1:477014/149‑1 (MQ=255)                                                                         cGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGa                                                                              <  2:85015/149‑1 (MQ=255)                                                                          gAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGAc                                                                             <  1:611352/149‑1 (MQ=255)                                                                            aGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAg                                                                           >  1:170957/1‑149 (MQ=255)                                                                             gTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGc                                                                          >  1:86227/1‑149 (MQ=255)                                                                              ttCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGa                                                                                           <  3:441974/131‑1 (MQ=255)                                                                              ttCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGa                                                                                           >  4:441974/1‑131 (MQ=255)                                                                                     tgtCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGcc                                                                  <  3:335509/149‑1 (MQ=255)                                                                                           aCTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATGTCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATACTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTgc                                                            <  3:423117/149‑1 (MQ=255)                                                                                                 gATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTgcccgccc                                                      <  3:194426/149‑1 (MQ=255)                                                                                                  ctCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATAGTTGACAGCAATCTGCCGGTTGcccgcccg                                                     <  2:377320/148‑1 (MQ=255)                                                                                                       gcgcgcTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTg                                                                                                        <  1:235088/93‑1 (MQ=255)                                                                                                       gcgcgcTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTg                                                                                                        >  2:235088/1‑93 (MQ=255)                                                                                                             tttctttAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCa                                          >  4:586208/1‑149 (MQ=255)                                                                                                             tttctttAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCa                                          >  4:269229/1‑149 (MQ=255)                                                                                                                              atCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCAc                         <  3:269229/149‑1 (MQ=255)                                                                                                                              atCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCAc                         >  2:429399/1‑149 (MQ=255)                                                                                                                              atCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCAc                         >  1:476365/1‑149 (MQ=255)                                                                                                                               tCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACg                        >  4:295157/1‑149 (MQ=255)                                                                                                                                gaaCCGGCACGGC‑GCGCAGCCATGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGc                       >  3:253361/3‑149 (MQ=255)                                                                                                                                  aCCGGCACGGC‑GCGCAGCCATGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCAt                     >  1:358639/1‑149 (MQ=255)                                                                                                                                   ccGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCAtt                    <  1:292199/149‑1 (MQ=255)                                                                                                                                   ccGGCACGGC‑GCGCAGCCATGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGa                                                                                                      <  1:136503/67‑1 (MQ=25)                                                                                                                                   ccGGCACGGC‑GCGCAGCCATGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGa                                                                                                      >  2:136503/1‑67 (MQ=25)                                                                                                                                   ccGGCACGGC‑GCGCAGCCATGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCAtt                    >  2:13754/1‑149 (MQ=255)                                                                                                                                      gCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATtgtg                 >  2:161349/1‑149 (MQ=255)                                                                                                                                        aTGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCa               >  2:127173/1‑149 (MQ=255)                                                                                                                                             tGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGa                              >  2:74332/1‑129 (MQ=255)                                                                                                                                             tGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGa                              <  1:74332/129‑1 (MQ=255)                                                                                                                                             tGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCATCTCg          >  2:517510/1‑149 (MQ=255)                                                                                                                                               cgc‑gcCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTAtt                                                                                   <  3:83306/74‑1 (MQ=255)                                                                                                                                               cgc‑gcCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTAtt                                                                                   >  4:83306/1‑74 (MQ=255)                                                                                                                                                     cGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCATCTCGGTTGCgg   <  2:437041/149‑1 (MQ=255)                                                                                                                                                     cGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCATCTCGGTTGCgg   <  4:207722/149‑1 (MQ=255)                                                                                                                                                      gTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCATCTCGGTTGCggg  <  1:571799/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TCAGATGCTGGACGACCTGAAGAAGAAAACCAACCCAGGTTACTCCGCAATTGGTCGTCTGAAAGGCCTGGCCGAAGTTCGCGGTGTCGAACTGGCGATCCAGCGCGCTTTCTTTAGCCGTTAATATCCACCGGCATGGCTGCGC‑GCCGTGCCGGTTCCTTCTTCCTTGCCGTCACTCTTTGAAGACGGATTCTGGAGTTTACGATGACTACTATTGTTGACAGCAATCTGCCGGTTGCCCGCCCGTCATGGGATCATTCTCGTCTGGAATCACGCATTGTGCATCTCGGTTGCGGG  >  NC_000913/4552697‑4552993

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