Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,889,238	C→A	A288S (GCC→TCC)	cysI ←	sulfite reductase, hemoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,889,238	0	C	A	100.0%	244.2 / NA	72	A288S (GCC→TCC)	cysI	sulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (34/38);  total (34/38)

GTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTT  >  NC_000913/2889103‑2889373                                                                                                                                        |                                                                                                                                        gTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGACATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCg                                                                                                                                     <  1:836687/140‑1 (MQ=255)   tCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGc                                                                                                                                   <  1:1600612/140‑1 (MQ=255)         ttGATCCCCGCGCGACGCTCCCCTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGc                                                                                                                             <  2:959679/140‑1 (MQ=255)          tGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGcc                                                                                                                            <  2:769287/140‑1 (MQ=255)          tGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGcc                                                                                                                            <  1:869888/140‑1 (MQ=255)           gATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCa                                                                                                                           <  1:1030799/140‑1 (MQ=255)           gATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCa                                                                                                                           <  2:1661687/140‑1 (MQ=255)                  gcgcgACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTTCCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATg                                                                                                                    <  2:893433/140‑1 (MQ=255)                    gcgACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCt                                                                                                                  <  2:753098/140‑1 (MQ=255)                      gACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTcc                                                                                                                >  1:326704/1‑140 (MQ=255)                      gACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTcc                                                                                                                >  1:463483/1‑140 (MQ=255)                               cTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTGCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGa                                                                                                       <  1:288265/140‑1 (MQ=255)                                  ccGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAg                                                                                                    >  1:225681/1‑140 (MQ=255)                                   cGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGc                                                                                                   <  1:665884/140‑1 (MQ=255)                                       ttAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCaaa                                                                                               >  2:533570/1‑140 (MQ=255)                                       ttAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCaaa                                                                                               <  2:356739/140‑1 (MQ=255)                                         aaaCGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAact                                                                                             <  2:463483/140‑1 (MQ=255)                                           aCGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCACTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAactca                                                                                           <  1:1343864/140‑1 (MQ=255)                                               ctcACCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTCCCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcg                                                                                       <  2:651322/140‑1 (MQ=255)                                               ctcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcg                                                                                       <  2:426878/140‑1 (MQ=255)                                               ctcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcg                                                                                       <  2:320421/140‑1 (MQ=255)                                               ctcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcg                                                                                       <  1:308412/140‑1 (MQ=255)                                                tcAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTcgc                                                                                      >  1:1486230/1‑140 (MQ=255)                                                   aCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGt                                                                                   <  1:1151778/140‑1 (MQ=255)                                                    cccccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTg                                                                                  <  1:405805/140‑1 (MQ=255)                                                    cccccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTg                                                                                  <  2:1654241/140‑1 (MQ=255)                                                    cccccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTg                                                                                  <  2:1486230/140‑1 (MQ=255)                                                       ccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCgg                                                                               >  2:174340/1‑140 (MQ=255)                                                        cACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCggg                                                                              >  1:1307779/1‑140 (MQ=255)                                                         aCGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGc                                                                             >  1:1294487/1‑140 (MQ=255)                                                         aCGCGCTCCAGCGTGTATTTGGTTTTGGCATTTCTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGc                                                                             >  2:1369930/1‑140 (MQ=255)                                                           gcgcTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGt                                                                           <  1:1046911/140‑1 (MQ=255)                                                              cTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAgg                                                                        >  2:1402154/1‑140 (MQ=255)                                                                 cAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGttt                                                                     >  1:327210/1‑140 (MQ=255)                                                                 cAGCGTGTATTTGGTTTTGGCATTTCTACGATCGGTTCGGTTCCCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGttt                                                                     <  1:1369930/140‑1 (MQ=255)                                                                      tgtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCttg                                                                >  2:1125400/1‑140 (MQ=255)                                                                      tgtATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCttg                                                                <  1:533570/140‑1 (MQ=255)                                                                         aTTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCttgttg                                                             >  2:887822/1‑140 (MQ=255)                                                                         aTTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCttgttg                                                             >  2:1174857/1‑140 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCttg                                                                <  2:390472/128‑1 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCttg                                                                >  1:390472/1‑128 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCg                                                    >  2:1298406/1‑140 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCg                                                    >  1:760212/1‑140 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCg                                                    >  2:772829/1‑140 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCg                                                    >  1:179459/1‑140 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCg                                                    >  1:1383318/1‑140 (MQ=255)                                                                                  tGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCg                                                    <  1:1402154/140‑1 (MQ=255)                                                                                        ttttACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGaa                                              >  1:1324580/1‑140 (MQ=255)                                                                                            aCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAgcc                                          <  2:225681/140‑1 (MQ=255)                                                                                                 cggttcggttACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAgcccaccg                                     <  2:1294487/140‑1 (MQ=255)                                                                                                    ttcggttACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAgcccaccgccc                                  <  2:1383318/140‑1 (MQ=255)                                                                                                    ttcggttACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAgcccaccgccc                                  >  1:249998/1‑140 (MQ=255)                                                                                                         ttACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCaa                             <  2:760212/140‑1 (MQ=255)                                                                                                          tACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAAc                            <  1:1125400/140‑1 (MQ=255)                                                                                                            ccccAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAgcccaccgcc                                   <  1:1618533/131‑1 (MQ=255)                                                                                                            ccccAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAgcccaccgcc                                   >  2:1618533/1‑131 (MQ=255)                                                                                                              ccAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGt                        >  2:655974/1‑140 (MQ=255)                                                                                                                 gTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTaa                     >  2:790500/1‑140 (MQ=255)                                                                                                                 gTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTaa                     <  1:790500/140‑1 (MQ=255)                                                                                                                  tCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTaaa                    <  1:1298406/140‑1 (MQ=255)                                                                                                                  tCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTaaa                    >  2:321022/1‑140 (MQ=255)                                                                                                                  tCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTaaa                    >  2:1354679/1‑140 (MQ=255)                                                                                                                   cACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAg                   <  1:1174857/140‑1 (MQ=255)                                                                                                                   cACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCAGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAg                   >  2:425730/1‑140 (MQ=255)                                                                                                                       cTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCcca               >  1:614790/1‑140 (MQ=255)                                                                                                                        tGAGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCccac              <  1:416233/140‑1 (MQ=255)                                                                                                                          aGTTGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCccacca            >  1:773112/1‑140 (MQ=255)                                                                                                                            ttGTCACGACGGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGc          >  2:108827/1‑140 (MQ=255)                                                                                                                                 acgacgGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGcc     >  1:4658/1‑140 (MQ=255)                                                                                                                                    acgGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCccacca            <  1:147053/130‑1 (MQ=255)                                                                                                                                    acgGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCccacca            >  2:147053/1‑130 (MQ=255)                                                                                                                                    acgGATTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGtt  >  1:1311292/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        GTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTT  >  NC_000913/2889103‑2889373

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