Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,943,712	T→C	W126R (TGG→CGG)	csdA →	cysteine sulfinate desulfinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,943,712	0	T	C	100.0%	259.7 / NA	69	W126R (TGG→CGG)	csdA	cysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/37);  total (32/37)

TGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCG  >  NC_000913/2943576‑2943846                                                                                                                                         |                                                                                                                                       tgatAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGc                                                                                                                                     <  2:55743/140‑1 (MQ=255)  gatAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCt                                                                                                                                    >  1:224762/1‑140 (MQ=255)          tATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCg                                                                                                                            <  2:728921/140‑1 (MQ=255)          tATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCg                                                                                                                            <  2:479245/140‑1 (MQ=255)          tATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCg                                                                                                                            >  1:916957/1‑140 (MQ=255)          tATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCg                                                                                                                            >  1:301538/1‑140 (MQ=255)            tCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGcc                                                                                                                          <  1:913341/140‑1 (MQ=255)            tCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGcc                                                                                                                          >  2:913341/1‑140 (MQ=255)             cGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGccc                                                                                                                         >  2:961911/1‑140 (MQ=255)                 tGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCcaaca                                                                                                                     >  2:430866/1‑140 (MQ=255)                   gACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACaaa                                                                                                                   >  2:1041140/1‑140 (MQ=255)                          ggCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGc                                                                                                            <  1:1041140/140‑1 (MQ=255)                            caccacTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCa                                                                                                          <  1:741321/140‑1 (MQ=255)                            caccacTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCa                                                                                                          <  2:601789/140‑1 (MQ=255)                            caccacTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCa                                                                                                          <  1:961911/140‑1 (MQ=255)                            caccacTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCa                                                                                                          <  1:688376/140‑1 (MQ=255)                                acTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAgt                                                                                                      <  2:812130/140‑1 (MQ=255)                                     aTCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGa                                                                                                 <  2:301538/140‑1 (MQ=255)                                         aTCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAAtt                                                                                             >  2:302097/1‑140 (MQ=255)                                            aaCATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGcc                                                                                          <  2:861107/140‑1 (MQ=255)                                            aaCATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGcc                                                                                          <  1:395703/140‑1 (MQ=255)                                            aaCATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGcc                                                                                          >  1:879887/1‑140 (MQ=255)                                               aTGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCt                                                                                       <  1:77442/140‑1 (MQ=255)                                                   tggCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAAt                                                                                   >  1:554651/1‑140 (MQ=255)                                                      cacaATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCATGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATgcg                                                                                >  1:25232/1‑140 (MQ=255)                                                      cacaATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATgcg                                                                                >  2:742986/1‑140 (MQ=255)                                                        caATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCa                                                                              >  2:621540/1‑140 (MQ=255)                                                         aaTGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAg                                                                             <  2:266729/140‑1 (MQ=255)                                                           tGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCg                                                                           >  1:1084281/1‑140 (MQ=255)                                                            gCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGa                                                                          <  2:642841/140‑1 (MQ=255)                                                            gCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGa                                                                          <  2:101471/140‑1 (MQ=255)                                                                tGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGc                                                                      >  1:1126495/1‑140 (MQ=255)                                                                 gcgcGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGcc                                                                     <  1:311005/140‑1 (MQ=255)                                                                   gcgtccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCa                                                                              <  1:937375/129‑1 (MQ=255)                                                                   gcgtccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCa                                                                              >  2:937375/1‑129 (MQ=255)                                                                   gcgtccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCgg                                                                   <  1:555976/140‑1 (MQ=255)                                                                   gcgtccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCgg                                                                   <  1:430866/140‑1 (MQ=255)                                                                     gtccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGAt                                                                 >  2:952560/1‑140 (MQ=255)                                                                      tccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATg                                                                >  2:989356/1‑140 (MQ=255)                                                                       ccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGt                                                               >  2:832275/1‑140 (MQ=255)                                                                       ccgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGt                                                               >  1:41570/1‑140 (MQ=255)                                                                        cgcgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTc                                                              >  2:46389/1‑140 (MQ=255)                                                                          cgtcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGa                                                            >  1:572621/1‑140 (MQ=255)                                                                            tcTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCTTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGAtt                                                          >  2:307146/1‑140 (MQ=255)                                                                             cTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGAttt                                                         <  2:1126495/140‑1 (MQ=255)                                                                             cTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGAttt                                                         <  1:989356/140‑1 (MQ=255)                                                                                  aCCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCATGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGc                                                    <  2:25232/140‑1 (MQ=255)                                                                                  aCCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGc                                                    <  2:879887/140‑1 (MQ=255)                                                                                    cGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCa                                                  <  2:408492/140‑1 (MQ=255)                                                                                    cGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCa                                                  <  2:916957/140‑1 (MQ=255)                                                                                    cGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCa                                                  <  2:871130/140‑1 (MQ=255)                                                                                    cGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCa                                                  <  1:656310/140‑1 (MQ=255)                                                                                     gggCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATttgt                                                       >  1:845376/1‑134 (MQ=255)                                                                                     gggCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATttgt                                                       <  2:845376/134‑1 (MQ=255)                                                                                      ggCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGa                                                <  2:554651/140‑1 (MQ=255)                                                                                      ggCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGa                                                <  1:302097/140‑1 (MQ=255)                                                                                               attattGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTAc                                       >  2:326940/1‑140 (MQ=255)                                                                                                        aGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAg                              >  2:187129/1‑140 (MQ=255)                                                                                                        aGCGTGGCAGAACACCACGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAAGGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAg                              >  2:990287/1‑140 (MQ=255)                                                                                                                    caccacGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGc                  <  1:326940/140‑1 (MQ=255)                                                                                                                       cacGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGtt               >  2:607387/1‑140 (MQ=255)                                                                                                                       cacGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGtt               <  1:1122512/140‑1 (MQ=255)                                                                                                                         cGCCAACCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTgg             >  1:281049/1‑140 (MQ=255)                                                                                                                             aaCCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCa         <  1:952560/140‑1 (MQ=255)                                                                                                                             aaCCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCa         <  2:572621/140‑1 (MQ=255)                                                                                                                             aaCCTCGTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCa         <  1:187129/140‑1 (MQ=255)                                                                                                                                   gTCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTc   >  2:857168/1‑140 (MQ=255)                                                                                                                                    tCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCg  >  1:722085/1‑140 (MQ=255)                                                                                                                                    tCCCCCGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCg  >  1:808327/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       TGATAAAACTATCGTCTGGACGCGCGGCACCACTGAATCCATCAACATGGTGGCACAATGCTATGCGCGTCCGCGTCTGCAACCGGGCGATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCG  >  NC_000913/2943576‑2943846

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