Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,202	A→G	D132D (GAT→GAC)	yjcS ←	alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,202	0	A	G	100.0%	256.8 / NA	74	D132D (GAT→GAC)	yjcS	alkyl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/39);  total (35/39)

TAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCG  >  NC_000913/4306065‑4306338                                                                                                                                          |                                                                                                                                         tAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTc                                                                                                                                        >  1:798851/1‑140 (MQ=255)     tAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTgg                                                                                                                                    <  1:727188/140‑1 (MQ=255)     tAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTgg                                                                                                                                    <  2:387171/140‑1 (MQ=255)       aCGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGcc                                                                                                                                  >  1:651493/1‑140 (MQ=255)        cGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCcg                                                                                                                                 >  1:1010828/1‑140 (MQ=255)        cGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCcg                                                                                                                                 <  2:387334/140‑1 (MQ=255)         ggCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCcgc                                                                                                                                <  1:306486/140‑1 (MQ=255)         ggCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGCTAGAGATGTCCTGGCcgc                                                                                                                                >  2:662318/1‑140 (MQ=255)         ggCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAAAGTAATGTTAGAGATGTCCTGGCcgc                                                                                                                                >  2:1324170/1‑140 (MQ=255)          gCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCcgcg                                                                                                                               <  1:653988/140‑1 (MQ=255)           caacaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCgcgc                                                                                                                              <  2:436649/140‑1 (MQ=255)           caacaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCgcgc                                                                                                                              <  2:210011/140‑1 (MQ=255)           caacaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCgcgc                                                                                                                              >  1:859668/1‑140 (MQ=255)             acaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTTCTGGCCGCGCAc                                                                                                                            <  2:1076020/140‑1 (MQ=255)             acaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCcgc                                                                                                                                <  2:706828/136‑1 (MQ=255)             acaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCcgc                                                                                                                                >  1:706828/1‑136 (MQ=255)              caaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCAcc                                                                                                                           >  1:1070416/1‑140 (MQ=255)                   ggTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGGGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATa                                                                                                                      >  1:285439/1‑140 (MQ=255)                   ggTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAGTGTTAGAGATGTCCTGGCCGCGCACCTGATa                                                                                                                      >  1:307095/1‑140 (MQ=255)                   ggTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATa                                                                                                                      <  2:1070416/140‑1 (MQ=255)                         tGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACAtttt                                                                                                                <  2:905251/140‑1 (MQ=255)                              aCGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCgg                                                                                                           <  2:651493/140‑1 (MQ=255)                              aCGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCgg                                                                                                           <  2:798851/140‑1 (MQ=255)                              aCGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCgg                                                                                                           <  1:628751/140‑1 (MQ=255)                              aCGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCgg                                                                                                           <  1:927360/140‑1 (MQ=255)                                   gCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACt                                                                                                      <  1:1040957/140‑1 (MQ=255)                                   gCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACt                                                                                                      <  2:859668/140‑1 (MQ=255)                                       gAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTCCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGa                                                                                                  >  1:486502/1‑140 (MQ=255)                                        aaGTGAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCATCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGaa                                                                                                 <  2:1408631/140‑1 (MQ=255)                                        aaGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGaa                                                                                                 <  2:310636/140‑1 (MQ=255)                                        aaGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGaa                                                                                                 <  2:256495/140‑1 (MQ=255)                                        aaGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGaa                                                                                                 <  1:867587/140‑1 (MQ=255)                                        aaGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGaa                                                                                                 <  1:835344/140‑1 (MQ=255)                                        aaGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGaa                                                                                                 <  1:380398/140‑1 (MQ=255)                                        aaGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAAAGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGaa                                                                                                 <  1:129272/140‑1 (MQ=255)                                            aaaGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACt                                                                                                      <  2:365536/131‑1 (MQ=255)                                            aaaGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACt                                                                                                      >  1:365536/1‑131 (MQ=255)                                            aaaGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAgg                                                                                             >  2:1250899/1‑140 (MQ=255)                                              aGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGcc                                                                                           >  1:1112208/1‑140 (MQ=255)                                                 tCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGa                                                                                        >  1:560090/1‑140 (MQ=255)                                                   aaGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGaaa                                                                                      >  1:801016/1‑140 (MQ=255)                                                         ggCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGt                                                                                >  1:459705/1‑140 (MQ=255)                                                          gCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGCTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGtt                                                                               <  1:662318/140‑1 (MQ=255)                                                          gCTTTTGCGGCAGGCCGCGTCACCAGCGGGGCGATAACAATAATTCCTTTCTCGCCCTCAAAGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGGACAGGCCAGAAATACCGtt                                                                               >  1:1127198/1‑140 (MQ=255)                                                               tGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGAttt                                                                          >  2:1229686/1‑140 (MQ=255)                                                                     aGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACt                                                                    <  2:1083235/140‑1 (MQ=255)                                                                     aGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACt                                                                    >  2:680880/1‑140 (MQ=255)                                                                         ggcgTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACg                                                                >  2:1258901/1‑140 (MQ=255)                                                                           cgTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGcc                                                              <  2:1010828/140‑1 (MQ=255)                                                                            gTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCca                                                             <  2:1127198/140‑1 (MQ=255)                                                                            gTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCca                                                             <  2:801016/140‑1 (MQ=255)                                                                                ccAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAgg                                                         <  1:1258901/140‑1 (MQ=255)                                                                                  aGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCt                                                       <  2:459705/140‑1 (MQ=255)                                                                                     gggTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTggg                                                    >  2:880975/1‑140 (MQ=255)                                                                                       gTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATAGATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGt                                                  >  1:1012235/1‑140 (MQ=255)                                                                                       gTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGt                                                  >  1:185293/1‑140 (MQ=255)                                                                                       gTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGt                                                  >  2:1096855/1‑140 (MQ=255)                                                                                           aTAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTAc                                              >  2:675725/1‑140 (MQ=255)                                                                                                 ataatGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTc                                        >  2:2944/1‑140 (MQ=255)                                                                                                  taatGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTcc                                       >  1:945418/1‑140 (MQ=255)                                                                                                     tGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGgc                                    >  2:447785/1‑140 (MQ=255)                                                                                                     tGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGgc                                    >  1:18043/1‑140 (MQ=255)                                                                                                            ctcGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGc                               >  2:993459/1‑138 (MQ=255)                                                                                                            ctcGCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGc                               <  1:993459/138‑1 (MQ=255)                                                                                                               gCCCTCAACGAACGTAGTGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCAt                          <  2:307095/140‑1 (MQ=255)                                                                                                               gCCCTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCAt                          <  2:486502/140‑1 (MQ=255)                                                                                                                 ccTCAACGAACGTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTa                        >  2:880902/1‑140 (MQ=255)                                                                                                                     aacgaacgTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAatat                    <  2:285439/140‑1 (MQ=255)                                                                                                                       cgaacgTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCa                  <  2:18043/140‑1 (MQ=255)                                                                                                                        gaacgTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCaa                 <  1:1006212/140‑1 (MQ=255)                                                                                                                           cgTAATGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAAtt              >  2:616398/1‑140 (MQ=255)                                                                                                                               aTGTTAGAGATGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTa          >  1:804230/1‑140 (MQ=255)                                                                                                                                     gagaTGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTc    >  1:165750/1‑140 (MQ=255)                                                                                                                                       gaTGTCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCcg  <  1:880902/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         TAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCG  >  NC_000913/4306065‑4306338

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