Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,874,407	T→G	D23A (GAC→GCC)	dgoK ←	2‑oxo‑3‑deoxygalactonate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,874,407	0	T	G	100.0%	238.0 / NA	63	D23A (GAC→GCC)	dgoK	2‑oxo‑3‑deoxygalactonate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (36/27);  total (36/27)

CCGACGTTGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGTCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCTGCTCTGCCGCATCGCCGTCCTGATG  >  NC_000913/3874260‑3874551                                                                                                                                                    |                                                                                                                                                 ccGACGTTGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGgcgc                                                                                                                                               >  2:750340/1‑151 (MQ=255)    aCGTTGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGccc                                                                                                                                             >  2:145222/1‑150 (MQ=255)    aCGTTGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCt                                                                                                                                            >  1:549496/1‑151 (MQ=255)     cGTTGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTg                                                                                                                                           >  2:576168/1‑151 (MQ=255)      gTTGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTg                                                                                                                                           >  1:64582/1‑150 (MQ=255)       ttGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTg                                                                                                                                           <  2:830753/149‑1 (MQ=255)        tgctgcCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGAt                                                                                                                                         <  1:737305/150‑1 (MQ=255)            gcCAACCATTCCTGCCATTCCCCCTGTCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATaaa                                                                                                                                      <  1:696513/149‑1 (MQ=255)             ccAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGcc                                                                                                                                   >  1:117195/1‑151 (MQ=255)                aCCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCa                                                                                                                                  <  2:634626/149‑1 (MQ=255)                 ccATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGgc                                                                                                                               <  1:160263/151‑1 (MQ=255)                 ccATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGgc                                                                                                                               >  2:319846/1‑151 (MQ=255)                        tGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGgcgcgc                                                                                                                           >  2:455810/1‑148 (MQ=255)                          ccATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGAtt                                                                                                                      >  1:234397/1‑151 (MQ=255)                            aTTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTgg                                                                                                                    <  1:802207/151‑1 (MQ=255)                            aTTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTgg                                                                                                                    >  2:116286/1‑151 (MQ=255)                              tACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTc                                                                                                                  >  2:591005/1‑151 (MQ=255)                               aCCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCg                                                                                                                 >  2:500959/1‑151 (MQ=255)                                ccACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGa                                                                                                                <  2:165519/151‑1 (MQ=255)                                     ggCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATcccc                                                                                                           <  1:141108/151‑1 (MQ=255)                                     ggCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATccc                                                                                                            <  2:856046/150‑1 (MQ=255)                                       cGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATcccc                                                                                                           >  2:424115/1‑149 (MQ=255)                                       cGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATcccc                                                                                                           >  2:161938/1‑149 (MQ=255)                                           ttCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCa                                                                                                      >  1:611578/1‑150 (MQ=255)                                            tctcTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAAt                                                                                                    >  2:907643/1‑151 (MQ=255)                                              tctTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTg                                                                                                  >  1:417584/1‑151 (MQ=255)                                               ctTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAAtt                                                                                                   <  2:72024/149‑1 (MQ=255)                                                 tCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTg                                                                                                  <  2:356957/148‑1 (MQ=255)                                                  cACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGAt                                                                                              <  1:19282/151‑1 (MQ=255)                                                     gCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGt                                                                                            >  1:555774/1‑150 (MQ=255)                                                     gCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTa                                                                                           <  1:175772/151‑1 (MQ=255)                                                     gCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTa                                                                                           <  2:171050/151‑1 (MQ=255)                                                        aGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCg                                                                                        >  2:816523/1‑151 (MQ=255)                                                        aGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGAAGCg                                                                                        <  1:826899/151‑1 (MQ=255)                                                         gtcggtcgTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCg                                                                                        >  1:67422/1‑150 (MQ=255)                                                            ggtcgTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCt                                                                                    >  2:258002/1‑151 (MQ=255)                                                                   aCTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCAtgt                                                                             >  2:86820/1‑151 (MQ=255)                                                                      tCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCAtgtg                                                                            <  1:406643/149‑1 (MQ=255)                                                                        tGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGAttt                                                                        >  2:864473/1‑151 (MQ=255)                                                                        tGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGAttt                                                                        >  2:652764/1‑151 (MQ=255)                                                                         gCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTc                                                                       <  1:4852/151‑1 (MQ=255)                                                                            aaCACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCtt                                                                    >  1:491660/1‑151 (MQ=255)                                                                            aaCACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCtt                                                                    >  2:687013/1‑151 (MQ=255)                                                                            aaCACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCtt                                                                    >  1:505748/1‑151 (MQ=255)                                                                                  gCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCt                                                              >  1:902319/1‑151 (MQ=255)                                                                                     gCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCg                                                           <  2:487068/151‑1 (MQ=255)                                                                                       cGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCgtg                                                         >  1:491805/1‑151 (MQ=255)                                                                                         gagaTTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTc                                                       >  1:843450/1‑151 (MQ=255)                                                                                           gaTTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCg                                                      >  2:632406/1‑150 (MQ=255)                                                                                            aTTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCg                                                      >  1:66173/1‑149 (MQ=255)                                                                                               ttCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTg                                                 >  1:785624/1‑151 (MQ=255)                                                                                                tCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTg                                                 <  1:108240/150‑1 (MQ=255)                                                                                                 ccGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGc                                               <  2:552215/151‑1 (MQ=255)                                                                                                    ttCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGAt                                            >  1:384258/1‑151 (MQ=255)                                                                                                      cAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCa                                          <  1:369940/151‑1 (MQ=255)                                                                                                        ggCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATg                                        <  1:712448/151‑1 (MQ=255)                                                                                                          cgcgTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGt                                      <  1:753403/151‑1 (MQ=255)                                                                                                                aCGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAgcgc                                <  2:368769/151‑1 (MQ=255)                                                                                                                     tGCTTCTGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCt                           >  2:672880/1‑151 (MQ=255)                                                                                                                           tGATTGCCTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCTGCTCTg                     >  1:768770/1‑151 (MQ=255)                                                                                                                                 ccTGCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCTGCTCTGCCGCAt               >  2:629779/1‑151 (MQ=255)                                                                                                                                    gCTCTCCAGGCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCTGCTCTGCCGCATc              <  2:263874/149‑1 (MQ=255)                                                                                                                                             gCAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCTGCTCTGCCGCATCGCCGTCCTGAt   <  2:475822/151‑1 (MQ=255)                                                                                                                                              cAGTGGGCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCTGCTCTGCCGCATCGCCGTCCTGATg  <  2:582900/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CCGACGTTGCTGCCAACCATTCCTGCCATTACCACTGGCGTTTTCTCTTCACGCCAGTCGGTCGTGACTTCTGCTAACACCGCAGCCGGAGATTTTCCGTTCAGGCGCGTGACGCCTGCTTCTGATTGCCTGCTCTCCAGGCAGTGGTCGCCCTGATAAAGCCAGGCGCGCAGATTGGTCGATCCCCAGTCAATTGCGATGTAGCGAGCTGTCATGTGATTTCCTTTAACCTTCGTGTCGAGCTGGCGATCATGGTAAGCGCCGCCTGCTCTGCCGCATCGCCGTCCTGATG  >  NC_000913/3874260‑3874551

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