Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,891,254	C→T	S215S (TCG→TCA)	cysJ ←	sulfite reductase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,891,254	0	C	T	100.0%	216.1 / NA	62	S215S (TCG→TCA)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/33);  total (29/33)

TTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTG  >  NC_000913/2891118‑2891389                                                                                                                                         |                                                                                                                                        tttCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGaa                                                                                                                                      >  2:510916/1‑140 (MQ=255)          ggCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCcgcg                                                                                                                             <  1:1181141/140‑1 (MQ=255)                 gTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGgcg                                                                                                                      <  2:777835/140‑1 (MQ=255)                   ttttttCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCa                                                                                                                    <  1:726596/139‑1 (MQ=255)                   cTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCa                                                                                                                    <  1:1134647/140‑1 (MQ=255)                    tttttCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCAc                                                                                                                   <  1:510916/140‑1 (MQ=255)                           gagTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCa                                                                                                                    >  2:554760/1‑132 (MQ=255)                           gagTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCa                                                                                                                    <  1:554760/132‑1 (MQ=255)                           gagTTACGCCCGGTAATTTTCTGGTTAACAGACAGGCTAGCCACCATCGTCGCGTCTTTGCTGTTCGGGTTTGTGTGGATGTCTTTTACCTGGCCAGTAGGGACGGTTTGTGAAGGGGCCGGTGCAGGCGCACTGGAttt                                                                                                            >  2:726596/1‑140 (MQ=255)                            agTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCa                                                                                                                    <  2:1211484/131‑1 (MQ=255)                            agTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCa                                                                                                                    >  1:1211484/1‑131 (MQ=255)                                  gCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCGGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCa                                                                                                     <  1:206513/140‑1 (MQ=255)                                  gCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCa                                                                                                     <  1:936329/140‑1 (MQ=255)                                      ggTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCaa                                                                                                 >  2:1182052/1‑140 (MQ=255)                                                 ggTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGgcg                                                                                      <  1:791783/140‑1 (MQ=255)                                                     aaCAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGcc                                                                                  >  2:705795/1‑140 (MQ=255)                                                      aCAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGccc                                                                                 >  1:536020/1‑139 (MQ=255)                                                       cAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCAc                                                                                <  2:41221/140‑1 (MQ=255)                                                            agGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCggg                                                                                        >  1:514549/1‑127 (MQ=255)                                                            agGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCggg                                                                                        <  2:514549/127‑1 (MQ=255)                                                            agGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCt                                                                           >  1:306493/1‑140 (MQ=255)                                                               cTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGgc                                                                        >  1:707204/1‑140 (MQ=255)                                                                 aGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGgcag                                                                      >  2:707149/1‑140 (MQ=255)                                                                 aGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGgcag                                                                      <  1:1182052/140‑1 (MQ=255)                                                                 aGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGgcag                                                                      <  1:846408/140‑1 (MQ=255)                                                                      ccaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGcc                                                                 <  2:446838/140‑1 (MQ=255)                                                                      ccaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGcc                                                                 <  1:705795/140‑1 (MQ=255)                                                                      ccaGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGcc                                                                 <  1:321631/140‑1 (MQ=255)                                                                       caGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCt                                                                <  1:606593/140‑1 (MQ=255)                                                                            gcgcGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTAt                                                           <  2:306493/140‑1 (MQ=255)                                                                                  cTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACa                                                     <  2:87251/140‑1 (MQ=255)                                                                                     tGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCg                                                  >  1:720884/1‑140 (MQ=255)                                                                                       cTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGc                                                <  2:707204/140‑1 (MQ=255)                                                                                       cTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGc                                                <  2:720884/140‑1 (MQ=255)                                                                                              ggCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGc                                                >  1:1039275/1‑133 (MQ=255)                                                                                              ggCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGc                                                <  2:1039275/133‑1 (MQ=255)                                                                                               gCTGGTGTGGATTTTATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGacac                                        >  2:451548/1‑140 (MQ=255)                                                                                               gCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGc                                                >  2:884704/1‑132 (MQ=255)                                                                                               gCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGc                                                <  1:884704/132‑1 (MQ=255)                                                                                               gCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGacac                                        >  2:431477/1‑140 (MQ=255)                                                                                               gCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGa                                           <  2:588156/137‑1 (MQ=255)                                                                                               gCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGa                                           >  1:588156/1‑137 (MQ=255)                                                                                                 tGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACgg                                      >  1:1205305/1‑140 (MQ=255)                                                                                                    tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCg                                                  >  2:411919/1‑125 (MQ=255)                                                                                                    tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCg                                                  <  1:411919/125‑1 (MQ=255)                                                                                                    tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGa                                           <  2:1195163/132‑1 (MQ=255)                                                                                                    tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGa                                           >  1:1195163/1‑132 (MQ=255)                                                                                                    tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCg                                   >  2:186513/1‑140 (MQ=255)                                                                                                    tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCg                                   >  2:973878/1‑140 (MQ=255)                                                                                                    tgtgGATTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCg                                   >  2:833831/1‑140 (MQ=255)                                                                                                         aTTTCATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAg                              >  1:949085/1‑140 (MQ=255)                                                                                                           tttaTTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGc                            <  1:451548/137‑1 (MQ=255)                                                                                                             cATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGt                          >  2:256480/1‑140 (MQ=255)                                                                                                             cATTTACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGt                          >  2:780437/1‑140 (MQ=255)                                                                                                                ttACCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTca                       <  1:189680/140‑1 (MQ=255)                                                                                                                  aCCGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTcacc                     >  1:466204/1‑140 (MQ=255)                                                                                                                   ccGCGCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTcacca                    <  1:186513/140‑1 (MQ=255)                                                                                                                    cgcgCCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTcaccac                   <  1:833831/140‑1 (MQ=255)                                                                                                                     gcgcCAGTAGCGACGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCaccacc                  >  2:770062/1‑140 (MQ=255)                                                                                                                                 aCGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGc         >  2:983256/1‑137 (MQ=255)                                                                                                                                 aCGGATTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGc         <  1:983256/137‑1 (MQ=255)                                                                                                                                     aTTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTg  <  2:949085/140‑1 (MQ=255)                                                                                                                                     aTTGTGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTg  >  2:277982/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        TTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTG  >  NC_000913/2891118‑2891389

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