Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,543,484	A→G	A15A (GCT→GCC)	cysP ←	thiosulfate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,543,484	0	A	G	100.0%	249.0 / NA	69	A15A (GCT→GCC)	cysP	thiosulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/37);  total (32/37)

TTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTTATATACCGTTT  >  NC_000913/2543353‑2543614                                                                                                                                    |                                                                                                                                   ttGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCa                                                                                                                             <  1:268094/139‑1 (MQ=255)   gTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAgc                                                                                                                           >  1:871894/1‑139 (MQ=255)          ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCgagtg                                                                                                                    >  2:649438/1‑139 (MQ=255)              tATCTTTTGTCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGtt                                                                                                                <  2:877736/139‑1 (MQ=255)                     tGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                                                         >  2:6408/1‑139 (MQ=255)                       cccATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAg                                                                                                       >  2:1890/1‑139 (MQ=255)                       cccATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAg                                                                                                       <  1:203500/139‑1 (MQ=255)                         cATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTa                                                                                                     <  1:6408/139‑1 (MQ=255)                                 cTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACg                                                                                             >  2:1044303/1‑139 (MQ=255)                                  tCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACgg                                                                                            <  2:92686/139‑1 (MQ=255)                                    aaaCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGcc                                                                                          <  2:77447/139‑1 (MQ=255)                                      aCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCAt                                                                                        >  1:952301/1‑139 (MQ=255)                                      aCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCAt                                                                                        >  2:127292/1‑139 (MQ=255)                                         gcggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTgc                                                                                     <  2:218502/139‑1 (MQ=255)                                          cggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTgcg                                                                                    <  2:119778/139‑1 (MQ=255)                                          cggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTgcg                                                                                    <  2:369466/139‑1 (MQ=255)                                          cggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTgcg                                                                                    <  2:531135/139‑1 (MQ=255)                                           ggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTgcgc                                                                                   <  1:836744/139‑1 (MQ=255)                                            gATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCa                                                                                  >  2:946227/1‑139 (MQ=255)                                            gATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCa                                                                                  <  2:871894/139‑1 (MQ=255)                                            gATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCa                                                                                  <  2:405180/139‑1 (MQ=255)                                                        cAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAAtt                                                                      >  2:105326/1‑139 (MQ=255)                                                              gCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGa                                                                >  1:611093/1‑139 (MQ=255)                                                               cTCGCGGGATACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGAc                                                               <  1:554841/139‑1 (MQ=255)                                                               cTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTgc                                                                                     <  2:505827/117‑1 (MQ=255)                                                               cTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTgc                                                                                     >  1:505827/1‑117 (MQ=255)                                                                tCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACt                                                              <  1:564804/139‑1 (MQ=255)                                                                tCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACt                                                              <  2:321310/139‑1 (MQ=255)                                                                           gTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGcc                                                   <  2:821444/139‑1 (MQ=255)                                                                              aTAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCcatc                                                <  1:485916/139‑1 (MQ=255)                                                                                     cTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAAc                                                                                              >  1:223442/1‑86 (MQ=255)                                                                                     cTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAAc                                                                                              <  2:223442/86‑1 (MQ=255)                                                                                      tGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTtat                                        >  1:400596/1‑139 (MQ=255)                                                                                        ttCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGTCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATaa                                      >  2:583371/1‑139 (MQ=255)                                                                                        ttCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTtata                                                                          >  2:384125/1‑103 (MQ=255)                                                                                        ttCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTtata                                                                          <  1:384125/103‑1 (MQ=255)                                                                                        ttCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATaa                                      <  1:1890/139‑1 (MQ=255)                                                                                          cagcagTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACg                                    <  1:605532/139‑1 (MQ=255)                                                                                            gcagTTCCGTTGCCTGTGCATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTA‑‑‑‑‑‑TAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGac                            <  1:649438/139‑2 (MQ=255)                                                                                                  ccGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACt                                                              >  1:323949/1‑105 (MQ=255)                                                                                                  ccGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACt                                                              <  2:323949/105‑1 (MQ=255)                                                                                                  ccGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGaa                            >  2:514276/1‑139 (MQ=255)                                                                                                     ttGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACg                                    <  1:161879/128‑1 (MQ=255)                                                                                                     ttGCCTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACg                                    >  2:161879/1‑128 (MQ=255)                                                                                                        ccTGTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGt                      <  2:735324/139‑1 (MQ=255)                                                                                                           gTACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAAc                   >  1:991747/1‑139 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATaa                                      <  2:579157/119‑1 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATaa                                      >  1:579157/1‑119 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  >  2:609262/1‑139 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  <  2:991747/139‑1 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  <  1:1044303/139‑1 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  >  2:408498/1‑139 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  <  1:795560/139‑1 (MQ=255)                                                                                                            tACATGGCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  <  1:408498/139‑1 (MQ=255)                                                                                                                  gCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAAt                                                       <  2:488950/96‑1 (MQ=255)                                                                                                                  gCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAAt                                                       >  1:488950/1‑96 (MQ=255)                                                                                                                  gCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTta            >  1:867176/1‑139 (MQ=255)                                                                                                                  gCCCGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTta            >  1:774705/1‑139 (MQ=255)                                                                                                                   cccGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGaa                            <  2:691981/122‑1 (MQ=255)                                                                                                                   cccGCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGaa                            >  1:691981/1‑122 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTaa                                                        >  2:487225/1‑91 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTaa                                                        <  1:487225/91‑1 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTaa                    <  1:859279/127‑1 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTaa                    >  2:859279/1‑127 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  >  1:14559/1‑129 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACg                  <  2:14559/129‑1 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTtatata        >  2:937920/1‑139 (MQ=255)                                                                                                                      gCCAGCAGCAGAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGCAAGTTAACGTCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTtatata        >  1:74719/1‑139 (MQ=255)                                                                                                                            agcagAGAGGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATNNNTAACGGACGTGAAAGGAGTAACGGTTTTATATACCGttt  >  1:563222/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                   TTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAAAGGAGTAACGGTTTTATATACCGTTT  >  NC_000913/2543353‑2543614

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