Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,129,619	C→T	G316G (GGC→GGT)	metB →	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,129,619	0	C	T	100.0%	179.5 / NA	53	G316G (GGC→GGT)	metB	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/32);  total (21/32)

TGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGT  >  NC_000913/4129483‑4129754                                                                                                                                         |                                                                                                                                         tGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCtgggtggg                                                                                                                                       <  1:495969/140‑1 (MQ=255) tGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCtgggtggg                                                                                                                                       <  1:811333/140‑1 (MQ=255) tGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCtgggtggg                                                                                                                                       <  1:1549221/140‑1 (MQ=255) tGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCtgggtggg                                                                                                                                       <  1:1282234/140‑1 (MQ=255)  ggTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggtgggc                                                                                                                                      <  1:593067/140‑1 (MQ=255)  ggTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggtgggc                                                                                                                                      <  1:634174/140‑1 (MQ=255)  ggTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggtgggc                                                                                                                                      <  1:449779/140‑1 (MQ=255)   gTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCt                                                                                                                                     >  1:659414/1‑140 (MQ=255)     gAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGt                                                                                                                                   <  2:307928/140‑1 (MQ=255)       aaaaaCTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCg                                                                                                                                 <  1:1400121/140‑1 (MQ=255)                aTCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTc                                                                                                                                  >  1:663854/1‑130 (MQ=255)                aTCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTc                                                                                                                                  <  2:663854/130‑1 (MQ=255)                   ccccGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTg                                                                                                                     <  2:1256719/139‑1 (MQ=255)                    cccGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACg                                                                                                                        >  1:1357565/1‑136 (MQ=255)                    cccGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACg                                                                                                                        <  2:1357565/136‑1 (MQ=255)                    cccGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGAGTGGGCTGTCGTTGTTTACGCTgg                                                                                                                    >  2:206636/1‑140 (MQ=255)                      cgtcgtTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCg                                                                                                                  <  2:659414/140‑1 (MQ=255)                       gtcgtTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCgg                                                                                                                 <  1:525983/140‑1 (MQ=255)                         cgtTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGaa                                                                                                               >  1:463913/1‑140 (MQ=255)                             gCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCAt                                                                                                           >  2:962460/1‑140 (MQ=255)                               cGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTa                                                                                                         >  1:650728/1‑140 (MQ=255)                               cGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTa                                                                                                         >  2:1456411/1‑140 (MQ=255)                               cGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTa                                                                                                         >  1:471993/1‑140 (MQ=255)                                   aaaTCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑gggg                                                                                                    <  2:1440278/140‑1 (MQ=255)                                     aTCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGa                                                                                                  <  1:465920/140‑1 (MQ=255)                                     aTCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGa                                                                                                  <  1:1227320/140‑1 (MQ=255)                                         ggggCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTgg                                                                                              <  1:1209440/140‑1 (MQ=255)                                              aTGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGt                                                                                         >  1:1186697/1‑140 (MQ=255)                                                    ttGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATAAGCAGCCGCGGCGTCGTTTCCTGGGTGGCCTGTCGTTGTTTACGCGGGCGGAATCATTA‑GGGGGGGGGGAAAGTTTAAtc                                                                                   >  2:258164/1‑140 (MQ=255)                                                       ccGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAAtctct                                                                                <  1:1244311/140‑1 (MQ=255)                                                         gcgcgcCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGGGTGGAAAGTTTAATCTCTCa                                                                              <  2:887585/140‑1 (MQ=255)                                                           gcgcCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCAcg                                                                            >  1:260483/1‑140 (MQ=255)                                                              ccAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCAcgccg                                                                         <  1:273213/140‑1 (MQ=255)                                                               cAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCAcgccgc                                                                        <  1:216756/140‑1 (MQ=255)                                                                   aaaaaGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAAcc                                                                    <  2:650728/140‑1 (MQ=255)                                                                       aGGTTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAACCTTTCACGCCGCAACCATGa                                                                <  2:1143773/140‑1 (MQ=255)                                                                       aGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGa                                                                <  1:790837/140‑1 (MQ=255)                                                                       aGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGa                                                                <  1:884451/140‑1 (MQ=255)                                                                         gCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGaca                                                              <  1:10478/140‑1 (MQ=255)                                                                            ttGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGa                                                                <  2:431397/135‑1 (MQ=255)                                                                            ttGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGa                                                                >  1:431397/1‑135 (MQ=255)                                                                             tGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGGTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCAcgccg                                                                         <  2:1562614/125‑1 (MQ=255)                                                                             tGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCAcgccg                                                                         >  1:1562614/1‑125 (MQ=255)                                                                             tGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATg                                                          >  1:1032721/1‑140 (MQ=255)                                                                                    aTGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCAt                                                   <  2:1186697/140‑1 (MQ=255)                                                                                    aTGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCAt                                                   <  1:1456411/140‑1 (MQ=255)                                                                                               tGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAg                                        >  1:251742/1‑140 (MQ=255)                                                                                                 aaCTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAgcg                                      <  2:463913/140‑1 (MQ=255)                                                                                                      gATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGtgc                                 >  1:55585/1‑140 (MQ=255)                                                                                                         ggCGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGtgctgc                              >  2:465059/1‑140 (MQ=255)                                                                                                           cGATGAGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTgccg                            >  1:1027014/1‑140 (MQ=255)                                                                                                                aGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGtg                                  >  1:492325/1‑129 (MQ=255)                                                                                                                aGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGtg                                  <  2:492325/129‑1 (MQ=255)                                                                                                                aGCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCggg                       <  2:471993/140‑1 (MQ=255)                                                                                                                 gCAGACGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCggg                       <  2:3655/140‑1 (MQ=255)                                                                                                                     aCGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTgccgcc                          >  2:96735/1‑132 (MQ=255)                                                                                                                     aCGCTGCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTgccgcc                          <  1:96735/132‑1 (MQ=255)                                                                                                                          gCGTCGTTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCgaga             <  2:733881/140‑1 (MQ=255)                                                                                                                            gtcgtTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACg           <  1:324754/140‑1 (MQ=255)                                                                                                                            gtcgtTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACg           >  2:563716/1‑140 (MQ=255)                                                                                                                             tcgtTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACgc          <  2:260483/140‑1 (MQ=255)                                                                                                                               gtTTCCTGGGTGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACgctg        <  1:465059/140‑1 (MQ=255)                                                                                                                                     tgggtgggCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGt  <  1:258164/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         TGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTA‑GGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAACCATGACACATGCAGGCATGGCACCAGAAGCGCGTGCTGCCGCCGGGATCTCCGAGACGCTGCTGCGT  >  NC_000913/4129483‑4129754

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