Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,875,459	C→T	G291R (GGG→AGG)	cysD ←	sulfate adenylyltransferase subunit 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,875,459	0	C	T	100.0%	225.0 / NA	64	G291R (GGG→AGG)	cysD	sulfate adenylyltransferase subunit 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/31);  total (33/31)

GGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGG  >  NC_000913/2875324‑2875592                                                                                                                                        |                                                                                                                                      ggTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGcc                                                                                                                                   >  2:72585/1‑140 (MQ=255) ggTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGcc                                                                                                                                   >  1:665183/1‑140 (MQ=255)  gTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCt                                                                                                                                  >  2:624260/1‑140 (MQ=255)   tCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTg                                                                                                                                 >  2:732989/1‑140 (MQ=255)      gAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTc                                                                                                                              >  2:1037614/1‑140 (MQ=255)            gcagcaGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTc                                                                                                                        >  2:346308/1‑140 (MQ=255)            gcagcaGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTc                                                                                                                        >  1:1493049/1‑140 (MQ=255)                cagGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATc                                                                                                                    >  2:1034129/1‑140 (MQ=255)                cagGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATc                                                                                                                    >  2:290927/1‑140 (MQ=255)                   gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATCCCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  1:423025/140‑1 (MQ=255)                   gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  2:209159/140‑1 (MQ=255)                   gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  1:28491/140‑1 (MQ=255)                   gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  1:755968/140‑1 (MQ=255)                   gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  1:624260/140‑1 (MQ=255)                   gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  1:444498/140‑1 (MQ=255)                   gCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  1:116742/140‑1 (MQ=255)                     ttttATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcgcg                                                                                                               <  1:259350/140‑1 (MQ=255)                      tttATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCgg                                                                                                              >  1:286156/1‑140 (MQ=255)                           gttgttGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTg                                                                                                         <  2:1210555/140‑1 (MQ=255)                             tgttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 >  2:1337090/1‑130 (MQ=255)                             tgttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCAcg                                                                                                                 <  1:1337090/130‑1 (MQ=255)                             tgttgCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGAc                                                                                                       >  2:836289/1‑140 (MQ=255)                                 gcgcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTg                                                                                                         >  2:390780/1‑134 (MQ=255)                                 gcgcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTg                                                                                                         <  1:390780/134‑1 (MQ=255)                                 gcgcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTc                                                                                                   >  2:346611/1‑140 (MQ=255)                                  cgcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCa                                                                                                  <  1:1060630/140‑1 (MQ=255)                                  cgcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCa                                                                                                  <  1:1256378/140‑1 (MQ=255)                                   gcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAGGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCAc                                                                                                 <  1:1083485/140‑1 (MQ=255)                                   gcAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCAc                                                                                                 <  1:1125350/140‑1 (MQ=255)                                      atcatcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACtgg                                                                                              <  2:665183/140‑1 (MQ=255)                                         atcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGa                                                                                                        >  2:1228343/1‑127 (MQ=255)                                         atcCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGa                                                                                                        <  1:1228343/127‑1 (MQ=255)                                              ggCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGaa                                                                                      <  1:499520/140‑1 (MQ=255)                                                 ttCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAAcc                                                                                   >  2:786284/1‑140 (MQ=255)                                                 ttCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAAcc                                                                                   >  1:784970/1‑140 (MQ=255)                                                   cGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAg                                                                                 >  1:242473/1‑140 (MQ=255)                                                    gACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGc                                                                                >  1:776671/1‑140 (MQ=255)                                                                     cGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTc                                                               >  1:1427383/1‑140 (MQ=255)                                                                         ttgttgTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGc                                                           <  2:1427383/140‑1 (MQ=255)                                                                          tgttgtGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCa                                                          <  1:346611/140‑1 (MQ=255)                                                                              gtgCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAgtgt                                                      >  1:165901/1‑140 (MQ=255)                                                                                gCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGttt                                                    <  1:1034129/140‑1 (MQ=255)                                                                                gCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGttt                                                    <  1:290927/140‑1 (MQ=255)                                                                                    gTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGc                                                >  2:1405511/1‑140 (MQ=255)                                                                                     tGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCa                                               <  2:165901/140‑1 (MQ=255)                                                                                            ttCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGa                                          <  2:851896/138‑1 (MQ=255)                                                                                            ttCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGa                                          >  1:851896/1‑138 (MQ=255)                                                                                                           ccccTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCa                         >  1:1365668/1‑140 (MQ=255)                                                                                                              cTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCg                      >  2:1633684/1‑140 (MQ=255)                                                                                                              cTGACGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCg                      >  2:119025/1‑140 (MQ=255)                                                                                                                 aCGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGcc                   >  1:193374/1‑140 (MQ=255)                                                                                                                 aCGTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGcc                   >  1:931159/1‑140 (MQ=255)                                                                                                                   gTTTTTTCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcag                 <  1:1405511/140‑1 (MQ=255)                                                                                                                      ttttCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGt                           <  2:1280924/127‑1 (MQ=255)                                                                                                                      ttttCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGt                           >  1:1280924/1‑127 (MQ=255)                                                                                                                      ttttCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCcagcag              >  1:1080818/1‑140 (MQ=255)                                                                                                                        ttCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc            <  1:346308/140‑1 (MQ=255)                                                                                                                        ttCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc            <  1:1633684/140‑1 (MQ=255)                                                                                                                        ttCAGCTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc            <  1:1565917/140‑1 (MQ=255)                                                                                                                             cTCCATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCg       >  2:65719/1‑140 (MQ=255)                                                                                                                               ccATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTa     <  1:72585/140‑1 (MQ=255)                                                                                                                                cATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc            <  1:821729/132‑1 (MQ=255)                                                                                                                                cATAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc            >  2:821729/1‑132 (MQ=255)                                                                                                                                  tAGACCTCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACgg  >  2:1517710/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                      GGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGG  >  NC_000913/2875324‑2875592

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