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%	143.6 / NA	45	W41* (TGG→TGA)	cadB	lysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/23);  total (22/23)

CCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAA  >  NC_000913/4359668‑4360144                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            ccgccgGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAAGTTCTCCGGCATAAGCAATTGGGCCTCCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCa                                                                                                                                                                                                                                    <  1:663339/251‑1 (MQ=255)        aaCAGGATCATTTAATACTGGGAAGAAGGGGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATaccac                                                                                                                                                                                                                             >  2:308513/1‑251 (MQ=255)              aTCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAAGAGAGATAATCCAACCTCAGATAGCAATACCACCGATAc                                                                                                                                                                                                                       <  2:87503/251‑1 (MQ=255)                 aTTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTg                                                                                                                                                                                                                    >  2:460940/1‑251 (MQ=255)                      aTACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAgg                                                                                                                                                                                                               >  1:223523/1‑251 (MQ=255)                        aCTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGt                                                                                                                                                                                                              <  1:494502/250‑1 (MQ=255)                                    gTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGtaata                                                                                                                                                                                                  >  1:538314/1‑250 (MQ=255)                                              aagataCATCGGTAATACCAATCGCCAGGTTACCAATCCCGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATActc                                                                                                                                                                                       >  2:86371/1‑249 (MQ=255)                                               agataCAGCAGAAATACCAATCGCACGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAACCCAACTGAAGGGGAACTTTCTCCAGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATGGAGCTCATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTACTAATGCAATACCGc                                                                                                                                                                                      >  2:705788/1‑251 (MQ=255)                                                gataCAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCt                                                                                                                                                                                     <  1:601505/251‑1 (MQ=255)                                                             tACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGAAATAATGCAATACCGCTCCCCATCATAtt                                                                                                                                                                         >  1:233742/1‑250 (MQ=255)                                                                         ggTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAacaaca                                                                                                                                                            <  1:379337/251‑1 (MQ=255)                                                                           ttACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACAcc                                                                                                                                                          >  1:339459/1‑251 (MQ=255)                                                                                   ccAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAgg                                                                                                                                                  <  2:648854/251‑1 (MQ=255)                                                                                          gCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAg                                                                                                                                           <  1:309875/251‑1 (MQ=255)                                                                                              ggTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCg                                                                                                                                       <  1:334996/251‑1 (MQ=255)                                                                                                   taaAGAACACCTGTATCAACACCAAATGTAGGGGCAATTTCTACGGCATAAGCAATTGGGCCACCTTGTTGCGGGATTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCATCAACACCGGTACAGGCAAATAGCCCGAtctt                                                                                                                                  >  2:66480/1‑251 (MQ=255)                                                                                                      aGAACACCTGTCTGAACACCAAATGCAGTGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGTTTTTTTGTTGTCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATTCCACCGATACTTGCTAGGTTCGCAGTTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGTCAAATATCCCGATcttctt                                                                                                                               <  2:537736/251‑1 (MQ=255)                                                                                                      aGAACACCGGTCTGAAAACCAAATGTAGGGGAAATTTCGCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTATTGCCAGGCGGGTATATACATACGCCAGCGACATTTCACCAATAATAGAGATACTACAACCTCAGATTGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGTAATACCGCTCCCCATCATATTACCGGCAACAAAACCGGTACAGGCAAATAGCCCGAGcttctt                                                                                                                               <  2:537729/251‑1 (MQ=255)                                                                                                                tCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatca                                                                                                                                                                             >  1:321209/1‑195 (MQ=255)                                                                                                                tCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGAGAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatca                                                                                                                                                                             <  2:321209/195‑1 (MQ=255)                                                                                                                          aaaTGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCtcttct                                                                                                           >  1:692676/1‑251 (MQ=255)                                                                                                                          aaaTGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCtcttct                                                                                                           <  2:81385/251‑1 (MQ=255)                                                                                                                              gtaGGGGAAATTTCTCCGTCATAAGCAATTGAGCAACCTTGTTGCGGGTTTTTTTTTGCCAGTCGGGCATATACATACGCAAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGAAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAAATCCGGTAAAGGCAAATAGCACGATCTTCTTGGGAGAACTCATGCTATTCTcca                                                                                                        >  2:494466/3‑249 (MQ=255)                                                                                                                              gCAGGTGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCt                                                                                                        <  1:153413/250‑1 (MQ=255)                                                                                                                              gCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCt                                                                                                        >  1:354326/1‑250 (MQ=255)                                                                                                                                              cGGCATAAGCAATTGGGCCACCTTGTTGCGTGTTTTTTGTTGCCAGTTGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCATCCTCAGAGAGCAATACCACCTATACTTTCTAGTTTCGCAGGTAATAATGCAAGACCGCTCCCCATCATATTACCGTCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAAtt                                                                                       <  1:736505/251‑1 (MQ=255)                                                                                                                                                                      gttgCGGGTTTTTTGGTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATATAGATAATCCAACCTCAGATAGCAAGACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATGGCATTTTTGAATTTGGAGTCCGGGTCATGATGTATaa                                                               <  2:429596/251‑1 (MQ=255)                                                                                                                                                                               tttttGTTGCCAGTCGAGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTCTTTCCt                                                      >  2:472203/1‑251 (MQ=255)                                                                                                                                                                                  ttgttgCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTgac                                                   >  1:296385/1‑251 (MQ=255)                                                                                                                                                                                             cGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAAACCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAAc                                         <  2:538314/250‑1 (MQ=255)                                                                                                                                                                                              gggCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTg                                       <  1:353112/251‑1 (MQ=255)                                                                                                                                                                                                 cATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTc                                                        >  2:487084/1‑231 (MQ=255)                                                                                                                                                                                                 cATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTc                                                        <  1:487084/231‑1 (MQ=255)                                                                                                                                                                                                             ccAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTc                                                        >  2:442983/1‑219 (MQ=255)                                                                                                                                                                                                             ccAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAAGAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTc                                                        <  1:442983/219‑1 (MQ=255)                                                                                                                                                                                                               aGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCTATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATTTATAACTATTGCCTGACCAGACCAAACTGGCGATAAGATTACTCAc                      <  1:460940/251‑1 (MQ=255)                                                                                                                                                                                                                 cGACATTGCACCAATAATAGAGATAATCCAATCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACGCAcaa                    <  2:570404/251‑3 (MQ=255)                                                                                                                                                                                                                  gACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGAAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGa                    >  2:378067/1‑250 (MQ=255)                                                                                                                                                                                                                          aCCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGAt           >  1:564363/1‑251 (MQ=255)                                                                                                                                                                                                                              ataataGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAAt       >  1:501001/1‑251 (MQ=255)                                                                                                                                                                                                                                  gagagATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCGTCTTGGCAGAACTCATGCTCTGCTCCTAATTTCAGTTTTGAATTTGGAGTCCGGGTCAGGATGTATAACTATTTCCTCACCAGCCTAAACTGGCGATAAGATTTCGCACGAAAAAAGGATTAATCCTa   <  1:705788/250‑1 (MQ=255)                                                                                                                                                                                                                                   agagATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACGATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTa   <  2:481677/250‑1 (MQ=255)                                                                                                                                                                                                                                      gaTAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTaa  >  2:559717/1‑248 (MQ=255)                                                                                                                                                                                                                                       aTAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATCAGATTACTCACGAAAAAAGGATTAATCCTaa  >  1:404052/1‑247 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            CCGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAA  >  NC_000913/4359668‑4360144

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