Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,553,284	T→C	C20C (TGT→TGC)	hemF →	coproporphyrinogen III oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,553,284	0	T	C	100.0%	248.1 / NA	66	C20C (TGT→TGC)	hemF	coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/33);  total (33/33)

ATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTC  >  NC_000913/2553153‑2553419                                                                                                                                    |                                                                                                                                        aTTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCt                                                                                                                                 <  1:1381753/140‑1 (MQ=255)   tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGa                                                                                                                               <  2:220971/140‑1 (MQ=255)    gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGAc                                                                                                                              >  2:297964/1‑140 (MQ=255)          ggCGAGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGAccgccgg                                                                                                                        >  1:1349448/1‑139 (MQ=255)            cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCg                                                                                                                      <  2:437620/140‑1 (MQ=255)             gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGa                                                                                                                     <  1:1183320/140‑1 (MQ=255)             gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGa                                                                                                                     <  1:1244837/140‑1 (MQ=255)             gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGa                                                                                                                     <  1:595946/140‑1 (MQ=255)               gATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATg                                                                                                                   >  1:727119/1‑140 (MQ=255)                     ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAg                                                                                                             >  1:1564159/1‑140 (MQ=255)                        tttCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAt                                                                                                          <  2:343702/140‑1 (MQ=255)                        tttCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAt                                                                                                          <  2:1435451/140‑1 (MQ=255)                        tttCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAt                                                                                                          <  2:617014/140‑1 (MQ=255)                        tttCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAt                                                                                                          <  1:314691/140‑1 (MQ=255)                          tCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAttt                                                                                                        >  1:1504875/1‑140 (MQ=255)                          tCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAttt                                                                                                        >  1:912650/1‑140 (MQ=255)                          tCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAttt                                                                                                        >  1:180907/1‑140 (MQ=255)                          tCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAAttt                                                                                                        >  2:1444587/1‑140 (MQ=255)                            cACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGt                                                                                                      <  1:42389/140‑1 (MQ=255)                               tGGTTCGACAACCAGAAAGCACATTCGAAAAAGTGTTAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGa                                                                                                   <  2:870177/140‑1 (MQ=255)                               tGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGa                                                                                                   >  1:1517809/1‑140 (MQ=255)                               tGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGa                                                                                                   <  1:25289/140‑1 (MQ=255)                               tGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGa                                                                                                   <  1:836384/140‑1 (MQ=255)                               tGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGa                                                                                                   <  1:1290240/140‑1 (MQ=255)                                      aCAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGt                                                                                            >  1:1520848/1‑140 (MQ=255)                                      aCAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGt                                                                                            >  2:1456043/1‑140 (MQ=255)                                           cAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCa                                                                                       >  1:1009197/1‑140 (MQ=255)                                           cAGAAAGCACATTCGAAAAAACGGTTAGTTATGAAACCCGACGCACACCTGGTTAAACAGCTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCa                                                                                       >  2:1279105/1‑140 (MQ=255)                                                     aTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCt                                                                             >  1:1492434/1‑140 (MQ=255)                                                     aTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCt                                                                             >  2:59450/1‑140 (MQ=255)                                                       tCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTgg                                                                           >  1:773119/1‑140 (MQ=255)                                                            aaaGCGTTAAGTTATGAAACCCGACGCAGACTAGGTTAAGCAGTTTCTGCCCAACCCTCAGGATACGAATTGCCAGCAGCTGACCGCCGTCGATGGTGCAGCATTTGTCGAAGATAGTTGGCAGCGCGAAGCTggaggcg                                                                      <  2:1161836/140‑1 (MQ=255)                                                            aaaGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTggcggcg                                                                      <  1:804129/140‑1 (MQ=255)                                                            aaaGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTggcggcg                                                                      <  2:1564159/140‑1 (MQ=255)                                                             aaGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTggcggcgg                                                                     <  1:647692/140‑1 (MQ=255)                                                               gCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGAT‑GCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGcggcggcgg                                                                  <  1:1117051/140‑1 (MQ=255)                                                                   tAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCggcgggcg                                                               >  1:611344/1‑140 (MQ=255)                                                                       ttATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGt                                                           <  1:1249462/140‑1 (MQ=255)                                                                          tGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCgg                                                        >  1:919698/1‑140 (MQ=255)                                                                          tGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCgg                                                        >  2:561018/1‑140 (MQ=255)                                                                          tGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCgg                                                        <  2:1213994/140‑1 (MQ=255)                                                                          tGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCgg                                                        <  2:1052316/140‑1 (MQ=255)                                                                          tGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCgg                                                        <  1:1624606/140‑1 (MQ=255)                                                                           gAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCggg                                                       <  1:930361/140‑1 (MQ=255)                                                                            aaaCCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCgggg                                                      >  2:176958/1‑139 (MQ=255)                                                                             aaCCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGgtg                                                     >  2:1403396/1‑140 (MQ=255)                                                                             aaCCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGgtg                                                     <  2:840167/140‑1 (MQ=255)                                                                                 cGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGc                                                 >  1:630068/1‑140 (MQ=255)                                                                                  gACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCg                                                <  1:1359135/140‑1 (MQ=255)                                                                                  gACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCg                                                <  2:414323/140‑1 (MQ=255)                                                                                    cGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTa                                              >  1:1130295/1‑140 (MQ=255)                                                                                          ccAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAAtggtg                                        >  1:845386/1‑140 (MQ=255)                                                                                             ggTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGgtggtg                                     >  2:633057/1‑140 (MQ=255)                                                                                             ggTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGgtggtg                                     >  2:465402/1‑140 (MQ=255)                                                                                              gTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGgtgt                                    >  1:1089407/1‑140 (MQ=255)                                                                                                   aCAGTTTCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCg                               >  2:1429321/1‑140 (MQ=255)                                                                                                       tttCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAAca                           >  2:123828/1‑140 (MQ=255)                                                                                                         tCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAAcagg                         <  2:727119/140‑1 (MQ=255)                                                                                                         tCTGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAAcagg                         <  2:919698/140‑1 (MQ=255)                                                                                                           tGCTCAACCTTCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAAcaggca                       <  2:845386/140‑1 (MQ=255)                                                                                                                     tCAGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGGAATGGGGGGGGTTTCGAACAGGCAGGCGTCAACt             >  2:26457/1‑140 (MQ=255)                                                                                                                       aGGATACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACttt           >  1:1159155/1‑140 (MQ=255)                                                                                                                           tACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGc       <  2:1492434/140‑1 (MQ=255)                                                                                                                           tACGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGc       <  2:538929/140‑1 (MQ=255)                                                                                                                            aCGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCa      >  2:1026131/1‑140 (MQ=255)                                                                                                                            aCGATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCa      >  1:842520/1‑140 (MQ=255)                                                                                                                              gATTTGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATg    >  2:149351/1‑140 (MQ=255)                                                                                                                               aTTTGCCAGCGGCTGACCGACGTCGATGGCGCAGAAGTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGGGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGt   <  1:549585/140‑1 (MQ=255)                                                                                                                                tttGCCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTc  <  2:630068/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        ATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGGCGCAGAATTTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTC  >  NC_000913/2553153‑2553419

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