Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,359,909	C→T	W41* (TGG→TGA)	cadB ←	lysine:cadaverine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,359,909	0	C	T	100.0%	134.2 / NA	43	W41* (TGG→TGA)	cadB	lysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/28);  total (15/28)

CATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCT  >  NC_000913/4359683‑4360142                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          cATTTAAGACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCGCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACt                                                                                                                                                                                                                    <  1:398388/250‑1 (MQ=255)      aaTACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTa                                                                                                                                                                                                               <  1:257384/250‑1 (MQ=255)                    ggTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGtaata                                                                                                                                                                                                <  1:515922/251‑1 (MQ=255)                                     cAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATgtgcattattacctg                                                                                                                                                                               <  1:172018/251‑15 (MQ=255)                                         ggTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATTGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGGTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatca                                                                                                                                                                           <  2:116282/251‑1 (MQ=255)                                          gTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATATAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACATTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatcat                                                                                                                                                                          >  1:59519/1‑251 (MQ=255)                                                         aGGTTACCAATCCAGTTAGCCTGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGcaccaa                                                                                                                                                            >  2:355768/1‑250 (MQ=255)                                                                cAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCTCCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAAGCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTa                                                                                                                                                    <  1:496318/251‑1 (MQ=255)                                                                cAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTa                                                                                                                                                    <  1:264325/251‑1 (MQ=255)                                                                         tAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCaaa                                                                                                                                            <  1:135504/250‑1 (MQ=255)                                                                               ggTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGccc                                                                                                                                      <  1:147980/250‑1 (MQ=255)                                                                               ggTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGAATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACAGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCg                                                                                                                                     >  1:26793/1‑251 (MQ=255)                                                                                  aataaAGAACACCTGTCTGAAAACCAAATGCAGGGGAATGTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAATAACACCGGTACAGGCAAATAGCCCGAtc                                                                                                                                  <  1:220075/251‑1 (MQ=255)                                                                                            acCTGTCCGAAACCCAAATGCAGGGGAAATTTCTCCTGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAg                                                                                                                        <  1:220599/251‑1 (MQ=255)                                                                                             ccTGTCTGAAAACCAAATGCAGGGGAAATTGCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAAAAGAGATAATCCAACCTCAGATAGCAATACCACCTATACTTGCTAGGTTCGCAGGtaata                                                                                                                                                                                                <  1:566649/178‑1 (MQ=255)                                                                                             ccTGTCTGAAAACCAAATGAAGGGGAAATTTATCCGGCATAAGCAGTTGGGAAACCTTGTTGCGGGTTTTTTGATGCCAGTCGGGAATATACATACGCCAGCGCCATTGCACCAATAATAGAGATAATCCAACCTAAGATAGCAATACCACCGCTCCTTGCTCGGTTCGCAGGtaata                                                                                                                                                                                                >  2:566649/1‑178 (MQ=255)                                                                                                gTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAAc                                                                                                                     <  1:105657/250‑1 (MQ=255)                                                                                                         ccAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCtctt                                                                                                           >  1:386292/1‑251 (MQ=255)                                                                                                          cAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCtcttc                                                                                                          >  1:160171/1‑251 (MQ=255)                                                                                                            aaTGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTtctc                                                                                                        >  2:161411/1‑251 (MQ=255)                                                                                                                   gggAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTTTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGACCTCATGCTCTTCTCCTAAttt                                                                                                 <  2:371133/251‑1 (MQ=255)                                                                                                                       aaTTTCTCCGGCATAAGCAATTGGGCCTCCTTGTTGCGGGATTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCAGATCTTCTTGGCAGAACTCATGCTCTTATCCTAATTTCAtt                                                                                             >  2:236617/1‑251 (MQ=255)                                                                                                                              ccGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGaa                                                                                       <  1:199433/250‑1 (MQ=255)                                                                                                                              ccGGCATAAGCAATTGGGCCACCTGGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGCCATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGAGCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTg                                                                                         <  2:420896/248‑1 (MQ=255)                                                                                                                                ggCATAAGCAATTGGTGCACCTTGGTGTGTGTTTTTTGTTGCCAGGCGGGCATATACATGCGCCAGCGACATTGCACCAATAATAGAGAGAATCCAACTTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAAGAATGCAATACCGCTCCCCATCATATTACCGGCATCAACACCTGTACAGGCAAATAGCTCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAAttg                                                                                    <  2:46983/251‑2 (MQ=255)                                                                                                                                    aaaGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAg                                                                                <  1:426698/250‑1 (MQ=255)                                                                                                                                            tGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTa                                                                                                                                                                      >  2:181417/1‑157 (MQ=255)                                                                                                                                            tGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTa                                                                                                                                                                      <  1:181417/157‑1 (MQ=255)                                                                                                                                               gCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATGTCATTTTTGAATTTGGAGTCCGGGTCatg                                                                     <  1:487111/251‑1 (MQ=255)                                                                                                                                                     ttgttgCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGtat                                                               >  1:440145/1‑251 (MQ=255)                                                                                                                                                         tgCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAAc                                                            <  1:569795/250‑1 (MQ=255)                                                                                                                                                            gggTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTAtt                                                        >  2:126644/1‑251 (MQ=255)                                                                                                                                                                     gttgCCAGTCGGGCATATACATACGCCAGCGACAGTGCACCAATAATAGAGATAAACCAACCTCAGATAGCAATACCACCGAGTCATGCTAGGTTCGCAGGTAATAAGGCAATACCTCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGATCTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTgacc                                                <  2:198533/250‑1 (MQ=255)                                                                                                                                                                              cgcgCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACCACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACt                                      >  2:445120/4‑251 (MQ=255)                                                                                                                                                                               gggCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAAGTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTg                                     <  1:334344/251‑1 (MQ=255)                                                                                                                                                                                           acGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGCATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTa                         <  1:161411/251‑1 (MQ=255)                                                                                                                                                                                                  cGACATTGCACCAATAATAGAGATAATGCATCCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAAGTAATAATGCAATACCGCACCCCATCATATGACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCAGACCAGACCAAACGGGCGATAAGATTACTCACGa                  <  2:378371/251‑1 (MQ=255)                                                                                                                                                                                                    aCATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGCTCTTCTTGGCAGAACTCATGCtctt                                                                                                           >  2:550849/1‑160 (MQ=255)                                                                                                                                                                                                    aCATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCtctt                                                                                                           <  1:550849/160‑1 (MQ=255)                                                                                                                                                                                                      aTTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGaaaaa              <  1:157106/251‑1 (MQ=255)                                                                                                                                                                                                               ataataGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAAt     >  1:129497/1‑251 (MQ=255)                                                                                                                                                                                                                 aataGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATcc   >  1:120853/1‑251 (MQ=255)                                                                                                                                                                                                                  ataGAGATAATCCAACCTCAGATAGCAGTACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACTCCGGTACAGGCAAATAGTCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATGTTTGAATTTGGAGTCCGTGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAAGCCt  <  2:389653/251‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          CATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCT  >  NC_000913/4359683‑4360142

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