Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	388,700	(G)5→4	intergenic (+54/+53)	tauD → / ← hemB	taurine dioxygenase/porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	388,696	0	G	.	100.0%	287.1 / NA	66	intergenic (+50/+57)	tauD/hemB	taurine dioxygenase/porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (42/24);  total (42/24)

ACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTA  >  NC_000913/388558‑388825                                                                                                                                           |                                                                                                                                  aCTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                  <  2:1031970/140‑1 (MQ=255)     tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggcg                                                                                                                              >  2:55529/1‑138 (MQ=255)     tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggcg                                                                                                                              >  1:671938/1‑138 (MQ=255)     tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggcg                                                                                                                              >  1:689690/1‑138 (MQ=255)     tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGCGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggcg                                                                                                                              >  1:986874/1‑138 (MQ=255)     tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGGAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAggggcg                                                                                                                              >  2:793118/1‑138 (MQ=255)        cAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTc                                                                                                                          <  2:766838/140‑1 (MQ=255)          aTGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCAc                                                                                                                        >  1:153781/1‑140 (MQ=255)           tGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAg                                                                                                                                   <  1:624913/129‑1 (MQ=255)           tGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAg                                                                                                                                   >  2:624913/1‑129 (MQ=255)           tGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCAcc                                                                                                                       <  2:1197237/140‑1 (MQ=255)             ccGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGt                                                                                                                     >  2:682508/1‑140 (MQ=255)               gATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTAc                                                                                                                   >  1:1119383/1‑140 (MQ=255)                     cTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGGGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCaa                                                                                                             >  2:1133546/1‑140 (MQ=255)                       gCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACa                                                                                                           >  1:1222724/1‑140 (MQ=255)                        ccACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCaa                                                                                                             >  1:357898/1‑137 (MQ=255)                        ccACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCaa                                                                                                             <  2:357898/137‑1 (MQ=255)                             gCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTaa                                                                                                     <  2:94009/140‑1 (MQ=255)                               gACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACt                                                                                                   <  2:671938/140‑1 (MQ=255)                                 cGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTcc                                                                                                 <  2:1119383/140‑1 (MQ=255)                                      aaTGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCtt                                                                                            >  2:62653/1‑140 (MQ=255)                                           aTCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTga                                                                                       >  2:1384482/1‑140 (MQ=255)                                            tCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTgag                                                                                      <  2:1091364/140‑1 (MQ=255)                                             cGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTgaga                                                                                     >  2:1147013/1‑140 (MQ=255)                                             cGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTgaga                                                                                     >  1:853252/1‑140 (MQ=255)                                              gggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGagag                                                                                    <  2:1387759/140‑1 (MQ=255)                                                   acgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAAc                                                                               >  2:1221686/1‑140 (MQ=255)                                                            ggggATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAAc                                                                      <  1:1133546/140‑1 (MQ=255)                                                             gggATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACg                                                                     <  1:1147013/140‑1 (MQ=255)                                                               gATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCa                                                                   <  1:441376/140‑1 (MQ=255)                                                                 tAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGa                                                                 >  2:325496/1‑140 (MQ=255)                                                                    aCCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtc                                                              >  2:1168155/1‑140 (MQ=255)                                                                          ttATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGG‑TAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttct                                                       <  1:55529/140‑1 (MQ=255)                                                                           tATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttct                                                       <  2:153781/140‑1 (MQ=255)                                                                           tATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAAtcttcttct                                                       <  2:853252/140‑1 (MQ=255)                                                                                ggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTGACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGcc                                                  >  2:682043/1‑140 (MQ=255)                                                                                ggcgggGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGcc                                                  >  1:1292385/1‑140 (MQ=255)                                                                                     ggTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATc                                             <  1:62653/140‑1 (MQ=255)                                                                                              agagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCaaa                                    >  1:1040039/1‑140 (MQ=255)                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGATGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                   >  1:1245525/1‑140 (MQ=255)                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                   >  1:1007958/1‑140 (MQ=255)                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                   >  1:1105215/1‑140 (MQ=255)                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                   >  1:234926/1‑140 (MQ=255)                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACCCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                   >  1:94829/1‑140 (MQ=255)                                                                                                agTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGt                                  <  1:1168155/140‑1 (MQ=255)                                                                                                agTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGt                                  >  1:1276453/1‑140 (MQ=255)                                                                                                 gTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTa                                 >  2:1153904/1‑140 (MQ=255)                                                                                                 gTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTa                                 >  1:808024/1‑140 (MQ=255)                                                                                                  tGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTGCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAg                                >  2:988123/1‑140 (MQ=255)                                                                                                    gACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCt                              >  1:1178905/1‑140 (MQ=255)                                                                                                    gACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCt                              >  1:19551/1‑140 (MQ=255)                                                                                                         tCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAg                                          <  2:737178/123‑1 (MQ=255)                                                                                                         tCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAg                                          >  1:737178/1‑123 (MQ=255)                                                                                                         tCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGaaaa                         >  1:7086/1‑140 (MQ=255)                                                                                                         tCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGaaaa                         >  1:586677/1‑140 (MQ=255)                                                                                                          ccccTCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAAt                        <  2:19551/140‑1 (MQ=255)                                                                                                              tCGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGa                    >  1:750244/1‑140 (MQ=255)                                                                                                               cGCCCCCTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGAt                   <  2:1222724/140‑1 (MQ=255)                                                                                                                  ccccTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccg                <  2:586677/140‑1 (MQ=255)                                                                                                                   cccTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTg                             >  2:968358/1‑126 (MQ=255)                                                                                                                   cccTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTg                             <  1:968358/126‑1 (MQ=255)                                                                                                                   cccTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgc               >  2:602465/1‑140 (MQ=255)                                                                                                                    ccTTGGGGAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgca              <  1:1221686/140‑1 (MQ=255)                                                                                                                        ggggAGAGGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccg          <  2:689690/140‑1 (MQ=255)                                                                                                                            agagGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACg      <  2:1292385/140‑1 (MQ=255)                                                                                                                              agGGTTAGTGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCt    >  1:408362/1‑140 (MQ=255)                                                                                                                              agGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCt    >  2:834201/1‑140 (MQ=255)                                                                                                                              agGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCt    >  2:1064999/1‑140 (MQ=255)                                                                                                                              agGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCt    >  1:691317/1‑140 (MQ=255)                                                                                                                              agGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCt    >  1:1287995/1‑140 (MQ=255)                                                                                                                              agGGTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCt    >  2:155882/1‑140 (MQ=255)                                                                                                                                ggTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTa  >  2:671783/1‑140 (MQ=255)                                                                                                                                ggTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTa  >  2:1190578/1‑140 (MQ=255)                                                                                                                                ggTTAGGGTGA‑GGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTa  >  2:182673/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                  ACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTA  >  NC_000913/388558‑388825

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