Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,875,172	T→C	D83G (GAC→GGC)	cysN ←	sulfate adenylyltransferase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,875,172	0	T	C	100.0%	214.0 / NA	60	D83G (GAC→GGC)	cysN	sulfate adenylyltransferase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/29);  total (31/29)

GCCAGTCGCCATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCG  >  NC_000913/2875039‑2875306                                                                                                                                      |                                                                                                                                       gcCAGTCGCCATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCccacca                                                                                                                                  >  1:1112192/1‑140 (MQ=255)  ccAGTCGCCATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACcag                                                                                                                                 <  1:938078/140‑1 (MQ=255)  ccAGTCGCCATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACcag                                                                                                                                 <  1:959127/140‑1 (MQ=255)   cAGTCGCCATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACcagc                                                                                                                                <  2:375788/140‑1 (MQ=255)    aGTCGCCATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACcagca                                                                                                                               <  2:46037/140‑1 (MQ=255)         ccATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGcc                                                                                                                          <  1:1305071/140‑1 (MQ=255)                    gtgtACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCtt                                                                                                               >  1:787353/1‑140 (MQ=255)                     tgtACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCttt                                                                                                              <  2:1347121/140‑1 (MQ=255)                      gtACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCtttt                                                                                                             >  2:424843/1‑140 (MQ=255)                       tACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTc                                                                                                            >  1:929730/1‑140 (MQ=255)                       tACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTc                                                                                                            >  1:541821/1‑140 (MQ=255)                         cTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGc                                                                                                          <  2:929730/140‑1 (MQ=255)                          tGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGcc                                                                                                         <  2:241212/140‑1 (MQ=255)                           gCTCGTGCCCTGGGGGGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGccc                                                                                                        >  2:144539/1‑140 (MQ=255)                              cGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTgg                                                                                                     >  1:219784/1‑140 (MQ=255)                               gTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTggg                                                                                                    >  1:1066839/1‑140 (MQ=255)                                 gCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTg                                                                                                  >  1:955466/1‑140 (MQ=255)                                 gCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTg                                                                                                  >  2:566468/1‑140 (MQ=255)                                 gCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGAAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCTGAGCCCGATCCAGCTTTTCGCCCTGGGTg                                                                                                  >  2:794434/1‑140 (MQ=255)                                     tGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGt                                                                                              <  1:566468/140‑1 (MQ=255)                                     tGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGt                                                                                              <  2:384913/140‑1 (MQ=255)                                      ggggTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTc                                                                                                            >  1:1008700/1‑125 (MQ=255)                                      ggggTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTc                                                                                                            <  2:1008700/125‑1 (MQ=255)                                      ggggTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTg                                                                                             <  2:1090495/140‑1 (MQ=255)                                            tCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCtt                                                                                       >  1:368021/1‑140 (MQ=255)                                             cGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTa                                                                                      >  1:310473/1‑140 (MQ=255)                                             cGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTa                                                                                      <  2:649569/140‑1 (MQ=255)                                              ggCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTAc                                                                                     <  1:53976/140‑1 (MQ=255)                                                   taataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCg                                                                                <  2:955466/140‑1 (MQ=255)                                                   taataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCg                                                                                <  2:362779/140‑1 (MQ=255)                                                       aaaTTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTAt                                                                            <  1:144539/140‑1 (MQ=255)                                                             aCGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCg                                                                      >  1:1252904/1‑140 (MQ=255)                                                             aCGCTTCTCGGAAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCg                                                                      >  2:320446/1‑140 (MQ=255)                                                               gCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGAt                                                                    <  1:424843/140‑1 (MQ=255)                                                                   ctcGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATgaga                                                                >  2:670309/1‑140 (MQ=255)                                                                       gTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCACACCAGCAGAGCCAGATCCAGCTGTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTg                                                            <  2:772053/140‑1 (MQ=255)                                                                             aaaTAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTc                                                      >  1:1048605/1‑140 (MQ=255)                                                                               aTAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGt                                                    >  1:968336/1‑140 (MQ=255)                                                                                tAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTa                                                   >  1:962655/1‑140 (MQ=255)                                                                                  gCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGa                                                 >  2:1164630/1‑140 (MQ=255)                                                                                  gCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGa                                                 >  2:530341/1‑140 (MQ=255)                                                                                  gCGGTAGGCCACGTCAAT‑GTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGAt                                                <  2:472773/140‑1 (MQ=255)                                                                                   cGGTAGGCTACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGAt                                                <  1:320446/140‑1 (MQ=255)                                                                                    ggTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGAtt                                               >  2:190422/1‑140 (MQ=255)                                                                                       aGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTgg                                            <  2:1048605/140‑1 (MQ=255)                                                                                        ggCCACGTCAATGGTGATGCCCTGTTCGTGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGc                                           >  2:243574/1‑140 (MQ=255)                                                                                             cGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGAt                                                <  2:196807/130‑1 (MQ=255)                                                                                             cGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGAt                                                >  1:196807/1‑130 (MQ=255)                                                                                             cGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTa                                      <  2:962655/140‑1 (MQ=255)                                                                                                   tGGTGATGCCCTGTTCGTGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgc                                <  1:243574/140‑1 (MQ=255)                                                                                                   tGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgc                                <  2:682351/140‑1 (MQ=255)                                                                                                       gATGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgcagca                            <  2:1112192/140‑1 (MQ=255)                                                                                                        aTGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTa                                      <  1:868747/129‑1 (MQ=255)                                                                                                        aTGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTa                                      >  2:868747/1‑129 (MQ=255)                                                                                                        aTGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGcagcag                           >  2:1001300/1‑140 (MQ=255)                                                                                                         tGCCCTGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAga                          <  2:787353/140‑1 (MQ=255)                                                                                                              tGTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGAcc                     >  1:191703/1‑140 (MQ=255)                                                                                                               gTTCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGCGTATCGTGCAGCAGACGACCa                    <  2:968336/140‑1 (MQ=255)                                                                                                                ttCGCGCTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCaa                   >  2:1044920/1‑140 (MQ=255)                                                                                                                     gcTCAGCTTGCAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCaga              >  2:1263006/1‑140 (MQ=255)                                                                                                                               cAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGc    >  1:21087/1‑140 (MQ=255)                                                                                                                               cAGGCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGGGCAGCAGACGACCAATCAGAGTACcttttc    >  1:1291384/1‑138 (MQ=255)                                                                                                                                 ggCCGCCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGGGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCg  >  1:1262618/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       GCCAGTCGCCATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCG  >  NC_000913/2875039‑2875306

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: ‑